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Preliminary Remarks 


Water diffusing into silica surfaces gives rise for several effects on diffusion behaviour 


and mechanical properties. The most interesting observations from our point of view may 
be listed here. 


(D Zn the absence of externally applied stresses: 


a) 
b) 
c) 
d) 
e) 
f) 
8) 


Density reduction with increasing water concentration [1]. 
Increasing concentration of molecular water at silica surfaces [2]. 
Decrease of measured water diffusivity with time [3]. 

Change of the shape of the water concentration with time [3]. 
Deformation by one-side water vapour soaked specimens [4]. 
Strength increase of water vapour soaked specimens. 


Increasing exponent n of the power law for subcritical crack growth with in- 


creasing soaking time. 


These points were covered in Part I. 


(II) Jn presence of externally applied loadings: 


a) 


b) 


c) 
d) 


e) 


Increasing diffusivity under externally applied tensile stresses and a decrease 


in compression [5]. 
Decrease of water solubility in tension and an increase in compression [5]. 
Increased water concentration in tension at high temperatures [6]. 


Strongly increasing strength at high tensile stresses even under very low water 


vapour pressure [7]. 
Residual deformation of thin fibers after 2-point bending tests [8]. 


Drastically increased water concentration on crack surfaces after subcritical 


crack growth experiments [9]. 


These topics are elaborated on in Part II. 


Preliminary Remarks 


(ID  Zncluding damage by hydroxyl generation: 


a) 


8) 


Damage behavior and Young’ modulus reduction 


Blunting of crack tip displacements in water-soaked silica [10]. 


Part III will deal with this topic. 


For an explanation of all the findings listed under (I) there is only volume swelling by 


hydrogen generation an appropriate tool. Any stress relaxation mechanism as suggested 


in [7] will hardly be able to explain disk deformation in stress-free disks. 


In order to give a transparent description, two main assumptions will be made: 


(1) 


(2) 


For an interpretation of the increase of diffusivity and water concentration under 
tensile loading, a further effect may play a role. Under tensile stresses, the nano- 
pores of the ring structure are possibly widened and diffusion and solubility may 
be enhanced. Such behaviour was observed in different materials and is interpret- 
ed in terms of the so-called free volume theory. This may also affect water diffu- 
sivity, generally described by the Doolittle equation [11], and solubility (see e.g. 
[12]). Due to the lack of results for silica, free volume effects must be ignored in 


the following considerations. 


Water molecules are believed to diffuse through silica glass in much the same way 
as the noble gases. Hence, the diffusion of water molecules is also expected to 
show little dependence on applied stress [9]. Following this reasoning, we assume 
that the molecular water in the glass is located in molecular size holes in the silica 
glass network where the water does not contribute to the volume of the glass. 
Consequently, molecular water cannot contribute stress effects. Thus, the concen- 
tration of molecular water in the glass depends only on the vapour pressure in the 


external environment, but not on the stress in the glass. 
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XI 


1 Basic results on diffusion and 
swelling 


In this Introductory Section, some of the results from Part I [1.1] may be reported which 
will be used in the following considerations. The relationships from Part I, required in 
this booklet, should be briefly reproduced. Details on derivations and applications can be 


seen in [1.1]. 


Water penetrated into silica reacts with the silica network according to 


=Si-O-Si= +H20 e =SiOH+HOSi= (1.1.1) 


with the concentration of the hydroxyl S=[=SiOH] and that of the molecular water C= 


[H20]. The equilibrium constant X of this reaction is depending on the temperature 0: 


S/C for @< 450°C 


= (1.1.2) 
S^/C for 92450°C 
In the range of 90?Cx0x300*C the equilibrium constant can be expressed by [1.2] 
sa let (1.1.3) 
RT 


where 7=0+273K is the absolute temperature, and R the gas constant. 


Due to diffusion, concentrations of molecular water, C, and hydroxyl water, S, decrease 
with increasing distance, z, from the glass surface. If the reaction given by eq.(1.1.1) is in 


equilibrium, the diffusion process is governed by the diffusion differential equation: 


OC © po] (1.1.4) 
Oz 


ôt oz 


where D is the effective diffusivity that takes into account the chemical reaction given by 


eq.(1.1.1) [1.3]. The diffusivity shows an Arrhenius type of temperature dependence. 


] Basic results on diffusion and swelling 


Solution of the diffusion equation requires an appropriate boundary condition, very often 


chosen as constant surface concentration of molecular water: 


C(z =0,t) = C(0) = C, = constant (1.1.5) 


For a semi-infinite body, this assumption results in the following equation as a solution 


of the diffusion equation, 7.e., eq.(1.1.4): 


C= Ced 


Z (1.1.6) 
24 Dt 


which is well established for liquid water as the environment [1.4]. 


In the case of water vapour as the environment, the behavior of water at the silica sur- 
face does not follow eq.(1.1.5). Experimental results show that the concentration of 
water is not a constant at the surface, but increases continuously with time, until it reach- 


es the solubility limit of water in the silica glass [1.5, 1.6]. 


Experimental results for C(z- 0, d) are in clear contradiction with the assumption of a 
constant surface value for the water concentration. There is, however, an incredible 
agreement with the thermal analogue of thermal shock behaviour under heat transfer 
boundary conditions, governed by a heat transfer coefficient. This calls for a surface 
condition (see e.g. [1.7], Section 4.7) 


E_Rc-c) at =0, (1.1.7) 


d D 


where again G is the concentration of molecular water reached at z-0 for o. 


A swelling effect in water-containing silica at high temperatures was early reported by 
Shelby [1.8]. This author showed own and literature results on density reduction by water 


uptake at high temperature 


Ap. 
Po 


> Py = FC, =0) (1.1.8) 


w 


where Cy is the weight fraction of water and y an empirical coefficient. The observed 


volume swelling strain £v reads in terms of hydroxyl concentration S$ 


] Basic results on diffusion and swelling 


é,=KxS (1.1.9) 
with 
x =0.97 [0.92,1.02] (1.1.10) 
and the numbers in brackets indicate the 95% confidence interval. 
The effect of the water-silica reaction is the generation of swelling strains. Swelling 
stresses are a consequence of the mechanical boundary conditions. A volume element in 
a thick plate that undergoes swelling cannot freely expand. If the diffusion zone is small 


compared to the component dimensions, expansion is completely prevented in the plane 


of the surface and can only take place normal to the surface plane. 


At a free surface, the stress state is plane stress and, consequently, also stresses caused by 


swelling are equi-biaxial (0.70) 


Ow =O wz = E > O wh = AE (1.1.11) 
S ; 3(1-v) i 9(1-v) 


Osw n= hydrostatic stress. In (1.1.11) Zis Young's modulus and v is Poisson's ratio. 


Water concentrations at silica surfaces below 500°C are available from investigations by 
Zouine et al. [1.4]. For these measurements the nuclear reaction analysis (NRA) was 
applied that subsumes molecular and hydroxyl water. From these measurements the wa- 


ter species Sand C'can be obtained in molar units by 


C,=C+1S=C(l+!k) (1.1.12) 
where the quantity X is the equilibrium constant describing the ratio of A=S/C, Fig. 1.1. 


Due to the swelling stresses, the diffusivity is a function of water concentration. The 
diffusivity for the case of stress-enhanced diffusion is given by the following equation 
[1.9] 


AV, 
D- Dyess} o, 24 (1.1.13) 


where Do denotes the value of the diffusivity in the absence of a stress. Tis the absolute 


] Basic results on diffusion and swelling 


temperature in X; AV, is the activation volume for stress-enhanced diffusion and R is 


the universal gas constant. 
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Fig. 1.1 Results on surface concentration of the molecular (blue) and hydroxyl (red) water species ob- 


tained computed from measurements by Zouine et al. [1.4]. 


When an externally applied stress Oappı and swelling stresses Osw are present simul- 


taneously, it is 


A 
D=D, e). + Sgn) ae (1.1.14) 


Also the hydroxyl] concentration S resulting from the reaction (1.1.1) is stress dependent. 


Under uniaxial tension in z-direction it holds [1.1,1.10] 


s= spex Sita ) (1.1.15) 


RT 


] Basic results on diffusion and swelling 


where A Ver is an effective activation volume. We obtained for the case of stress-en- 
hanced swelling [1.1, 1.10] 


AV... — 14.4cm?/mol +10% (1.1.16) 


eff 
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2 Proof of swelling via crack 
terminating angles 


2.1  Terminating angle 


In earlier papers [2.1,2.2] we described a method for the identification of stresses in the 
surface region by the observation of crack-terminating angles of Vickers indentation 
cracks. The procedure was applied to chemically strengthened soda-lime glass [2.3] and 


surface layers after water soaking undergoing ion exchange [2.4.2.5]. 


The reaction eq.(1.1.1) results in a volume swelling [2.6] and since the free expansion in 
the thin layer is restricted in negative swelling stresses [2.7]. Such stresses at a free spec- 
imen surface must affect the crack front contours of cracks ending at these surfaces. This 
Section deals with the determination of shielding stress intensity factors from the termi- 


nating angles. 


a) 
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Fig. 2.1 Crack fronts terminating at the specimen surface under an angle q, a) crack retard in a zone of 
compressive stresses, b) crack advance by tensile stresses. Arrows indicate crack growth 


direction. 


2 Proof of swelling via crack terminating angles 


Figure 2.1 shows a crack growing from left to right in a bar with residual stresses in thin 
surface layers. The crack front terminates at the free surface under an angle q. If 
compressive stresses (expansive strains) occur at the surfaces, the actual crack front in a 
crack growth test under superimposed external load must stay behind (Fig. 2.1a). In 
contrast, tensile stresses caused by shrinking effects must result in an advance of the 
crack (Fig. 2.1b). 


2.2 Stress intensity factors 


The problem of a crack-front intersecting a free surface has often been studied in 
theoretical fracture mechanics and was also discussed in [2.1]. From curves reported in 
[2.8] it can be concluded that the '2-singularity of stresses characterizing the stress 


intensity factor is only possible for a crack terminating angle of 
o = 90° — 38.8?v (2.2.1) 


For most brittle materials with v=0.25, the crack terminating angle is about 80^, i.e. 
deviating by 10° from the normal. For silica with v=0.17 we have to expect: @=83.4°, 


that may be abbreviated as qo. 


Stress intensity factors for DCDC-specimens were determined by 3-dimensional FE 
computations on materials with Poisson ratios in the range of 0<v<0.25. In [2.1] the 
results were reported for v=0.25 as fits to soda-lime glass. Similar results for the special 
case of silica with v=0.17 are represented in Fig. 2.2a for crack-terminating angles of 
@=45°, 60°, 83.4°, and 90? in form of the fracture mechanics geometric function =F, 
defined for the DCDC-specimen by 


(2.2.2) 


p- K 
|p| JR 


In (2.2.2) p is the pressure at the end surfaces and R is the radius of the drill hole. The 
FE-results are introduced as the circles. We computed stress intensity factors along the 
crack front for different terminating angles. For estimations on a wide range of 
terminating angles two theoretical limit cases were included. For ọ=0 it must hold 
K(o)>» and for @=180°: K(@)=0 (see e.g. Fenner and Abdul Mihsein [2.9]). The stress 


2.2 Stress intensity factors 


intensity factors as a function of the terminating angle are given in a normalized 


representation in Fig. 2.2b recommended for interpolations in the region @<90°. 


For the applied stress intensity factor Kappi(@), the result of curve fitting reads for v=0.17 


(silica) 


Ka (P) ~ 1 odd 
K (90°) 0.00658 p + 7.66 x 107? g* 


and for v=0.225 (soda-lime glass) 


K appi (p) = 1 DEUS 
K opi (90°) 0.00606 o - 8.54x 102 p° 


These dependencies are shown in Fig. 2.2b as the curves together with the data points for 


v=0.17 as the circles. 


(90°)/K(9) . 


o K 


F 0.2 
uet v=0.17 
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Fig. 2.2 a) Stress intensity factors in the surface region of cracks terminating under different angles ọ in 
terms of the geometric function Faccording to eq.(2.2.2), b) surface stress intensity factor nor- 


malized on the value for p=90°, curves: interpolation equations (2.2.3a, 2.2.3b), circles: FE re- 


sults for v=0.17. 


2 Proof of swelling via crack terminating angles 


When residual stresses are present at the surfaces at which a crack terminates, they affect 
the crack tip stress field by generating an additional stress intensity factor. Since for 
compressive stresses the crack tip is partially shielded from the externally applied 


stresses, it may be called “shielding” stress intensity factor, Ki. 


Following the superposition priciple of linear-elastic fracture mechanics, the total stress 
intensity factor Kip is the superposition of the applied one, Kapp, and the shielding stress 
intensity factor, Ki. 


K +K, =K, (2.2.4) 


appl sh tip 


Under stable crack growth conditions, the total stress intensity factor as the sum of must 
equal the fracture toughness Kic. Then the terminating angle at an untoughened surface, 
Ka=0, must simply fulfil Kip- Ki. The fracture toughness of silica is at room temperature 
Kip =Kic = 0.75-0.8 MPaYm as reported by Wiederhorn [2.10]. Under conditions of sub- 
critical crack growth with crack-growth rates of about 10^ m/s, the stress intensity factor 
is roughly Kip = 0.57 MPa\m [2.11] that has to be introduced on the right-hand side of 
eq.(2.2.4). Then it results 


Kip _ K(g) (2.2.5) 
Ki, Kp) 


with the terminating angle @o in the absence of shielding. Knowledge of terminating 
angles c, qo, and Kip from the experiment allows the computation of Ki via eqs.(2.2.4) 
and (2.2.5). 


2.3 Experimental determination of shielding stress 
intensity factors 


We studied crack terminating angles on DCDC-specimens [2.12] made of the silica glass 
ENOSNB (GVB, Herzogenrath, Germany) containing 99.98% SiO». DCDC-specimens 
were machined and annealed for 1h at 1150?C in order to remove residual stresses. Then 
the specimens were heat-treated in water vapour at 250?C for 92h under saturation pres- 


sure. Afterwards, fast DCDC-tests under increasing load were performed with partial 
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2.3 Experimental determination of shielding stress intensity factors 


unloading for the generation of crack-arrest markings on the fracture surface. Final frac- 
ture was obtained by introducing a needle with an excessive diameter into the DCDC- 
hole. Figure 2.3a shows the result for a specimen without heat treatment. The terminating 
angle of @o=83°, expected from eq.(2.2.1), is tentatively introduced by the straight lines 
showing sufficient agreement with the crack contours. A specimen tested after a 192h, 
250°C vapour-treatment is given in Fig. 2.3b. This image shows clearly deviating ter- 
minating angles of only @=25° as is indicated by the straight line. In addition, we tested 
DCDC-specimens soaked in liquid water for 48h and 192h at 250?C. 


Table 2.1 compiles the relevant test data and the experimental findings. The crack termi- 
nating angles are listed in Column 4 as the average angle. The layer thicknesses for the 
liquid-water soaked specimens are compiled in Column 5. The layers for the water- 
soaked specimens may be slightly smaller, because some dissolution of silica in water 
cannot be excluded. The approximation signs may indicate this. Column 6 gives the 
stress intensity factors for the crack extension under subcritical and spontaneous crack 


growth conditions. 


Test conditions t oO) b Kappi.o Kappi Kn 
(um) | (MPavm) | (MPavm) | (MPa\m) 
Silica = 0 83.4 0 0.57 0.57 0 
vapour 192h | 25 17.7 2.94 -2.37 
liquid 48h | 41.0 | «89 1.77 -1.20 
liquid 192h | 29 | =17.7 0.8 3.44 —2.64 
Soda Chem. 13.5 0.75 9.16 —8.41 
-lime hardened 
water 0 107 0 0.75 0.75 0 
7d 71 0.34 2.68 -1.93 
35d 43 0.52 5.44 —4.69 
42d 45 0.52 5.17 —4.42 


Table 2.1 Experimental results; Kippio =K(@o) is the applied stress intensity factor in silica specimen without 


heat treatment [2.2], results on soda lime glass from measurements in [2.1]. 


11 


2 Proof of swelling via crack terminating angles 


In the case of liquid water tests for the 250?C/192h-tests, we evaluated the first crack 
contour as is visible after spontaneous crack propagation starting from the hole. The 
applied stress intensity factors Kapp are given in Column 7. Strongly negative shielding 
stress intensity factors Ksn»<O are observed in all cases, indicating strong compressive 
stresses in the water diffusion layer, Column 8. In all cases, for silica as well as for soda- 


lime glass the materials undergo volume expansion in the water-affected surface layer. In 


[2.13] it is roughly shown in which way the compressive swelling stresses can be deter- 


mined. 


Fig. 2.3 Crack terminating angles in DCDC-specimens, a) without hot-water soaking (theoretical angles 
of 9783.4? tentatively introduced as the tangents), b) after 192h, 250°C soaking in water vapour 
at saturation pressure, c) after 192h 250?C soaking in liquid water.; One arrest line in a) and b) 


stand out as black contour, crack growth direction indicated by the arrows. 
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3 Effect of soaking on inert strength 


3.1. Inert strength 


In earlier publications [3.1-3.4], we explored the idea that water can toughen silica glass 
by diffusing into the glass structure from the crack tip. This process sets up a negative 
stress intensity factor that shields the crack tip and enhances the strength of the specimen. 
Because diffusion rates increase with temperature, crack tip shielding should also 
increase as the temperature is increased, as should the specimen strength. Using the 
model presented in [3.2], we explored the mechanism of crack tip shielding and showed 
that the calculated changes in strength resulting from water exposure at 88 °C agreed 
sufficiently with experimental values measured on high-silicate glass Vycor by Ito and 
Tomozawa [3.5], and Hirao and Tomozawa [3.6]. The effects on inert strength at this 
rather low soaking temperature are only about 10%. In this paper, we will address the 


same phenomenon, but at higher temperatures of 8=260°C. 


Strength measurements at 250?C are in principle known from literature. Li and Tomoza- 
wa [3.7] soaked silica bars for up to 4 days soaking time and tested the strengths in dy- 
namic bending tests under subcritical crack growth conditions. To the authors' 
knowledge, so far no really inert strength measurements are available for water-soaked 


silica. 


3.1.1 Water diffusion in silica 


Water diffusion into the surface of silica glass has been studied experimentally by a 


number of investigators, and shown to depend on temperature according to 


D, = A,exp[-Q, / RT] (3.1.1) 
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3 Effect of soaking on inert strength 


where Qw is the activation energy, Tis the absolute temperature, and R is the gas con- 
stant. As reported in reference [3.8] for silica in the temperature range 0 °C to 200 °C: Qw 
= 72.3 kJ/mol, logio 4o — -8.12 (Ao is in m?/s). 


The diffusion distance 5, an appropriate measure for the thickness of the diffusion zone 
(where the water concentration is roughly half of that at the surface) is given by 


b=./D.t (3.1.2) 


w 


(time). 


3.1.2 Volume swelling 


The relation between volume swelling strain £v and the hydroxyl concentration $ (in 


weight units) is given by eq.(1.1.9) 
€,=KxS , «=0.97 (3.1.3) 
The hydroxyl water S concentration can be computed from Section 4.1 in [3.9] 


0.000780 exp(0.00868 9) (3.1.4) 


S — A7 
e | 
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1 +0.031ex 
coste [12/005 


A volume element near a surface that undergoes swelling cannot freely expand. If the 
diffusion zone is small compared to the component dimensions, expansion is completely 
prevented in the surface plain and can only take place normal to the surface. The swelling 


stress at the surface, Osw,o is then given by eq.(1.1.11). 


3.2 Total stress intensity factor for a semi-circular 
surface crack 


As in [3.4] we begin our analysis by estimating the magnitude of the total stress intensity 


factor, Kio, experienced by the critical crack during the strength test. The total stress 
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3.2 Total stress intensity factor for a semi-circular surface crack 


intensity factor consists of three contributions: Kapp, which is primarily the result of the 
applied stresses and the geometry of the crack; Ks, which results from water penetration 
into the fresh fracture surfaces and the other surfaces of the fracture specimen. In the 
absence of subcritical crack growth the condition for crack extension is Ki 2: Kic. 

In order to provide a transparent analysis we restricted the possible aspect ratios of cracks 


to the commonly chosen semi-circular surface crack. The depth is denoted as a (Fig. 3.1) 


and the aspect ratio as a/c=1. The plane of the crack is assumed to be perpendicular to the 


specimen length axis. 


The applied stress intensity factor given by eq.(3.2.1) and eq.(3.2.2) were derived and 
discussed in [3.10] for the case of straight specimen surfaces. Since the initial natural 
surface cracks are very small compared to the specimen thickness W, a/W<<1, it holds 


for the stress intensity factor at the deepest point (A) 
Kaya = 04,1173 a (3.2.1) 
and for the surface points (B) 


Kapi,B 23 c.,,129 Ja $ (3.2.2) 


where in bending tests Gappi is the outer fiber tensile stress. This solution can also be used 
for cylindrical specimens if the crack depth is small compared with the cylinder radius, 
a/R<<1. The crack is assumed to be at the midpoint of the test specimen, with one axis 
perpendicular to the specimen length axis and the other principal axis in the plane of the 


specimen surface. 


As noted above, two contributions make up the shielding stress intensity factor: one 
coming from the diffusion zone originating from the crack faces, the second originating 


from the external surfaces of the specimen. 


The shielding stress intensity factors at the deepest point A and the surface points B, Ku 
and Kss, are for b««R [3.4] and 0.15<b/a<1.25 


sw,0 


K,,,z1.17Jao LU + 03172] (3.2.3) 
d a a 
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3 Effect of soaking on inert strength 


Ky, 2129 a 0,,, tan 27 + 0002| (3.2.4) 
' : a a 


Fig. 3.1 Single failure relevant crack in a semi-infinite body exhibiting a diffusion and swelling zone, a) 


side view, b) top view on the crack surface. 


3.3 Toughness and inert strength 


3.3.1 Apparent Toughness 


The total stress intensity factors including shielding are given by 


K oa = Kona +K (3.3.1) 


appl,A 


K „=K .+K (3.3.2) 


sh,B appl,B 


where Kappa and Kapps are obtained from eq.(3.2.1) and eq.(3.2.2); Kena and Kis, from 
eq.(3.2.3) and eq.(3.2.4). The applied stress intensity factors for inert tests have to be 


computed using the initial crack dimensions æ and æ. 
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3.3 Toughness and inert strength 


The shielding effect is responsible for an apparent increase of the fracture toughness 
which is, in general, identified with the applied stress intensity factor at failure, Kappicr. 


From eqs.(3.3.1, 3.3.2) the apparent fracture toughness, in the following denoted as K. ; 


simply results as 


ER EK AR. (3.3.3) 


Since Ksn<0, it holds Re >Kic. On the other hand it becomes clear from the fact that the 
inert strengths are proportional to the applied stress intensity factor, oc œ Kappicr, that 
strength and apparent toughness increase by the same factor. In the following considera- 
tions, we therefore concentrate on the strength, exclusively. The apparent toughness for 


the surface points and the deepest point of a semi-circular crack is shown in Fig. 3.2a. 


^ 
Kic | Oc 
2.5} 500F 
(MPavm)!} (MPa) | b) 
| a) H 192h 


2r | soaking time 


96h 


t=24, 48, 96, 192h 


15r 
300 
1k 
L 200 
x [ point A 
[ —— pointB 
—————————— E 
100 150 200 250 
0 co 
Fig. 3.2 a) Apparent fracture toughness after hot-water soaking for a semi-circular crack of 30um depth, 


b) calculated inert strengths as a function of water-soaking temperature 0 and soaking time, pre- 


dicted for inert strengths of o-=100, 150, and 200 MPa for freshly abraded test specimens. 
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3 Effect of soaking on inert strength 


3.3.2 Inert strength 


Using the temperature dependent swelling strain together with the temperature effect on 
the diffusivity, the inert strengths could be computed as a function of 0 in the same way 
as outlined in [3.4]. The inert strength, oc, is given by the condition 

c, LC -Kha Ky es (3.3.4) 


ii74a ^ 129Ja 


Strength predictions are shown in Fig. 3.2b. For this purpose initial inert strengths were 
assumed to be o-=100, 150, and 200 MPa. 


3.4 Strength of unsoaked and soaked specimens 


We studied the strength behaviour of the silica glass ENO8NB (GVB, Herzogenrath) con- 
taining 99.98% S1O» [3.11]. Cylindrical bending bars of 45 mm length were cut from 
silica rods of 4 mm diameter maintaining the original surface from manufacturing. Then 
all specimens were annealed for 1h at 1150?C in vacuum in order to eliminate pre- 
existing residual surface stresses. To avoid any water contact, the series foreseen for 
strength measurements in an inert environment were immediately after cooling down 
stored in fresh silicon oil as proposed by Sglavo and Green [3.12] as an inert medium. In 
total, two glass deliveries were considered, denoted in the following as Batch (I) and 


Batch (II). The soaking conditions and surface zone parameters are given in Table 3.1. 


Batch Ref. Time Water Vapour b S Ev Osw,0 
(h) CC) CC) (um) | (wt%) | (%) | (MPa) 

(1) [3.11] 192 263 - 21.7 0.85 0.82 -240 
(D [3.13] 192 250 250 17.7 0.74 0.72 -209 


Table 3.1 Soaking conditions and properties of the water layer at the surface. 


Inert bending strength tests were carried out in liquid nitrogen in a 3-point testing device. 
The results for unsoaked specimens (open circles) are shown in Fig. 3.3a in Weibull 


representation. The distribution of the unsoaked specimens shows the usually expected 
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3.4 Strength of unsoaked and soaked specimens 


Weibull straight-line with a Weibull modulus of m= 12.3 and a characteristic strength of 
oo = 171.3 MPa for Batch (I) both obtained by application of the Maximum Likelihood 
procedure [3.14]. For Batch (II) the strength distributions are shown in Fig. 3.3b. For the 
untreated material it was found m= 6.6 and oo = 250 MPa. It should be noted that for 
Weibull-distributed initial strength values, the strength under swelling stresses cannot 
result in a Weibull distribution, too. The reason is due to the fact of varying apparent 
fracture toughness values Ky according to eq.(3.3.3), whereas the derivation of the 
Weibull distribution needs a constant Ki. (for details see e.g. Section 10.3 in [3.15]). The 


parameters derived before may be denoted as apparent Weibull parameters. 


Even the strengths of the unsoaked specimens appear rather large. The reason for this is 
of course the 3-point bending loading with its reduced effective surface compared to 4- 
point bending tests and on the other hand the circular cross section of the specimens, that 
reduces the effective surface once more. The strengths could be described roughly by 
Weibull distributions except of Batch (I) after water soaking. Table 3.2 shows the related 
Weibull parameters. Since Batch (I) shows clearly deviating behaviour, the median val- 
ues as the strength parameters are compiled in Table 3.3. From the strengths in the un- 
treated state and the fracture toughness of Ki.=0.8 MPavm [3.16], the initial crack depth 


ao can be concluded as shown in Fig. 3.4a and 3.4b. 


Specimens Batch (I) Batch (II) 
60 (MPa) m 60 (MPa) m 
untreated 171.3 12.3 250 6.6 
Water soaked - - 395 4.5 
Vapour soaked - - 340 8.1 


Table 3.2 Strength results, represented by the Weibull parameters. 


Specimens Batch (I) Batch (II) 
untreated 160.6 MPa 232 MPa 
Water soaked 367 MPa 373 MPa 
Vapour soaked - 316 MPa 


Table 3.3 Strength results, represented by the median values. 
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3 Effect of soaking on inert strength 
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Fig. 3.4 Distributions of the initial crack depth ao. 
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3.5 Strength predictions 


3.5 Strength predictions 


Measurements in liquid nitrogen showed a clear increase of the inert strength for both the 
water-soaked and the heat-treated specimens in saturated water vapour. The results of 
Fig. 3.3 are plotted once more in Fig. 3.5 together with the predictions according to 


eqs.(3.2.1-3.3.4). 


The predictions were made on the basis of the Weibull distribution for the unsoaked 
specimens. First, the strength of the 7^ specimen was computed for a total number of N 


specimens via 


om -o(n—— | (3.5.1) 
1-(n-0.5)/N 


For each selected specimen, the eqs.(3.2.1-3.3.4) were applied resulting in the predicted 
strength of soaked specimen corresponding to the same failure probability F This pro- 


cedure yields the red data in Fig. 3.5. The red arrows stand for the strength increase. 
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Fig. 3.5 Measured strengths (solid circles and squares) and predictions by the shielding stress intensity 


factor (red circles). 
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3 Effect of soaking on inert strength 


In the case of Batch (I), the lowest four strength data after soaking obviously differ clear- 
ly from the prediction showing typical features of a bimodal strength distribution. Due to 
their high strengths, the soaked specimens shattered in more than two fracture pieces 
(mostly 4-6 fragments). Therefore, it was not possible to carry out a fractographic study 


on the broken test pieces. 


3.6 Analysis of nonlinear strength distributions 


In the following considerations on Batch (I), it is assumed that the soaked specimens 
exhibit a bimodal strength distribution caused by competing of two different strength 


populations. Then the Weibull distributions of the strength populations “1” and “2” are 


given by 
F, -1- exp (2 | (3.6.1) 
L O9 
and F, =1-exp | 9c (3.6.2) 
L Op 


with the two different characteristic strengths oo: and oo and the related Weibull expo- 
nents mı and mp. Superposition of these strength results in the total failure probability 


(for details see e.g. Section 8 in [3.15]) 


F =F,+F,-FF,=1-exp & E | (3.6.3) 
On Oo 


The Weibull parameters obtained by curve-fitting according to eq.(3.6.3) are compiled in 


Table 3.4. The curve fit according to eq.(3.6.3) is introduced in Fig. 3.6 by the solid 


curve. The dash-dotted straight lines give the asymptotes representing the individual 


strength populations. 


The predicted strength made by using the shielding stress intensity factor is close to the 


measured ones of “population 1”. 
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3.7 A possible reason for the bimodal strength distribution 


Strength population *1" strength population “2” 
m 601 (MPa) nn 60» (MPa) 
21.4 437.5 [404, 471] 2.7 394 [306, 482] 


Table 3.4: Weibull parameters for the water-soaked specimens of Batch (I), 90% CI in brackets. 
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Fig. 3.6 Tessellation of the strengths for the hot water soaked specimens in two strength distributions. 


3.7 A possible reason for the bimodal strength 
distribution 


From the nearly linear shape of the strength distribution of unsoaked specimens in the 


Weibull plot, it can be concluded that for these results the initial cracks responsible for 


failure were Weibull distributed, too. The bimodal strength distribution of the soaked 


specimens calls for an interpretation by a bimodal distribution of the failure relevant 


flaws. In order to discuss the observed curve shape, let us assume coexistence of the sur- 
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3 Effect of soaking on inert strength 


face defects (cracks) and inner flaws (pores) as schematically illustrated in Fig. 3.7a. The 
blue line indicates the strengths for the surface cracks, whereas the red line gives the 


strengths for the inner flaws. 


Under normal circumstances, the specimens will fail at the surface, because the strength 
is lowest there. Failure starting from internal defects is very seldom. This is possible only 


when the two curves intersect at a very small failure probability F. 


Due to soaking, the surface is strengthened by swelling stresses. Consequently, the blue 
line shifts to higher strength values as indicated by the arrow in Fig. 3.7b. The defects 
below the surface zone of about 20um remain unaffected and can cause failure. The 
measurable strength is the minimum of the shifted blue and the red line. This results in 
the kinked black curve as found in our measurements. Since no fractographic examina- 
tion of the fracture surface was possible for the specimen fractured in liquid nitrogen, we 
inspected the specimens for internal flaws with a confocal microscope. Figure 3.8 shows 
subsurface defects that cannot be introduced by handling of the specimens. Their radii 
are about 5-10um. At least by this inspection technique we can confirm the existence of 


internal flaws competing with the surface defects on the lateral surface. 


Inin(1/(1-F)) InIn(1/(1-F)) 
1r 1f ! 
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Fig. 3.7 a) Strength for failure starting at the surface flaws and at inner defects, b) shift of the surface 


strength by soaking resulting in a bimodal strength curve. 
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3.7 A possible reason for the bimodal strength distribution 


[1]21um]| 


Fig. 3.8 Images of subsurface defects obtained with a confocal microscope. 
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Fig. 3.9 Results of disk tests ([3.17] and Section 2 in [3.9]) and from inert strength measurement in liquid 
water and vapour under saturation pressure, compared with predictions based on hydroxyl con- 


centrations in liquid water. 
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3 Effect of soaking on inert strength 


Finally, Fig. 3.9 shows the shielding stresses from inert bending strength measurements 
together with the predictions via S. The strengths for water soaked specimens are only 
slightly above the prediction and the strength with possibly internal failures somewhat 
below the predictions. The data point for the vapour-tests clearly falls below the predic- 


tion as was found already for the evaluation of crack-terminating angles in Section 2. 
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4 Strength in silicon oil, humid air and 
water 


4.1 Strength measurements on silica 


Very often strength measurements are carried out in humid environments. This may lead 
to subcritical crack growth effects on strength. In this Section the effect of swelling 
stresses on strength is discussed for subcritical crack growth conditions. An interesting 
effect is the occurrence of discontinuities in the load-displacement curves for strength 
tests performed in silicone oil as is generally assumed to be an inert environment, 


preventing the access of water to the crack tip [4.1 ]. 


Strengths were measured on silica glass ENOSNB (GVB, Herzogenrath) containing 
99.98% SiO». Bars 3x4x45 mm? were cut out of a plate 350x350x4 mm? and surface 
machined by grinding with a grinding wheel, D91-C75. To eliminate residual surface 
stresses introduced by grinding, we annealed our specimens in vacuum for 1h at 1150?C. 
The series for strength measurements in an inert environment was immediately stored in 
fresh silicon oil as an inert medium after cooling to room temperature, following the 
procedure of Sglavo and Green [4.2]. In our tests, we used the original silicon oil 
(Wacker AK 100, Wacker-Chemie AG, München), without an additional drying 


procedure. Each test series was carried out with fresh oil. 


Four-point bending strength tests with the normal rectangular bending bars were made in 
silicone oil, as recommended by Sglavo and Green [4.2]. These investigators dried silica 
glass rods at 120?C for 2h before being completely immersed and fractured in a silicone 
oil bath. This procedure was considered in [4.2] to be an inert environment, preventing 


the access of water to the crack tip. 


In all the tests the loading rate was 50N/s (42 MPa/s) resulting in test durations of 2-3 
seconds. The average strength of 15 tests, resulted in o.=103.5 MPa (SD 9.2 MPa). 
Under the assumption of half-penny shaped initial cracks, the average depth a of the 
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4 Strength in silicon oil, humid air and water 


introduced crack during grinding, can be calculated from Kr and the measured strength, 


Oc. 


A, E (4.1.1) 
oY 


By using Y&1.3 and K,=0.8 MPa\m we obtain a=35 pm. 


Two series of specimens were then soaked for 24h in liquid water at 250°C in an 
autoclave. Drying the specimens in vacuum removed surface moisture. After 2 days of 


drying at 60°C in vacuum, the average strength in silicon oil was found to be o,=119.7 
MPa (SD 9.9 MPa). 


Load 
soaked 
annealed 
Displacement 
Fig. 4.1 Load displacement plots in silicon oil. The water-soaked specimens showed a load pop-in that was 


missing in the tests on unsoaked bars. 


The load-displacement curves obtained in the strength tests in silicon oil were completely 
different from those obtained from the specimens that were not exposed to saturated 
steam at 250 °C. Only the unexposed annealed specimens showed the usual straight 


displacement versus load behaviour until failure. All specimens exposed to water at 
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4.1 Strength measurements on silica 


250 ?C showed pop-ins as schematically indicated in Fig. 4.1, accompanied by clearly 
audible cracking noise. The fracture surface of one such test specimen is shown in Fig. 
4.2a and Fig. 4.2b. Figure 4.2b represents the fracture origin. At this magnification, the 
initial crack of depth a=31um is clearly visible. The same holds for the increased crack 


after the load pop in. 


IAM-KM 188un 
1-Aug-2812 Detector- Si EHT-28 .88 
Mag 


18un 
Detector 


Fig. 4.2 Fracture surface of a water-soaked specimen measured in silicon oil, a) overall view, b) 


identification of the initial crack size a=31 um and the extended crack after the pop-in (a-48um). 
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The arrest contour gives a new depth of 48um from which later final fracture starts. The 
crack ends at the surface are not visible, therefore the evaluation is limited here to the 


deepest points of the cracks, represented by semi-circular surface cracks, a/c= 1. 


The stress intensity factors related to the failure stress can be computed via FE-results of 
Newman and Raju [4.3]. Since the initial natural surface cracks are very small compared 
to the specimen thickness W, a/ W<<1, the solution holds for the stress intensity factor at 


the deepest point (A) of the semi-circle: 


Kappi, = O pp11:173Va (4. 1 .2) 


From the strength of of=119.7 MPa and the outer crack contour a=48 um at which 
catastrophical failure starts, we obtain Kappa = 0.97 MPavm, that is slightly above the 
fracture toughness of Ki = 0.8 MPawvm [4.4]. This results in a shielding term at the 
deepest point of the critical crack of Ks = 0.97—0.8 =—0.17 MPaYm. 


The applied stress intensity factor at which the crack extended from its original depth of 
a= 31um occurred at about 80-85% of the strength, i.e. at the stress where the load pop- 
in was observed. This results in an applied stress intensity factor of Kappa = 0.61-0.63 
MPavm that is clearly below Aj. Consequently, crack extension from a to a. must be 
caused by subcritical crack growth. The fracture behaviour, Fig. 4.2b, will be addressed 


below in the discussion Section 4.2. 


4.2  Pop-in behaviour 


4.2.1 Subcritical crack growth in silicon oil tests 


Measurements in silicon oil resulted in real inert strengths only in the case of the 
annealed un-soaked specimens. In strength tests on soaked specimens, also carried out in 
silicon oil, we found clear evidence for the occurrence of subcritical crack growth, even 


though the environment had prevented any water supply. 


Now let us look for a water source in the absence of any supply from the environment. 
Water that diffused into silica during the water soaking procedure reacts with the silica 


network according to eq.(1.1.1). For 250 °C, the ratio S/C is about 2.7. From Fig. 4.4a in 
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[4.5] we get by a slight temperature extrapolation to 250°C a total water concentration at 
the surface of C,~0.7 wt-percent. Consequently, the surface concentration is C«0.3 wt-% 
of molecular H;O. This shows that molecular water, necessary for bond splitting during 
subcritical crack growth is in the soaked specimens always available ahead a crack tip. 
The distribution of molecular water after 24h soaking at 250°C is illustrated in Fig. 4.3. 
Under these soaking conditions it results from eqs.(3.1.1) and (3.1.2): 5=6.3 um shown 
in Fig. 4.3 by the dashed contour. The C-profile along the x-axis as the prospective crack 
plane is concluded from FE-results in [4.6]. The content of molecular water is negligible 


for x» 3b. 


For strength tests in water, the cracks grow also strongly accelerated due to the effect of 
escaping from the initial swelling zone (Section A1 in [4.5]). However, in contrast to the 
silicon oil tests, there is now no crack arrest possible since water is always available to 
the advancing crack tip. Consequently, in these tests accelerated cracking is directly 
followed by final fracture without any further increase of loading necessary. The pop-in 


is therefore missing. 


15 20 (um) 


Fig. 4.3 Distribution of molecular water in the surrounding of a crack tip after 24h water soaking at 250°C. 
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4.2.2 Interpretation of crack extension accompanied by load 


pop-ins 


A side-view of the semi-circular surface crack is shown in Fig. 4.4a. From all surfaces 
which are in contact with water vapour during high-temperature soaking, water diffuses 
into silica. This holds also for crack surfaces. For reasons of transparency in the 
computation of the shielding stress intensity factor these contours are approximated by a 


zone of constant thickness b ending in a half-circle. 


If the crack grows under increasing load by an amount of Aa, it will escape from the 
initial swelling zone, as illustrated by the dashed crack extension. The related shielding 
stress intensity factor is schematically plotted in Fig. 4.4b as a function of crack 
propagation Aa. When Aa-b is reached and the crack tip leaves the swelling zone, the 
shielding stress intensity factor strongly decreases as was shown in [4.7] for ion- 
exchange layers in soda-lime glass. Superposition of the applied stress intensity factor 
Kappi © Va (dashed curve in Fig. 4.4c) and the shielding stress intensity factor Ky, results in 
the total stress intensity factor Ku Kup as given by eq.(3.3.1). 


The development of Ki, with increasing crack length, a+Aa, and increasing load is sche- 
matically shown in Fig. 4.4c by the solid curve. Under increasing load the first positive 
total stress intensity factor, Kmin = - Ka, is reached at Aa=0 for the curve indicated by Ormin 


which is the minimum possible strength in a dynamic strength tests 


-Ko (4.2.1) 


O fmi "I. [— 
= 117a, 


From the shape of the Kap-A^Aa-curve, it is evident that a crack extension by subcritical 
crack growth must accelerate rapidly in the region where Ki increases rapidly, namely at 
Aa b. 


On the other hand, we have to take into consideration the fact that simultaneously with 
the loss of shielding also the available water for subcritical crack growth disappears after 
a crack extension of about Aaxb. This results in an abruptly arresting crack so far Kiip<Kic. 
Final cracking can occur only after a sufficient increase of the loading. This trivially must 


result in the observed load-displacement curve of Fig. 4.1. 
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4.2.3 Interpretation of the shieldig stress intensity factor 


The shielding stress intensity factor at the deepest point of a semi-circular surface crack 


soaked from the side surface and the crack surfaces, Fig. 4.4a, is given by [4.8] 


E 
K pa Z 7a xt ean osoa E + sav?) (4.2.2) 
—Vv a a 


b 8) Ksna b) 


TE» 


2 
Aalb 


b Aa 


Fig. 4.4 a) Crack propagating out of the initial swelling zone, b) shielding stress intensity factor during 
crack escape, c) increase of the total stress intensity factor Ki during crack growth for three 
applied loads (schematic). Below the applied stress Oappi=Ofmin defining a lower limit for 


measurable strengths, no subcritical crack growth is possible because of Kior<0. 


For the crack with depth of a=48 um and a volume strain at 250°C of £,2 0.72296 given 
by eq.(3.1.4), it results from (4.2.2) with E-72 GPa and v=0.17: 
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K,, = -0.48 MPa m (4.2.3) 


From the experiment we found a shielding term at the deepest point of the critical crack 
of only Ku=-0.17 MPa\m, i.e. 36% of the expected result, eq.(4.2.3). The reason is 
rather trivial. Equation (4.2.2) holds for a crack of depth a that is fully surrounded by the 


swelling zone as is schematically shown in Fig. 4.4a. 


Two contributions make up the shielding stress intensity factor. One is coming from the 
diffusion zone originating from the crack faces, Ka, the second originating from the 


external surfaces of the specimen, Ki. 


For a transparent interpretation, let us apply the K-separation according to [4.8]. The total 


shielding stress intensity factor is the sum 


Ky = =K,ıt Ka» (4.2.4) 
Konz 
ooVna -Ksna(1-v) 
| e,ENb 


Fig. 4.5 a) Influence of aspect ratio a/c and relative depth of the swelling zone b/a on the shielding stress 
intensity factors Kx at the deepest point of a semi-elliptical surface crack [4.8], b) shielding stress 


intensity factor Kn, for a crack growing out of the initial swelling zone, [4.9]. 
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4.2 Pop-in behaviour 


The stress intensity factor by the side surface is shown in Fig. 4.5a as a function of b/a 
and the crack aspect ratio a/c (2c= width of the crack). In our case results for a/c=.1 and 
b/a=6.3um/48um =0.13 (circle in Fig. 4.5a) with 6,, =—210 MPa: 


Ka2=-0.0844 MPavVm. 


When the crack can grow subcritically by an amount of Aa, it will escape from the initial 
swelling zone, as illustrated in Fig. 4.5b. The stress intensity factor Ki; for the swelling 


zone developing from the crack faces, is shown in Fig. 4.5b as given in [4.9], Fig. 15.6. 


The shielding stress intensity factor Kn, for an extension of Aa=48 um-31 um = 17 um 
and Aa/b=2.7 decreases (red arrow in Fig. 4.5b) to a normalized stress intensity factor of 


0.036. Consequently, it results 


K,,, = —0.036 < vE Jb = -0.0566MPaym (4.2.5) 
i -v 


The total shielding term after 17um crack extension results from eq.(4.2.4) as 
Ka= 0.0566 MPa\m — 0.0844 MPaYm =-0.141 MPaYm 


This is about 83% of the experimental value of —0.17 MPa\m. 
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5 Effect of swelling on dynamic 
strength 


5.1 Dynamic strength 


In Section 3 we considered the effect of volume swelling by hydroxyl generation on the 
inert strength of hot-water soaked silica. The strength behaviour was discussed on the 
basis of the total water concentrations from Zouine ef al. [5.1] and Wiederhorn et al. [5.2]. 
The hydroxyl concentration at silica surfaces is time dependent as can be seen from the 


studies by Davis and Tomozawa [5.3] and Oehler and Tomozawa [5.4]. 


Very often strength measurements are carried out in humid environments. This leads to 
subcritical crack growth effects on strength. In this Section, the effect of swelling stresses 
on strength is discussed for subcritical crack growth conditions according to [5.5]. 
Special attention will be drawn on the evaluation of the power law exponent n for 


subcritical crack growth as had very early been developed for glass by Charles [5.6]. 


For strength tests carried out in humid environment at various loading rates, so-called 
dynamic strength tests, one can show theoretically that the swelling effect caused by the 


reaction of water with silica must result in apparently increased crack-growth exponents. 


5.2 Computation of dynamic strength for edge-crack 
like defects 


5.2.1 Surface defects modeled by edge cracks 


The most transparent description of crack-growth problems can be given by using a 
solution for one-dimensional edge cracks. Under an applied remote stress o; the related 


applied stress intensity factor is 
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K „ =1.122,/za o (5.2.1) 


appl appl 


The shielding stress intensity factor for an edge crack of depth a is derived in Section 
A1.5 in [5.7] based on results in [5.8] 


2/3 
K,, =1.122/zao,,, ona " 4 0.832 a | (5.2.2) 
a a 


with the thickness of the swelling zone b given by eq.(3.1.2) and the swelling stress Oww.o 
at the surface. The crack-tip stress intensity factor results by superposition of the applied 


and shielding stress intensity factors: 


Kip = Kay + Ky, where K, «0 (5.2.3) 


appl 


5.2.2. Cracks fully embedded in the swelling zone 


If the crack is fully embedded in the swelling stress zone, 1.e., if a << b, it results from 
eqs.(5.2.1 — 5.2.3): 


Kp =1.122(0, 591 + 0,0470 - (5.2.4) 
Let us describe the subcritical crack growth rate, v, by a power law relation 


da á " 
v= dt = AK;, = A(K appi +K,,) (5.2.5) 


In the absence of swelling, 6,4070 and Ka =0, the following well-known equation holds 
(see e g. [5.9]): 


2-n 
Ge; = Bo", 6(n } 1), B = SEE v (5.2.6) 
, , 1.122? z A(n- 2) 


where oz is the inert strength, the strength in the absence of subcritical crack growth. In 


presence of swelling stress Osw,o, we obtain the following equations: 


(c, +0,,0)" = Bo; c(n«1) (5.2.7) 
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o,= (Bo? (n ü: Pyare giae nOD. (5.2.8) 
Oro 
>09, = fo Oo (5.2.82) 


The inert strength is in soaked specimens 


= K,, = Ky, (5.2.9) 
* 1122/za, 
and in un-soaked specimens it holds 
Be Kr (5.2.10) 
^— 14122478, 


where æ is the initial crack depth. It follows that the swelling stress is simply given as 
Oow 0 = Oeo ~ Fe" (5.2.11) 


Two conclusions can be drawn from eq.(5.2.8a) and eq.(5.2.11): 


Strength increase 


Since the inert strength after soaking exceeds that without soaking, an increase of 


strength in humid environment results equal to the increase of the inert strength 
0,-0,,70,-0,, (5.2.8b) 


Apparent power-law exponent n’ 


Whereas for the un-soaked material the exponent in the power-law relation, z, can simply 
be obtained from the slope of the o;— f( O ) plot, this is no longer possible for the water- 


soaked material as can be easily shown. Let us consider eq.(5.2.8) in the form 
aop =A +, Lo (5.2.8c) 
with the abbreviation 


A = (Bo? (n4 1)» (5.2.12) 
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A straight-line evaluation of eq.(5.2.8c) must result in an apparent crack exponent z 
defined by 


6, =- Ogg t AGM >, oc Mm (5.2.13) 
or in logarithmic representation 
1 = d(logo;) _ d[log (AG LE 0,,.9)] (5.2.14) 
n*1 d(logo) d(logo) 


By using logarithmic derivations, d(log x)=(1/x) d x, we obtain the following equation: 


n'=n+(n+1) Ee . (5.2.15) 
O fo 
O,—O 
n'x n4 (n41) — —LÀ 
Oro 


Since o,0<0 it follows that 7° >n. Thus, the apparently increased n’-value is an artifact 


of the logarithmic representation of a shift of a straight line by a constant value of -0;,o. 


5.2.3 Crack partially embedded in the swelling zone 


The treatment for crack sizes competing with the swelling zone size is rather complicated 
and nontransparent, since it requires numerical computations. One reason for this 
complication is the fact that the swelling stresses are now no longer constant but decrease 
from a maximum value at the surface continuously with increasing distance from the 


surface, roughly given as 
2 5.2.16 
Oy = Foy 9 Erfe — (5.2.16) 
: 2b 


As a consequence of the large z-value for silica (7230), eq.(5.2.5) makes clear that the 
largest stress increments do;ppı per crack length increment da occur at the initial crack 
depth ax a. This follows from dOap= 6 dt, and da=AK"dt via do.p/daxl/K”. The geo- 
metric function can be assumed to be constant Ab/a)=Rb/ao). In the presence of 


swelling, we obtain then the following equations: 
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[c, +F(2) o, I = ao (5.2.17) 


(5.2.18) 


sw,0 


>g, =A" F(*3g. 


The geometric function for an edge crack occurring in (5.2.2) was determined by FE- 


2/3 
F(b/a)z am foses f? + Ld | (5.2.19) 
a a 


Again a straight-line evaluation of eq.(5.2.18) results in an apparent crack exponent zr 
defined by 


computations as 


E b > 1/(n+1) = 1n'+1) 
e; SFG) yt AG >0,x6 (5.2.20) 


or equivalently 


1 d(logo,) dllog(4ó """— F(2)o,, )] 


(5.2.21) 
n+1 d(logo) d(logo) 
Again using logarithmic derivations, we obtain the following equation: 
-F(2Jo. = 
uet er Tin, (5.2.22) 


f,0 97,0 
5.2.4 Semi-circular and semi-elliptical surface cracks as 


defects 


The derivation of strength expressions for other crack shapes is identical as has been 
outlined for edge cracks. The geometric function Ab/a) has of course to be replaced by 
those for the special type of crack. For the semi-circular crack it holds at the deepest 


point (A) and the surface points (B): 


Kaas z= dao, F, t) (5.2.23) 
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FF) -1.17am 0600 ie + a2] (5.2.24) 
j ao A, 

F,(2)= tan 327 je m a (5.2.25) 
° A, A, 


This solution for Kas was given in [5.10] and has been used in Section 3.2. Since the 
stress intensity factors at the deepest point and the surface points are different it is 
impossible that the crack will maintain its shape during extension. The equilibrium shape 


will become a semi-ellipse with an aspect ratio of a/cz1 (c-half crack width). 


5.3 Strength predictions compared with experimental 
data 


As an example of application eqs.(5.2.19-5.2.22), a result from literature may be disc- 
ussed. Li and Tomozawa [5.11] investigated the effect of soaking time on the dynamic 
fatigue behavior of silica glass for soaking at 250 °C in saturated water vapour. All their 
strength measurements were performed at room temperature, in air at a relative humidity 
of 15% to 30%. The material tested was a grade of industrial silica, TO8, by Heraeus- 
Amersil Inc. The strength data after different soaking times are re-plotted in Fig 5.1a. The 
dynamic strength measurements are shown in Fig. 5.1b in form of normalized strength 
increments as functions of exposure time / (squares in Fig. 5.1b). The curve is tentatively 
introduced as a guideline to the squares representing a soaking-time dependency of £^ 


since Koc Vb at". 


A further interesting result found by Li and Tomozawa [5.11] was the observation that 
the crack growth exponent of the v- Ki. curve increased with exposure time. The circles in 
Fig. 5.1c show the reported exponents as a function of the strength increase for a stress 
rate of do,pı/dt= 10 MPa/s. The perpendicular bars represent the standard deviations. The 
straight line represents the predictions by eqs.(5.2.19-5.2.22). For comparison, we have to 
take into account the large data scatter in Fig. 5.1c, and, consequently, a large uncertainty 
of the experimental m’-values. For a comparison of measured and computed z-values in 


Fig. 5.1, we determined the confidence intervals of the values 1/(7’+1) from Fig. 5.1a. 
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Fig. 5.1 a) Dynamic bending strength tests by Li and Tomozawa [5.11] (for clarity the symbols and scatter 
bars show slightly displaced abscissa values), b) normalized increase of dynamic strength versus 
heat-treatment time at 250?C; Squares: results from [5.11], curve: tentatively introduced guide line 
to the squares representing a dependency ocv b«c #4, c) n’-values reported in [5.11] (circles), 
straight line: eq.(5.2.22), b) and c) at doapp\/dt= 10 MPa/s. 
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From the limits of these intervals we computed the related limits for m. Whereas the 
intervals for 1/(7’+1) are symmetric, those for n become strongly non-symmetric. The 
50%- and 90%-confidence intervals (CI) are compiled in Table 5.1. These intervals are 
also plotted in Fig. 5.1c. Li and Tomozawa [5.11] reported for soaking times of 0 and 4 
days the same interval as we found for 50% confidence levels. Confidence levels of 5096 


are not very often used for data assessment. For this purpose, mostly 90% or 95% levels 


are applied. 
Time | or(MPa) | ox-or (MPa) | r 50% CI 90% CI 
Od 59.7 0 32.7 | [29.7, 37.9] | [21.8, 63.7] 
ld 149.3 89.6 40.4 | [33.4, 50.9] | [23.06, 147] 
3d 183.4 123.7 80.4 | [64.2, 107] | [42.1, >1000] 
4d 207.0 147.3 9] | [75.1, 116] [51.4, 380] 


Table 5.1 Evaluation of the strength data from Fig. 5.1a with confidence intervals for crack-growth exponents. 


On the basis of the 90%-CI, we conclude from Fig. 5.1c that the theoretical prediction 
and the experimental results by Li and Tomozawa [5.11] are in agreement, i.e. the crack 
growth exponent increases as expected from eq.(5.2.22). Since the geometric function is 
Rbla)<l, the swelling stress at the surface must be |o.,|2147.3 MPa. Figure 5.2 shows the 
predictions introduced into Fig. 5.1c. This representation states that the predicted straight- 
lines according to eqs.(5.2.20) and (5.2.22) are in good agreement with the measurements. 
This can be seen especially for the strengths after 1 day soaking because in this case the 


standard deviations were given in [5.11]. 


An equivalent representation of the agreement between experiment and prediction is 
given in Fig. 5.3 where the measured strengths are plotted again with a linear ordinate 
scaling. The black curve represents the power-law equation (5.2.5) with m=32.7. In a 
linear plot this curve is of course no longer a straight line. The small curvature is hardly 
visible. The red curves are identical with the black one but shifted by constant stress 
values Ao; (indicated by the arrow in Fig. 5.3). From this plot, it can be concluded that 
the shifted curves agree with the measurements within the natural data scatter (standard 
deviation bars given for 3 days soaking exclusively). This fact again confirms the de- 
pendency (5.2.18) 
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>0,=0,,+const. 
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Fig. 5.2 Predicted straight lines (red) according to eqs.(5.2.20) and (5.2.22) introduced by the red lines. 


The derivations in this section were performed for the power-law description by eq.(5.2.5) 
because of its simple mathematical handling. From chemical reaction theory (see e.g. 
[5.12]) it follows for the v( K)-relation (see also Section 10.3) 


V(K) =v, exp[B K] (5.3.1) 


with material parameters vo and B. A first estimation of the parameters for eq.(5.3.1) can 
be obtained analytically when a crack-growth law is given as a power law relation and 


the exponent is large 


v(K)- AK" (5.3.2) 
From (5.3.1) we get the derivative 
dlnv es dlnv 1 = (5.3.3) 
dK +dK K 
= 
d(InK) 
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dlnv 
- BK (5.3.4) 
dink 
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Fig. 5.3 Data of Fig. 5.2 represented with linear ordinate scaling (note that curves are no longer exactly 
straight). Red curves agree with the black one shifted by a constant stress Aor. 
On the other hand we obtain from eq.(5.3.2) 
div _, (5.3.5) 
dlnK 
ie. BKzn (5.3.6) 
and 
v, = AK" exp[- BK] (5.3.7) 


When evaluating experimental results, a mean value of K on the v-K curve has to be 


used. As a consequence of (5.3.6) the effective slope 7? is equivalent to a coefficient B’ 
B'K =n' (5.3.8) 


Summary on strength results: The present section dealt with the influence of swelling and 


shielding in silica caused by the reaction between SiO» and water. 
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5.3 Strength predictions compared with experimental data 


We considered the effect of swelling stresses on the crack-growth exponent of the power- 
law for subcritical crack growth. It could be shown that the apparent exponent zx of a 
power law description according to 

voKk 


n 
appl 


must decrease with increasing swelling. This theoretical consequence of swelling in silica 


is in rather good agreement with experiments from literature. 
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6 Swelling anisotropy 


6.1 Motivation for anisotropic swelling 


Strongly different stress states can result in the same hydrostatic stress term on. An ex- 
ample for this fact is illustrated in Fig. 6.1 where in all cases the hydrostatic stress dis- 


appears completely, 05-0. This is trivial when all stress components disappear, Fig. 6.1a. 


In the case of Fig. 6.1b, where the transversal stresses are half of the longitudinal stress 
with inverse sign, it is again on=0 according to on-(o0x*0y*o;)/3. The same holds for 


torsion loading, 6,—-o;, Fig. 6.1b’. 
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Fig. 6.1 Three stress states acting on a Si-O-bond exhibiting the same disappearing hydrostatic stress on=0: 
a) isotropic swelling is expected, for b), b’) swelling in z-direction should dominate as illustrated 
by c), case b’) describes torsion loading, c) water molecule reacting with a stressed bond oriented 


in zdirection, producing two hydroxyls, 2.5:. 
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6 Swelling anisotropy 


It is hardly imaginable that the high tensile stresses in bond direction should have no 
effect on bond fracture. It has to be expected that the reaction in tensile direction should 


produce swelling even for the hydrostatic stress state on=0. 


Figure 6.2 shows the load-displacement curves for bond breaking by reaction (1.1.1). 
Perpendicular to the loading direction the necessary energy for splitting the bond is U 


and must be delivered by thermal fluctuations. 


b) F 


mn. 


Fig. 6.2 Energy conditions for bond breaking. 


U-Fappi 04 


The portion M (subscript x means the bonds oriented in x-direction in which no stress 


acts) of opened bonds results from the Boltzmann equation as 


RS erf- =| (6.1.1) 


The same equation holds for the y-direction 
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6.2 Anisotropic swelling in tension and torsion test 


U 
N = (6.1.2) 
ye exo 


The potential barrier U belongs to the case when no stresses are present or when the bond 
is perpendicular to the stress, Fig. 6.2a. Under action of a force Fin bond-direction, Fig. 
6.2b, the potential barrier is reduced to U-FS:. Consequently, the occurrence of opened 


bonds in z-direction is increased 


N,« ef- er (6.1.3) 


where Vis an activation volume. In the following considerations, it will be shown that 


the strongest anisotropy effect has to be expected for the torsion loading case. 


6.2  Anisotropic swelling in tension and torsion test 


6.2.1 Uniaxial loading 


The case of uni-axial loading was outlined in detail in [6.1] and Sections 6 and 7 of [6.2]. 
In the following considerations, we will handle also the torsion test by superimposing 


uniaxial stresses in y- and z-directions as illustrated in Fig.6.3. 


Fig. 6.3 Loading cases: a) uniaxial tension, b) uniaxial compression, c) cylinder under a torsion moment 
M. 
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6 Swelling anisotropy 


The hydroxyl concentration .S is under purely uniaxial stress according to eqs.(1.1.15, 
1.1.16) 


S=S, exp | (6.2.1) 


where A Ver is an effective activation volume of 


AV,, 214.4 cm "/mol +10% (6.2.2) 


eff 
and So the initial hydroxyl concentration in the absence of load. 


For small arguments of the exponential function, eq.(6.2.1) can be linearized by first 


order expansion 


AV, AV. 
ses. zi u - sus Lu (6.2.3) 


The increase of the hydroxyl concentration by the stress, A.S=,S—So, is simply 


26 ef (6.2.4) 


Taking into account the proportionality between S and the volume swelling strain £v 


according to eq.(1.1.9) 
&, -KXS,k = 0.97 (6.2.5) 


eq.(6.2.4) yields for the stress-affected increase As, —£v—£vo = KAS 


op Ver (6.2.6) 


appl 


AE, = Ev RT 


The individual swelling strain components were found to be interrelated by [5.1] 


O 2g? = 0.2469 (6.2.7) 


SW,X swy -— SW,Z 


€ 


where the upper index (Z) stands for "loading in z-direction". Equation (6.2.7) gives rise 


for the definition of a “Poisson’s ratio Vsw for swelling strains 
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6.2 Anisotropic swelling in tension and torsion test 


_ Ewy (6.2.8) 


Its value would result in vsw=0.24. Finally, the stress-affected part of volume swelling is 


obtained as 


Ae, 2&0 +EP +8 20,5260 (6.2.9) 


v SW,X SW,y SW,Z SW,Z 


so that from eqs.(6.2.6) and (6.2.9): 


E2, =1,926,, Canh Var (6.2.10) 


SW,Z 


6.2.2 Torsion loading 


The results for the uniaxial loading case may be used to describe more complicated load- 
ing by superposition. Under torsion loading, the maximum and minimum principle 
stresses appear in a coordinate system turned by 45° with respect to the length axis. In 


this system, Fig. 6.3, the “applied” stresses are 


[oj Oo o =0 (6.2.11) 


appl,y = appl,z ? appl,x 


In order to show the principal effect of torsion stresses on the swelling strains and hy- 
droxyl concentration, let us assume a volume element without restrictions on the boun- 
daries, so that swelling strains can freely expand or contract without generation of swell- 
ing stresses. In the absence of load, this volume element may contain an initial hydroxyl 


concentration So. 


For the applied stress in z-direction, the swelling strains are given by eq.(6.2.7). Under 


the compressive stress in y-direction, the subscripts have to be changed 


MN 2 y gO =y gO (6.2.12) 


SW,Z SW,X sw sw, y sw” sw,z 


€ 


The total swelling strains from (6.2.7) and (6.2.12) then read by superposition: 


E, 7 £O, c eO, =(1+v,)e® (6.2.13) 


SW,Z SW,Z SW,Z SW,Z 
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6 Swelling anisotropy 


as) [y = (z) 
Esw,y Esw,y H Esw,y (a H Vow)Eswz 
=e) o) _ 
E wx = E wx + E wx = 0 


The total change of volume swelling As, is simply 


AE, = Egy x + EV y + 6,,70 


SW,y SW,Z 


(6.2.14) 


(6.2.15) 


(6.2.16) 


As an important consequence, we have to conclude that under torsion loading, no volume 


increase occurs. From eq.(6.2.5) we have to conclude furthermore that no increase of the 


hydroxyl concentration can occur, i.e. AS=0 > S=So=const. 


The e/astic strains result from Hooke's law 


1 V 1-v 
£, CE? 5 > €, CE Ow 
m l+v 
€, "E? ——0, >E, =p Capt 
e€ =0 


The ‘otal strains as the sums of elastic and swelling contributions read 


€, total E E, at € oz 
€ y total = €, ud Esw,y 
E 0 


x,total ^ 


(6.2.17) 


(6.2.18) 


(6.2.19) 


(6.2.20) 


(6.2.21) 


(6.2.22) 


A small volume element A V is considered, that contributes to the torsion moment or the 


drilling angle. The share to the applied moment may be denoted as AM,appi, the contribu- 


tion by swelling as AMisw. Since the swelling strains have the same signs as the elastic 


strains, the deformations in a torsion test on a cylindrical specimen under constant load 


must increase by swelling. In the first case, the torsion angle, increases by an amount 


AQtsw 
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6.2 Anisotropic swelling in tension and torsion test 


AQ, sw - 2 = (1+ Vow) 1.928, E AV, (6.2.23) 
AQ, appl £, (19 v) RT 


Inversely, the stresses must decrease due to swelling, if the cylindrical specimen is load- 
ed under displacement conditions where the drilling angle, @, is kept constant. For the 
latter case, the ratio of the torsion moment by swelling, Ay, and due to the applied 


moment, AMspı, are simply related by 


AM ww __ AD (6.2.24) 


6.2.3 Expectations for torsion tests 


It is self-evident that only the water-affected zones directly under the surface can con- 
tribute to the change of the torsion moment or twist angle, i.e. water must be present at 
the considered location. The water diffusion layer b increases with time ¢ according to 
eq.(3.1.2). 


M=const 


time t 


Fig. 6.4 Time dependence of twist angle (red curve) in a constant moment test and time dependent mo- 


ment for twist angle kept constant (black curve). 
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6 Swelling anisotropy 


The main expectations are: 

1. the remaining S after the test must be the same as the value So before the test 

2. ina test under constant twist angle (black curve in Fig. 6.4), the swelling stresses in 
the surface layer must decrease, resulting in a reduced torsion moment M=M-AM 
with AMecéww,2X bct 

3. in a constant moment test (red curve in Fig. 6.4), the reduced moment AM must be 
compensated, i.e. an additional twist angle Age t must appear, so that the instanta- 
neously applied angle qo is increased to @otAQ. 


4. After unloading a remaining twist angle Ag should be observable. 


6.3 FE-Model for a single reaction event 


Due to the reaction (1.1.1), a SiO-bond may fracture in the absence of externally applied 
stresses. The coordinate system is oriented so that the zaxis agrees with the direction in 
which the repulsive forces between the silanols act (Fig. 6.5a). The related “nano-pore” 
defined by the dashed circle in Fig. 6.5a was modelled as a sphere of radius œ in an “in- 
finite body" [6.3]. 


y 
) a e. 
I 
a) YY 


es Oo 


Fig. 6.5 a) Bond breaking by water/silica reaction; the two hydroxyls yield an expansive displacement 8o in 
z-direction; the volume element in which the reaction occurred is represented by the dashed cir- 


cle, b) Finite Element model of a “nano pore" loaded by axial displacement ôo. 
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6.4 FE-Results 


In a FE-study, we realized the “infinite body" by a cylinder of radius 500 R and half 
height 500 R, Fig. 6.5b. For the computations, we used ABAQUS Version 6.9 on a mesh 
of 1246 elements and 3880 nodes. Due to symmetry, only a 4" of the total body had to be 
modelled. As the only material parameter of relevance we used the Poisson ratio of 


v=0.17. Expansion of the pore is prescribed by a displacement 8o in z-direction. 


6.4 FE-Results 


Figure 6.6 represents the displacement distribution for a cylinder contour with A=p=16 
R, so that the distance from the pore is sufficiently large. On the other hand, the contour 
is small compared with the modelled body. Saint-Venant's theorem ensures that no 
boundary effects and details of the pore shape can affect the results noticeably. The dis- 
placements are normalized on the maximum displacements present on the symmetry 


lines. 


The volume changes in axial and radial directions, A V; and A V, are given by 


p h-p 
AV, -2z[&,()rdr > AV, -2zp [6,(z)dz (6.4.1) 
0 0 
Since the cylinder volume is 
V, = p ah = pa (6.4.2) 
ph 
the increase of the volume swelling strain results in 
em ann (6.4.3) 
pu 
Finally, the anisotropy ratio £,, , / A€&, is obtained by: 
£w. | AV, (6.4.4) 


Ae, AV, +AV. 


For the cylinder contour of p = A= 16 R, the evaluation of eq.(6.4.4) results in 
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6 Swelling anisotropy 


Font 9 06 (6.4.5) 


AE, = Epyx Ey + Eas = Eus + Few, (6.4.6) 


a “Poisson’s ratio V ” for swelling strains can be defined via 


p= Emr (6.4.7) 


E 


z 


0.2 0.4 


l ı L 


Fig. 6.6 Displacement distributions 8; and ô at the cylinder contour in a distance of A=p=16 R. 


From eqs.(6.4.5) to (6.4.7) we obtain 


Emz — fona = E l — =2.06> 7 = 0.257 (6.4.8) 


resulting in a Poisson’s ratio as is usual for ceramics. 
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6.4 FE-Results 


The transverse swelling strains are 


Eswr _ Eswx  Üwy _ 0.53 (6.4.9) 


It has to be emphasized here once more, that the results from FE were obtained under the 
assumption that all stress-enhanced swelling strains are oriented in direction of the load. 


This is of course an upper limit case. 


Influence of the chosen geometry: In order to show the effect of the finite body that had 
to be modelled for the FE-computations, we computed the average strains on cylinder 


contours, € , in different distance from the pore. Figure 6.7 shows the ratios £,/ AE, and 


EE. 


= -E,/€, 
E, AE, E 
21r 
zk 4 
._—— 1 
M - 0.25 
19 r = 
02 
1.8 F 
1.7 l 1 l 1 ! i ! : l 1- 0.15 
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h/R=p/R 


Fig. 6.7 Strain ratios £, / A&, and — e. lE, vs. A/R or p/ R. 


It can be seen that the results are rather independent of the chosen cylinder contour. For 
E, I ^£, the maximum variation within the range of 168 < Ap < 100 R is less than 2%. 
This weak dependency clearly indicates that the theorem of Saint- Venant is sufficiently 
fulfilled, i.e. that 


63 


6 Swelling anisotropy 


R««h, p<<500R 


When this condition holds, the validity of Saint-Venant's theorem suggests that special 


details of the pore-shape cannot have a noticeable influence on 6,1 Ag, and Ele 
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7 | Stress-enhanced isotropic swelling 


7.1 Quasi-hydrostatic stress state at crack tips 


In gases and liquids, the loading is always hydrostatic. In a solid the situation is more 
complicated, since the individual stress components ox, Oy, Oz are in general independent 
of each other and are not necessarily equal. Especially in uniaxial tension or compression, 
the hydrostatic stress deviates clearly from the tensile stress. In contrast to this, the stress 
state ahead of crack tips is more hydrostatic since the tensile stress in the prospective 
plane and the tensile stress normal on this plane are identical, ox-0y, and the tensile stress 


6; is very close to the tip: 0,=v(0x+0,). 


In the case of purely hydrostatic loading 6.-20y-0;, swelling must be isotropic since no 
direction is preferred. In this Section the behaviour of cracks is considered and therefore 
isotropic swelling behaviour with identical swelling strains &swx=&sw,y=&swz may be as- 


sumed. 


In order to allow a most transparent analysis, we assume in the following considerations 


the applicability of: 
1. Jsotropic swelling behavior, €swx=€sw,y=Esw.z, independent of the stress direction. 
2. Replacement of the pressure p by the negative Aydrostatic stress term on. 


Whereas in uniaxial tension the ratio of the hydrostatic stress to the maximum stress (the 


stress oy normal to the crack plane) is only 
0,/0,=3 (7.1.1) 
it holds along a prospective crack plane 


[07 
=3dtv) z 078 (7.1.2) 


y v=0.17 
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7 Stress-enhanced isotropic swelling 


This value is much closer to the purely hydrostatic stress state described by on/oy=1 and 
should therefore be applied to crack problems. This possibility of a captivatingly trans- 
parent description will be chosen here accepting the probability of underestimation of 
swelling effects. In experiments in which uniaxial stresses rather than hydrostatic ones 
are applied during the reaction, this assumption was re-examined in Section 6 to account 


for the possibility that local swelling is not isotropic [7.1]. 


In the following Section, the experiment by Tomozawa et al. [7.2] will be discussed 


under the simplified assumptions mentioned before. 


7.2 Application of simple crack-tip models 


7.2.1 Sharp crack approach 


Different crack-tip models are commonly used in fracture mechanics for the description 
of stresses ahead of a crack. Linear-elastic fracture mechanics as a continuum mechanics 
theory assumes the existence of a mathematically sharp crack. In this case, the stresses 


near the crack tip become singular, 6>% for a crack-tip distance of 70. 


The stresses in the crack plane ahead the crack tip, caused by the total stress intensity 


factor Kip, are 
o, =2vo (7.2.1) 
with 


EET. (1.2.2) 


Generally, it is assumed that crack extension under inert conditions occurs, when the 
normal component oy equals the ideal strength oo. When Ki. is the fracture toughness, it 


then simply holds 


LET (7.2.3) 
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7.2 Application of simple crack-tip models 


The fracture toughness of silica is Ki. — 0.75-0.8 MPaVm [7.3] and the theoretical strength 
is about [7.4] 


o= = GPa (7.2.4) 
m 
as is in agreement with strengths up to 25GPa measured by Brambilla and Payne [7.5] on 


extremely thin silica fibers of about 60 nm radius. 


The highest tensile strengths for thicker silica glass fibers in ultra high vacuum at room 
temperature are about ©: = 12.6 GPa [7.6]. Strengths for optical fibers are lower than the 
theoretical strength, co, due to the unavoidable presence of surface defects in these fibers 
[7.7]. 


The critical distance r« at which the normal stress component o; for K=Kic equals the 


theoretical strength, is according to (7.2.2) 


2 
r= ze = 0.19 nm (7.2.5) 
24740, 


This value is smaller than the "microstructure" given by the SiO» ring diameter. This 
makes clear that failure within such small dimensions cannot be described by a mathe- 
matically sharp crack. Consequently, the K-description via continuum mechanics must 


fail, when the tip distance competes with the microstructure of the material. 


7.2.2 Slender notch approach 


In a micro-structurally motivated approach, the crack tip region (Fig. 7.1a) is considered 
as a slender notch with root radius p in the order of the average radius of the SiO» rings, 
Fig. 7.1b. This description may be applied in the following considerations. For such a 


notch, Wiederhorn et al. [7.8] suggest a crack-tip radius of p = 0.5 nm. 


Stresses at slender notches were given by Creager and Paris [7.9]. The stress component 


normal to the crack plane (the tangential stress o=oy) is in distance € from the notch root 


2K pts (7.2.6) 


t aeo 28) p +28 
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7 Stress-enhanced isotropic swelling 


The other stress components are 


o 2 22K 5 (7.2.7) 
*oXm-(o-24) p+ 26 
_ 2K (7.2.8) 


Az (p +25) 


Consequently, the hydrostatic stress term from the applied stresses results as 


2(1+v)K (7.2.9) 


The location at which failure occurs can be computed by application of the “Theory of 


[07 


Critical Distances” (TCD) as was developed by Taylor [7.10,7.11]. In a “non-local ap- 
proach”, the TCD postulates that failure occurs when a distance dependent effective 


stress Oet exceeds the tensile strength oo. 


The effective stress is in this approach the stress value oy in a distance of €=Z/2 from the 


notch root, where the length Z is given by 
2 
pub e daas (7.2.10) 
Z7 Os 


ie. twice of the critical crack-tip distance in the continuum mechanics description by 
eq.(7.2.5). Introducing (7.2.10) and (7.2.4) into (7.2.6) and setting €=Z/2 results in 


p = +5) 


2 


K? = 0.62nm (7.2.11) 
270, 


This says that in the light of the TCD-approach the guess by Wiederhorn et al. [7.8] is a 
very good estimate. Their value of p=0.5 nm would be exact for a theoretical strength of 
00=25.6 GPa. 


The hydrostatic stress at K=Kic is 


O appl,h x (v5 -1) n Op (7.2.12) 
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7.2 Application of simple crack-tip models 


For a test under subcritical crack growth conditions, the necessary stress intensity factor 
is smaller than fracture toughness Kic. From the average crack rate of væ 10°’m/s (estima- 
tion in [7.12]) and crack-growth curves in water by Wiederhorn and Bolz [7.13], the 
stress intensity factor results to be Kapp: = 0.48 MPavm. 


The total stress intensity factor is the sum of the applied, Kapp, and shielding term Ksh 


K=K HR Ky <0 (7.2.13) 


app 


The shielding stress intensity factor was estimated in [7.12] to be Ki, =—0.1 MPavm at v 
= 10°’m/s. Then the total stress intensity factor would be K=0.38 MPa‘\m. 


For the reduced stress intensity factor during subcritical crack growth, K<Kic, eq.(7.2.12) 


reads 
(+v) K 
O appl,h = (45-1) 30 (7.2.14) 
3 Ky 
Fig. 7.1 a) Crack in silica terminating in a nano-pore, b) equivalent slender notch with a finite notch root 


radius p, grey molecules are mechanically inactive. 


Swelling stresses under isotropic swelling conditions at slender notches were determined 
in [7.14] by using the FE-method. Results for the hydrostatic stress term under plane 
strain conditions are shown in Fig. 7.2 by the circles. It can be concluded from this repre- 
sentation that the hydrostatic swelling stress normalized on the volume strain hardly 


changes over the swelling zone. 
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7 Stress-enhanced isotropic swelling 


Outside the swelling zone, the hydrostatic stress term disappears practically. In addition, 
the effect of the swelling zone height 5 is small. The square indicates the theo-retically 


known value for an infinitely thin surface layer, given by 


Ee 
o, 2-25. ns (7.2.15) 
1-v 
with the coefficient 
20.97E 
=- =18.7GPa (7.2.16) 
9 1-v 
0.05r 
x/b 
0.5 
1 1L 
-0.05r 
b/p 
-0.1f e 0.2 
@ 0.5 
-0.15r © 1 
-0.2r 
-0.25r i 
E 
Fig. 7.2 Hydrostatic swelling stresses, normalized on the volume expansion strain for 3 different relative 


zone heights b/p. 


The hydroxyl swelling stress at the surface (x0) can be approximately expressed as 


af?) Ele c sera | (7.2.17) 
9(1—v) p 
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7.3 Hydroxyl concentration under swelling stresses 


The other individual swelling stress components in the directions of a rectangular coordi- 


nate system are given in Section A1.3. 


7.3 Hydroxyl concentration under swelling stresses 


Since now the hydroxyl concentrations S, So are known from Sections 5 to 8 in [7.15] 
and applied hydrostatic stress oappin from 7.2.2, we can estimate the molar volume. 
In presence of externally applied and internal swelling stresses, Oappı and Osw, the total 


stresses are given as their sum 


O — Opi + Oso Th = Opin + Os (7.3.1) 
This leads to 
S =S,exp (Gupta + Ponas (1.3.2) 
? RT 


The hydrostatic swelling stress is proportional to the hydroxyl concentration 
=-nS (7.3.3) 


with a coefficient of proportionality n, depending on the stress state of loading. Conse- 


quently: 


S -S, eo [Em SWer ) (7.3.4) 


This is an implicit equation with respect to S that occurs on the left side as well as in the 


exponent. Its explicit form is 


SAT nag 15% of J (7.3.52) 


where PLog stands for the product-logarithm function (see e.g. [7.16]). Since this func- 
tion may not be available on all computer systems, the inverse representation can be used 


resulting equivalently in 
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7 Stress-enhanced isotropic swelling 


2484 a(S] (7.3.5b) 


S 0 


O appl „h 


In terms of the hydrostatic swelling stress, eq.(7.2.15), we analogously obtain 


Ca = RT PLog Owh Vs exp O appl, h V's (7.3.6) 
V. RT RT 


where Osw.ho is the swelling stress for S=.. Having in mind that the stress state ahead of a 
crack is not exactly a hydrostatic one and swelling cannot be completely isotropic it 


would be better to call the result an effective volume V, >A V“, On the other hand, it 


has to be expected that A V* may depend on the chosen crack-tip model, too. 


The effective volume can be determined from the measurement of hydroxyl concentra- 


tion as 


aver 2. RT | S | (7.3.7) 
© appl,h = ns S, 


Unfortunately, the authors of [7.2] didn't provide correct stress intensity factors for their 
tests (details in the Appendix of [7.17]). The resulting A V“" is plotted in Fig. 7.3a versus 
the theoretical strength for Kappi 0.48 MPaVm according to [7.13] and constant shielding 
stress intensity factors of Ka-0 to 20.3 MPa\m. 


Figure 7.3b shows Ver for 3 applied stress intensity factors Kip for which the stress in- 
tensity factors Ki; were computed by use of eq.(10.6.6). For the most probable value of 
theoretical strength, oo = 23 GPa, large values of the effective volume AV“ have to be 


expected for the low-temperature region «450?C. 


At temperatures T> 500°C, the water/silica reaction is of the order 2 


=Si-O-Si= +H20 o 2[SiOH] (7.3.8) 


At temperatures T < 450°C, hydroxyls are immobile. Only those hydroxyls can react in 
the reverse reaction step, which are directly neighboring. Consequently, the reaction is of 
first order. When [2.5]-[2 SiOH] is the concentration of the immovable hydroxyl, the 


reaction equation reads now 
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-Si-O-Si- +H>0 © [2 SiOH] (7.3.9) 


The activation volumes related to the two completely different reactions may now be 


quite different. 


efff 
AV b) 
100 p 100 

i eff 
TH AV 

L 50 

(cmĉ/mol) + 
0.43 MPaym 
(cm?/mol) | 
10 L 0.46 MPavm 
10r 
5 L | 
5r. 
Ks, (MPavm) l Kappi=0.48 MPa vm 
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Oo (GPa) Oo (GPa) 
Fig. 7.3 Effective volume A V“ derived for the slender notch crack-tip model as a function of the ideal or 
theoretical strength oo of 00=23 GPa; a) for constant shielding stress intensity factors at 
Kıppı=0.48 MPa m, b) for constant applied stress intensity factors according to eq.(10.6.6). 
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8 Effects of swelling at crack tips 


8.1 Consequences of shielding on v-K-curves 


The shielding stress intensity factor, Ks, caused by an expanding zone with a constant 
volume strain over the zone size was given by McMeeking and Evans [8.1]. The shield- 
ing term for variable hydroxyl concentrations was described in [8.2] and will be extended 


to variable Young's modulus in Sections 11.6.1 and 12.2. Generally it holds 
Ky = [(Sx E)hto)do (8.1.1) 
0 


with the Green's function 


pean dec e (8.1.2) 


and y«0.22 [8.1]. The effect of the shielding stress intensity factor on the v-K-curves 
was outlined in [8.2] for constant modulus Æ. Equations (8.1.1) and (8.1.2) result in the 


same relation as given in [8.1] with an effective coefficient yer 


E&E 
Ky, =Y ey Do (8.1.3) 


From theoretical considerations by Wiederhorn [8.3], the subcritical crack growth law is 
of the form 


def 


v 2 v, exp(bK,, / RT) - v, exp(BK yp) (8.1.4) 


tip 
defining B=b/RT. Equation (8.1.4) would suggest for small stress intensity factors that a 
subcritical crack growth rate of vo would occur in the absence of any loading. The model 
of single bond rupture by Michalske and Freiman [8.4] gives rise for the back reaction. 


Instead of (8.1.4) it was suggested in literature (see also [8.5]) to use 
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v -2v,sinh(bK,, / RT) (8.1.5) 


The effective stress intensity factor acting at the crack tip, Kip, represents the singular 


stress field at the crack tip. From the principle of superposition, Ki, is given by 


K., =K (8.1.6) 


tip appl 


+K 


sh 


where Kappi is the externally applied stress intensity factor and Kas is the stress intensity 


factor due to shielding. 


The effect of a shielding stress intensity factor on the v- K-curve of silica is a shift of the 
V- Kip-curve to the left (to lower stress intensity factors) as has been shown in [8.2] on the 
basis of the v-Kapp-curve from Double Cantilever Beam (DCB)-measurements by 
Wiederhorn and Bolz [8.6]. 


8.2 Threshold behaviour under shielding conditions 


8.2.1 Approximation by rate-independent shielding stress in- 
tensity factor 


In the case of silica, swelling causes a negative shielding stress intensity factor. This fact 
gives rise for a threshold of the v-Kappi-curve as was shown in [8.7]. This can be simply 
concluded from eq.(8.1.5). In order to give a transparent description of the general con- 
sequences of shielding first a constant shielding stress intensity factor is assumed. For 
this purpose a shielding stress intensity factor of Ku=—'4 Kic was chosen (note that — Ki is 
a positive number). In Fig. 8.1 the dash-dotted curve represents the v- Kip dependency as 
the subcritical crack growth curve of the material in the absence of shielding, Ka/Ki.—0. 
Compared to this curve, the shielding-affected blue curve for Ksn=—!'4 Kic is shifted to the 
right. A threshold value at 1⁄4 Ki. is visible. 


From the data of Fig. 2.8b in [8.8] it can be concluded that not exclusively swelling of 
glasses can be expected. Whereas the Al-glass and the Ca-containing glass of Fig. 2.8b 
should show negative shielding stress intensity factors due to volume expansion similar 
to silica, the volume shrinking of the B20s-glass and the K-containing glass must result in 


positive shielding stress intensity factors, where now the notation "shielding" does loss 
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8.2 Threshold behaviour under shielding conditions 


its original sense. The related v- Kiyp-curve for Ka-!4 Kic (red) shows a finite subcritical 
crack growth rate that does not disappear at Kapp=0, no surprising result in the light of 


shrinking. 


Fig. 8.1 Effect of constant shielding on the v- K curves for negative, disappearing, and positive shielding 


stress intensity factors. 


8.2.2 Change of the swelling zone shape during crack growth 


Muraoka and Abe [8.9] carried out static tensile tests on silica fibers with small semi- 
elliptical surface cracks. By microscopic inspection, the crack sizes were determined and 
via time measurements, the crack growth rates were computed. In Fig. 8.2a, the crack 
growth rates are plotted versus Kap (circles and black curve) after correction of the stress 
intensity factor solution from v=0.3 to v=0.17 [8.10]. The blue line represents results by 
Wiederhorn and Bolz [8.6] and is introduced as a reference curve. In terms of the applied 
stress intensity factor, Kapp, this curve can be expressed in the range of 3x10! 'm/s« v«10 
?m/s by 


v =v, exp(B'K (8.2.1) 


appl ) 


with the parameters vo 74.5 10? m/s, B" =91.2/MPavm for room temperature. 
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Fig. 8.2 a) v-K-curve of subcritical crack growth for silica obtained from static lifetime tests by Muraoka 
and Abe [8.9] (circles and black curve), blue line: results by Wiederhorn and Bolz [8.6], red data 
and line: re-evaluated data by Sglavo and Green on cone cracks [8.11], b) Development of the 


swelling zone and the applied stress intensity factor for tests with constant load (schematic). 


Figure 8.2a shows v-K-data on cone cracks by Sglavo and Green [8.11] as the red 
squares. The absolute values of original stress intensity factors indirectly obtained from 
fracture toughness Ki. were corrected by a Finite Element solution including partial crack 
closure and crack-face friction as derived in [8.12]. In this case, a threshold stress intensi- 
ty factor at about 0.35 MPavm is visible. From Fig. 8.2a it can be concluded that the 
measurable crack-growth curves WKappi) are not a material property but are also influ- 
enced by the testing procedure. When a crack starts growing in a test under static load, 
the initial zone shape is first heart-shaped. Such a zone does not show a shielding stress 
intensity factor as has been early stated by McMeeking and Evans [8.1]. Consequently, 
the v-Kappi-curve must start on the curve for Ki-0, i.e. y=0, indicated in the schematic 
Fig. 8.2b by the circle. With increasing crack length, the shielding zone develops in 
length direction [8.1] resulting in an increasing coefficient y until the saturation value of 
about w=0.22 is reached. Then any further crack growth goes along the red curve, which 


may show the result from a crack-growth test under monotonously decreasing stress 
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intensity factor as is for instance feasible in tests on cone cracks introduced by a Vickers 
indentation. Such a crack shows a crack extension Aa that is large compared to the height 
œ of the swelling zone, A2»50, so that the shielding zone is fully developed. As has been 
shown by McMeeking and Evans [8.1] the shielding coefficient is not a constant but 


needs a crack extension of at about Aa250 before w=0.22 is sufficiently reached. 


8.2.3 Threshold stress intensity factor 


Figure 8.3 shows W(Kappi)-results for silica in water at room temperature by Sglavo et al. 
[8.13] (solid circles) together with the results by Wiederhorn and Bolz [8.6] (dash-dotted 
line). The threshold value for subcritical crack-growth was obtained in [8.13] as 
Ka=0.31-0.32 MPavm. The open circles represent the measurements of subcritical crack 
growth at 700°C in lab air of relative humidity 40-60 % by Aaldenberg and Lezzi [8.14]. 


Eq.(8.2.2) 


107r [8.11, 8.13] 


2 3 
1000/T (°K*) 
PEEN SN CS es Dort 


0.6 0.7 0.8 
K  (MPavm) 


Fig. 8.3 Subcritical crack growth data at room temperature obtained from the ISF-procedure by Sglavo 
and Green [8.11, 8.13] (solid circles), dash-dotted line: Wiederhorn and Bolz [8.6], measure- 
ments by Aaldenberg and Lezzi [8.14] at 700°C in lab air (open circles). Inset: Threshold vs. 


temperature. 
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Fitting these results by equations (8.1.5) and (8.1.6) with Ki —- Kw results in: 


For 20°C: w=2.5x10" m/s, B=90 (MPa\m)", Ka =0.32 MPaYm, 


and 700°C: w=2x10° m/s, B=46.8 (MPaVm)"!, Kin=0.53 MPaVm. 


The related fitting curves are introduced in Fig. 8.3. Ignoring the different humidity in 
[8.11] and [8.14], the threshold stress intensity factor as a function of temperature 0 


(plotted in the inset in Fig. 8.3) may be described by an Arrhenius dependency 


K, = Ks Zr (8.2.2) 


with the parameters Kino=0.638 MPa vm and Q-1.83 kJ/mol. 


8.2.4 Variable shielding 


Deviations from the assumption of a constant shielding stress intensity factor Ki have to 
be expected. The steepness of the v- Kapp curve is influenced by two opposite effects. An 
increase Of Kapp! causes an increase of the shielding-stress intensity factor due to in- 
creased S-concentration. On the other hand, the increase in the crack velocity v greatly 
shortens the time available for diffusion in the crack tip region and reduces the size of the 


shielding zone. It remains unclear which influence dominates. 


1) According to the derivation by Wiederhorn, the stress intensity factor K is the physi- 
cally active total stress intensity factor Kip governing the singular crack-tip stress field. In 


a Ig( V)- Kup-plot this must result in a straight line: 
log(V/V,) =a, +b, Kọ (8.2.3) 


2) As the subcritical crack growth curve for silica in liquid water at 25°C (Wiederhorn 
and Bolz [8.6]) shows, the crack velocity as a function of the applied stress intensity 
factor is also very straight over a large region of crack rates, Fig. 8.4a. This can be ex- 


pressed as 


log(wv,) =a, +b, K (8.2.4) 


appl 
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Fig. 8.4 a) v-Kappi curve of silica in liquid water by Wiederhorn and Bolz [8.6] at a temperature of 25°C, 
b) comparison of v-Küp and v-Kappı —curves for a constant shielding stress intensity factor, Ksn 


=const., c) comparison of v- Kip and v- Kapp! -curves for variable Ksn(Kic)=0, d) most general case. 
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Equation (8.1.6) requires that even the shielding stress intensity factor must fulfil a simi- 


lar relation, namely, 


log(V/V,) =a, +b, K (8.2.5) 


If Ki depends on v, there must exist an intersection of the straight lines according to 
(8.2.3) and (8.2.4). Figure 8.4b shows the case of a constant shielding stress intensity 


factor 


K,, *f(v)--K, (8.2.6) 


When Ki, reaches fracture toughness Krc, the crack extension becomes very fast. Veloci- 
ties in the order of sound velocity occur. Consequently, diffusion effects of surface water 
into the crack-tip region will disappear as well as the related shielding term, i.e. 
Ka Ki.)0. 


This limit situation is indicated by the circle in Fig. 8.4c. In this special case, the 


shielding stress intensity factor holds in terms of Küp 
K, 
K,, ® i ES a hk (8.2.7) 


In terms of Kapp 


Kys- M K, (8.2.8) 
I Ki i 
and 
"E Kop Ku K, (8.2.9) 
K,, = Ki, 


as illustrated in Fig. 8.4c. The most general case is illustrated in Fig. 8.4d. In this case, 


the apparent fracture toughness is in the range 


K.<K 


c appl,c 


€ K, - K,, (8.2.10) 
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8.3 Swelling strains at the tip of an arrested crack 


8.3 Swelling strains at the tip of an arrested crack 


8.3.1 AFM results on surface overlapping 


Atomic force microscope (AFM) measurements on fracture surfaces from [8.15] were 
used to estimate the swelling strains [8.7]. From Fig. 8.5a one can see that there is extra 
material (the dark grey area) at the tip of an arrested crack held in water of room temper- 
ature at Ki = 0.25 MPa:m"? for 80 days. This material was not there before the crack was 


formed in the glass. 


Not any crack growth was observed during the long holding time. For causing crack 
growth, the stress intensity factor had to be increased above K= 0.37-0.4 MPavm. If the 
crack were just loaded and unloaded, the crack profile would be relatively flat and to 


experimental precision, the two fracture surfaces would overlap perfectly. 


The excess material ahead of the actual crack tip can be interpreted as a water diffusion 
zone where swelling by the volumetric strain & took place. After complete cracking of 
the specimen, each half of this region can expand normally to the new crack face. This 
results in a volume expansion and in an overlapping of the surface profiles. Figure 8.5b 
gives a few measurements of the overlapping displacements with a maximum value of 
about 6 nm. Our very first evaluation of the displacement results was performed based on 


the few equidistant displacements (squares introduced in Fig. 8.5b). 


Comparing these few displacements with finite element results provided an average vol- 
ume swelling strain of &,=13.8% [8.2, 8.7]. The related average hydroxyl concentration is 


Ey 


S= =14.2 wt% (8.3.1) 
0.97 


and the average of Young’s modulus results from Section 11 as 
E =5900 MPa (8.3.2) 


This is a strong reduction to only ~8% compared with the undamaged modulus of 
E9-72000 MPa. 


83 


8 Effects of swelling at crack tips 


0 1.00 


0 kK 1 1 4 i 
200 100 0 
X (nm) 


Fig. 8.5 Determination of glass expansion due to water penetration of a crack tip: a) AFM scan of a crack 
tip that was held under load in water at room temperature, Kapp = 0.254 MPa-m!” for 80 days, b) 
profile of the swelling region on the fracture surface, c) Finite Element analysis of a heart-shaped 


swelling zone [8.7]. 


8.3.2 Influence of reduced Young's modulus due to damage 


In our opinion the swelling visible from the AFM measurements is caused by stress- 
enhanced generation of silanol groups in the water-silica reaction, eq.(1.1.1) with the 
equilibrium of the reaction under tensile stress shifted to the right. The increased hydrox- 
yl content results in the observed volume. In order to check the influence of a probably 


reduced Young's modulus in the swelling zone on the displacements, the displacements 


were additionally computed for the case that in the swelling zone the modulus is reduced 
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8.3 Swelling strains at the tip of an arrested crack 


to 10% of that outside the swelling zone. Figure 8.6 shows a comparison of the related 
displacements. For reasons of clearness the results are slightly shifted in height direction. 
Differences between the two profiles are negligible. This makes clear that the evaluation 


of the swelling strains is not affected by damage accompanied by module reduction. 
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Fig. 8.6 Comparison of the displacements v; for a strongly reduced Young's modulus in the swelling zone 
(blue symbols) with the results for constant modulus (data points shifted slightly in height direc- 
tion for better visibility). 
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9.  Crack-tip shielding 


9.1 Effect of local swelling at surface cracks 


The stresses ahead a crack tip are according to eq. (7.1.2) nearly hydrostatic. Conse- 


quently, we obtain with the hydroxyl concentration 5 for on=0 


te (9.1.1) 
RT 


where œ is the universal gas constant, T the temperature in °K, and AV is the reaction 
volume. This relation implies that in the high crack-tip stress field nearly all water Cy is 
present in form of hydroxyl S. As could be deduced from density measurements, the 


hydroxyl generation results in a volume expansion £v according to eq.(1.1.9). 


vue 


c 
a 


Fig. 9.1 a) Swelling zone at the tip of an arrested crack under mechanical loading, caused by stress- 


enhanced diffusion, b) zone for a crack grown by Aa. 


For a sharp crack, the singular hydrostatic near-tip stresses are given as 


6, 2149) 5 cos(o/2) (9.1.2) 


J2zr 


where rand @ are the polar coordinates with the origin at the crack tip, Fig. 9.la. As a 


consequence of eq.(9.1.1), very high hydroxyl concentrations and swelling strains must 
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occur in the crack-tip region. This allows the assumption that near a crack tip the water is 


present predominantly in the form of hydroxyl water. 


9.2  Shielding stress intensity factors from measured 
water profiles 


9.2.1 Shielding zone at a crack tip 


Due to restrictions in free expansion, the swelling strains £y result in swelling stresses 
which give rise for the “shielding” stress intensity factor Ks<O. At first loading, the 
swelling zone at the crack tip z(01,P) is heart-shaped as illustrated in Fig. 9.1a. 


The height of the contour for constant hydrostatic stress on is due to eq.(9.1.2) 


PX 2 (%) (9.2.1) 
4/37 On, 


The shielding stress intensity factor caused by a swelling zone of height o is [9.1] 


K, --v nE Jo (9.2.2) 


where Eis Young's modulus and v Poisson’s ratio. The zone of an arrested crack (not 
grown before loading) is heart-shaped, Fig. 9.1a, resulting in the coefficient y=0. When 
a crack has grown at least for Aa = 50, Fig. 9.1b, a value of y-0.22 is reached as was 
shown by McMeeking and Evans [9.1] and confirmed in [9.2]. The coefficient y was 
computed in [9.1] also for different shapes at the zone end as illustrated in Fig. 9.2a, (A)- 
(C). Figure 9.2b shows the related coefficients y for the general representation by 
eq.(9.2.2) as the solid circles. We interpolate these data by an interpolating curve through 
the data points that can be expressed simply, see Appendix A1, by 


"E ox bam 2 To ) (9.2.3) 
90 


As an application of this relation, the coefficient y may be estimated for a deviating zone 


shape. Computations of swelling zones are mostly carried out for diffusion in a motion- 
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9.2 Shielding stress intensity factors from measured water profiles 


less coordinate system with a fixed crack tip, i.e. for an arrested crack. Whereas for a 
growing crack the transversal diffusion normal to the crack plane is hardly affected by 
the moving crack tip, the diffusion and the crack propagation compete in crack direction 
[9.3, 9.4]. In [9.5] it was found for crack rates v>10" m/s under the assumption that the 


diffusion normal to the crack-plane direction is not affected by the crack rate ( v,=0) 


8 n 8 


— te 
ADDE" @ AO 


(9.2.4) 


For n/o=0.8x8/V27=1.23, it results from eqs.(9.2.2, 9.2.3) a coefficient w=0.23, 


introduced in Fig. 9.2b as the open circle. 


(A) (d 


(B) 


(C) 


(D) 


—- fo <+ 


Fig. 9.2 Cracks grown for Aa>5o, a) swelling zones with different shape at the zone end, b) shielding 
stress intensity factors by eq.(9. 2.2) (circles) and interpolating curve according to eq.(9.2.3) for 
the zone ends given in a), open circle represents the case of a growing crack with competing rates 


of crack rate and diffusion, case (D). 
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9.2.2 Shielding stress intensity factor for varying hydroxyl 
concentration 


In Section 9.2.1 we considered zones with constant swelling volume ev. In the general 
case, we have to expect hydroxyl concentrations and swelling strains that decrease with 
distance from the crack. A Green's function procedure for such cases had been developed 
in [9.6]. Equation (9.2.2) holds for the case of a step-shaped constant strain &=& inside 
and &-0 outside the swelling zone. For the more general case, we have to subdivide the 


zone in parts of thickness dw’ at a distance, œ, from the crack plane (Fig. 9.3). 


Ey C) 
do’ o 
Fig. 9.3 a) Zone with constant Sand ev, b) differential layer of thickness dw’, c) continuously varying 


swelling strain. 


The stress intensity factor caused by this zone of strain &(@’) results from (9.2.2) as 


E 1 
d(K 10,222 do' (9.2.5) 
( sh) 2 b-vAo 


Now the stress intensity factor for the varying £v results from 


© 


K, - [ e, hio) do (9.2.6) 
0 
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9.3 Water profiles below crack faces of subcritically grown cracks 


with the Green's- or “weight function" A 


jsp t i. (9.2.7) 
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Fig.9.4 Water profiles measured by Lechenault et al. [9.7] on DCDC specimens fractured in heavy water, 
a) water profiles expressed by the neutron reflectivity b, b) normalized on the concentration at the 


surface. 


9.3 Water profiles below crack faces of subcritically 
grown cracks 


9.3.1 Results by Lechenault et al. [9.7] 


Subcritical crack growth tests were carried out on DCDC specimens of silica by 
Lechenault et al. [9.7]. The water entrance into the fracture surfaces formed by the pas- 
sage of the crack exposed to deuterium oxide D2O was evaluated with a neutron reflec- 
tion technique to measure the penetration of the deuterium oxide into the silica glass. The 
authors found a satisfactory fit to the reflection data by assuming that the water concen- 


tration was constant at the surface up to a distance of L, followed by an exponential de- 
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crease in concentration for distances greater than Z. Consequently, Lechenault et al. [9.7] 


fitted their results by the expression 


1 for z<L 
exp[-(z-L)/n] for z>L 


cuco] (9.3.1) 


and found the parameters Z=4.3 nm, n=3.5 nm for a region of low crack rates (v~10%* 
m/s) at K= 0.61 MPa-m"? and L-4.6 nm, n=2.3 nm for higher crack rates ( v~4x10° m/s) 
at K= 0.77 MPa-m!”. These water profiles are illustrated in Fig. 9.4. 


9.3.2 Result by Tomozawa et al. [9.8] 


The concentration of water below the crack surfaces was measured by Tomozawa et al. 


[9.8] using the nuclear reaction analysis (NRA). The profile of the hydrogen concen- 


tration is shown in Fig. 9.5a. In order to allow an approximate extrapolation of measured 


data to the surface, we fitted the data cu by a complementary error function 


Cu = Cap erfd Z |+ 0.005x10” / cm? (9.3.2) 


Background 


taking into account the background level of 0.005x10?/cm?. The coefficient cuo is the 


surface concentration. 


The red bar in Fig. 9.5a indicates the 90%-CI of the surface concentration. The depth bat 
which the concentration decreased to =50% of the surface value in Fig. 9.5a is about 5 
nm. By fitting (9.3.2) to the measurements, the best set of parameters was found in [9.9] 


to be 
ca.0=1.30 [1.055, 1.55] 10?/cnv), 524.95 [3.28, 6.61] nm (90% CD). 


with the mean squares sum of 11.65 (10?/cm?. The maximum hydroxyl concentration 


at 20, S(0), could be computed in [9.9] with the result 


S(0) 216.5[14.4, 18.4] (wt96) (9.3.3) 


92 


9.4 Computation of shielding stress intensity factors 


For a second description of the water distribution let us be guided by the results of 
Lechenault et al. [9.7]. We fitted the measurements by Tomozawa et al. [9.8] also via 


eq.(9.3.1) with the parameters 
cuo71.12 107/cm?, 73.0 nm, 122.96 nm 


and a mean squares sum of 11.78x(10?/cm?Y.. This result is shown in Fig. 9.5b. Having 
in mind the rather strong scatter of the surface values, the two representations in Figs. 


9.5a and 9.5b are equal. 


9.4 Computation of shielding stress intensity factors 


9.4.1 Evaluation for undamaged material 


Subtracting the background level and introducing eq.(9.3.2) into eq.(9.2.6) or (8.1.1), and 
assuming at the tip ,S/.Sozcu/cuo, yields 


ke [Sto hto)do'- -y421n LT br) (9.4.1) 
0 =y 


with the Euler Gamma function I that for the argument % reads I (3/4)=1.2254, so that 
we can finally write with /2/zT (3/4) = 0.9777 21 


Kay DE (9.4.2) 


1-v 


sh 


Comparing eqs.(9.2.2) and (9.4.2) allows to identify approximately: bzo. For the repre- 
sentation of the water profile by eq.(9.3.1), the evaluation of eq.(9.2.6) yields 


Ka = fesoho)do y SOF Wr + Wesen [Ee fE) (9.4.3) 


From the surface concentration in Fig. 9.5a, given by eq.(9.3.3), the shielding stress in- 


tensity factor results as 


K,, 2-0.21 MPa/m (9.4.4) 
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For the description by eq.(9.3.1) we obtain: 
K,, =-0.20 MPa-,/m (9.4.5) 


where =73% of the total value come from the constant part of the water profile. 
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Fig. 9.5 Hydrogen concentration in the surface region of a growing crack via NRA-measurements by 


Tomozawa et al. [9.8], a) representation by eq.(9.3.2), b) expressed via eq.(9.3.1). 


9.5  Crack-tip shielding and threshold behaviour 


When a grown crack has been unloaded and is again loaded, the crack tip will only see 
loading by a stress intensity factor if Kapp>-Ksn. Consequently, the shielding stress inten- 
sity factor causes a threshold value for the subcritical crack growth as is schematically 
illustrated in Fig 9.6b. The stress intensity at which Kota and crack growth rate disappear 
in a crack-growth test under monotonously decreasing loading may be denoted as 
Kapp Ktip=0)=Ko=or Kin where “th” stands for “threshold”. Then it holds for a crack-arrest 


test under decreasing load 


94 


9.5 Crack-tip shielding and threshold behaviour 


Korat = Kapp + Key 70 (9.5.1a) 


total T app! 


or K, =K,, =—-Ky, (9.5.1b) 
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Fig. 9.6 
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a) Shielding stress intensity factors for the unloading case of Ktotai = Kappit Ka <0, b) expected 


threshold values of the v-Kappi curve, schematic. 


It should be emphasized that this threshold behaviour can only occur if the cracks have 
already grown subcritically (Fig. 9.1b). Cracks that are loaded for the first time can grow 
immediately, since the shielding zone has a heart-shaped contour (Fig. 9.1a) for which 
Ka=0, [9.1]. The shielding stress intensity Kx is plotted in Fig. 9.1 for a zone of height œ 
according to McMeeking and Evans [9.1]. These data can be simply described by Evans 
and Faber [9.10] as 

Ky 


= 2 arctan(Aa / Q) (9.5.2) 
m 


oo 


The applied stress intensity factor for a constant applied load depends on the crack exten- 


sion Aa via 
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K appi = 1, 28 _ 1, 22/0 (9.5.3) 


appl,0 


where Kappio stands for Kappı,aa=0. 


d 
n 


Fig. 9.7 Shielding stress intensity factor Ks, normalized on its saturation value Ka, reached for Aa/@—00 
obtained by McMeeking and Evans [9.1]. 


For a crack growing at Kio 70, the applied stress intensity factor increases and the 
shielding stress intensity factor becomes stronger negative. The total stress intensity fac- 
tor must first increase and then decrease. After a minimum, the total stress intensity fac- 


tor rises again. 


The arrest condition is illustrated in Figs. 9.8a and 9.8b for short cracks of different ini- 
tial lengths a/ and varying saturation values X, of the shielding stress intensity factor. 
The shielding stress intensity factor can depend in principle on the externally applied 
load and the actual crack rate. For reasons of transparency we assume in Figs. 9.8a and 
9.8b constant values for o and K,. Crack arrest occurs in those cases for which the 
curves intersect the line Ktoa=0, indicated by the open circles. The dashed extensions are 
hypothetic total stress intensity factors which make no physical sense since stress singu- 
larities at a crack tip are only present for Ktoa>0. From Fig. 9.8b it becomes obvious that 
subcritical crack growth for small cracks can take place below that threshold stress inten- 
sity factor that occurs in the case of cracks under decreasing loading (here Ko-0.2 
MPavm). 


96 


9.5 Crack-tip shielding and threshold behaviour 


0.15 Kappı,o (MPavm) L 
Kota a) K total F b) 
tota - 
0.15} 
0.1 F 
(MPavm) (MPaymy 
0.05 o1r 
0.05- 
ZW Eus oe 8 Aao 
t \ = 
a —- 9075 i 
e CX Wi ie ; TE ii MEEN EENEN ETETA ETENA D 
EON [ 0.1 0.2 0.3 
L " m F 
iin NN 0.05 Kappı (MPavm) 


Fig. 9.8 Total stress intensity factor for K»=-0.2 MPavm, open circles: arrest condition Kiow=0; a) Kui 
vs. Aa for a; 36, b) Kia vs. Kappi for different initial crack lengths æ and Ko=0.2 MPaNm. 
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10 Anomalous subcritical crack growth 


10.1 Anomalous and normal subcritical crack growth 


Similar to other glasses, silica also shows the effect of subcritical crack growth. The 
fundamental effect, namely bond break processes, requires very high stresses. They are 
available in the singular stress field of cracks at crack-tip distances in the order of a few 
molecule diameters. Crack growth in silica glass has been well described by fracture me- 


chanics techniques representing the mechanical load by the stress intensity factor. 


Subcritical crack propagation of silica is usually presented in form of v( K)-curves where v 
is the crack-growth rate and K the externally applied stress intensity factor. Results are 
reported as a function of relative humidity and temperature in [10.1-10.4]. For silica, the 
astonishing effect of decreasing crack-growth rate v at an increased temperature was ob- 
served for constant partial water pressure in the humid environment. This surprising result 
observed in v- K experiments by Suratwala and Steele [10.5] is called anomalous subcriti- 


cal crack growth behavior. A comparison for different glasses is given in Fig. 10.1. 


Whereas for silica the crack velocity decreases with temperature, it increases for all other 
glasses [10.5]. This conclusion is of course somewhat problematic since curves (1) and (2) 
are not for silicate glass and curves (3-6) obtained in vacuum cannot describe any influ- 
ence that water might have. In addition, it is mentioned in [10.2] that silica in the vacuum 
tests failed not before toughness was reached. This implies a very steep v-K-curve and 
consequently subcritical crack growth at a value of K-0.8 Ki. could simply not be 


reached. 


Suratwala and Steele [10.5] assume that the negative temperature dependence will stem 
from a change in the slow crack growth resistance of the glass upon exposure to tempera- 
ture and stress. In the following considerations we will try to explain the effect by a sim- 
pler reason. In [10.6] we discussed the apparently incompatible temperature trends quali- 


tatively. 
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“Normal” glasses 


Fused silica (1), (2) Phosphate glasses 


(3)-(6) Silicate glasses 


Dn 
100 200 300 400 
Temperature (*C) 


Fig. 10.1 — Subcritical crack growth rate vs. temperature for different glasses at an individual stress intensity 
factor of K=0.8 Kic. Curves (1, 2, and silica): Measurements in N at 2700 Pa partial water vapour 
pressure. Curves (3-6): measurements in vacuum [10.2]. Curves (1, 2): Phosphate glasses [10.7, 
10.8]. Details can be found in [10.5]. 


10.2 Experimental crack growth results from literature 


In order to discuss anomalous crack growth, two basic results from literature are com- 
piled. In Fig. 10.2a v- K-curves for constant vapour pressure of p=2100-2200 Pa are plot- 
ted [10.5]. There is a clear decrease of crack velocity with increasing temperature visible 
and is the reason for the notation “anomalous subcritical crack growth". Measurements on 
silica in liquid water were published by Wiederhorn and Bolz [10.1] as shown in Fig. 
10.2b for the temperature region 2°C<7<90°C. From this figure it is evident that the 
crack-growth rates increase with increasing temperature. At first glance, this seems to be 


in contrast to the results of Suratwala and Steele [10.5]. 


In this Section it will be shown that the main difference between these contrary results is 
the fact of strongly different water vapour pressure in the two tests. Whereas in the tests of 
Fig. 10.2a the pressure was kept constant at 2.1-2.2 kPa, the saturation pressure varies in 
the water tests of Fig. 10.2b strongly, namely by a factor of 2100 (0.6 kPa at 2°C and 70 
kPa at 90?C). 


Consequently, the water concentrations at the crack tips are clearly different. This will be 


outlined in the following theoretical considerations. 
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Fig. 10.2 a) v-K-curves for constant vapour pressure of p=2.1-2.2 kPa by Suratwala and Steele [10.7], b) v- 
K-curves in liquid water by Wiederhorn and Bolz [10.1]. 


10.3 Expectation from reaction rate theory 


In [10.9] Wiederhorn et al. suggest for the subcritical crack growth rate vas a function of 


stress intensity factor K, and temperature 
-Q+bK 
V -ACk De] 29K) (10.3.1) 


where again Cis the concentration of molecular water [H20], Tthe temperature in °K, R 
the gas constant, Q the activation energy and Ai and b are proportionality constants. In 
eq.(10.3.1) the quantity &(7) is the rate coefficient for the single-step reaction showing 


the temperature dependency 


k.(T) «c T EE d (10.3.2) 


(A G-free energy of activation). From eqs.(10.3.1) and (10.3.2) we get with @=QtAG 
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Under a water vapour pressure p, and the solubility for molecular water at p-1 kPa, Cirea, 


the water concentration reads 
C=pPC (10.3.4) 


and equation (10.3.3) results in 


V =A, PCa T ep SP) (10.3.3a) 


When evaluating the temperature dependence, it should be noted that both pressure and 
solubility also vary with temperature. These influencing factors will be considered in the 


following Sections. 
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Fig. 10.3 — v-K-curves in liquid water by Wiederhorn and Bolz [10.1], transformed to a common saturation 


pressure of 3.2 kPa (saturation pressure at 25°C) according to eq.(10.3.3a). 
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10.4 Effect of vapour pressure 


The crack-growth rates from Fig. 10.2b are transformed by eq.(10.3.3a) to a constant 
vapour pressure of p=3.2 kPa (saturation pressure at 25°C), Fig. 10.3, exhibiting clearly 
the same temperature trend as found by Suratwala and Steele [10.5]. It should be noted 
that in the plot only the ordinate has been reduced to constant pressure and the abscissa 
values remain unaffected. From this plot one could claim that even in liquid water 


anomalous subcritical crack growth is present. 


10.5 Effect of water solubility 


Since the v- K-curves should depend on the water concentrations at glass surfaces, results 
of water concentrations are shown in Fig. 10.4 as a function of temperature. These data 
represent the concentrations of the molecular and hydroxyl water species derived in Sec- 
tion 4 of [10.10] on the basis of measurements by Zouine et al. [10.11] under saturation 
pressure. According to eqs. (4.1.3) and (4.1.4a) in [10.10] it results for the molecular 


water Cin mass units 


cu, (10.5.1) 
1+4k 
and for the hydroxyl concentration $ 
17 C, (10.5.2) 


S$-—.—t— 
18(1 1 
eur + s 
2 k 
(Cy = total water concentration, k= equilibrium constant). Figure 10.4 shows the concen- 


trations of the two water species Socp and Cocp normalized on the water vapour pressure 


of p=1 kPa. The straight lines can be expressed by 


zl exp|- 22.79 + nM (10.5.3) 
GNE ex} 154.4 29.09 7 zd (10.5.4) 
RT 
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Fig. 10.4 Solubility of water species at silica surfaces under saturation pressure derived from data by Zouine 


et al. [10.11], water concentrations normalized on the saturation pressure. 
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Fig. 10.5 — Crack-growth rates normalized on the molecular water concentration and temperature T, a) data 
from Fig. 10.2a, b) data from Fig. 10.2b. 
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The crack velocities of Fig. 10.2 are plotted in Fig. 10.5 in the normalized representation 
as V(Cx T) vs. applied K for the partial water vapour pressure of p=3.2 kPa (i.e. the satu- 
ration pressure at 25?C). It can be seen that the temperature tendency of Fig. 4b is now 
again as normal crack-growth behaviour. This makes clear that even for constant vapour 
pressure the anomalous crack-growth behaviour is at least partially caused by the reduced 


surface water concentration with increasing temperature. 


10.6 Effect of shielding stress intensity factor 


10.6.1 Constant shielding 


In our opinion, at least a part of the observed trend is caused by crack-tip shielding. In the 
preceding considerations we implicitly identified the externally applied stress intensity 
factor Kapı with the K-value acting at the crack tip, Kip. It was early outlined in [10.12, 
10.13] that a shielding stress intensity factor must develop as a consequence of volume 
increase by the hydroxyl generation (for references see e.g. [10.14]). 

Computations by McMeeking and Evans [10.15] result in a shielding stress intensity 
factor for a grown crack which is surrounded by a zone of thickness o undergoing a 


volume strain £v 


Sy ae Jo (10.6.1) 


with y=0.22 [10.15]. Eis Young's modulus and v Poisson’s ratio. In case of glass, œ is 
the thickness of the water diffusion zone around the crack tip. 

It has been shown in [10.16] that water-induced swelling at crack tips generates an intrin- 
Sic stress intensity factor Ka«O that shields a crack from an externally applied loading 
Kappi so that the stress intensity factor Kip at the crack tip is reduced. The effective stress 
intensity factor acting at the crack tip, Atip, represents the singular stress field at the tip. 
From the principle of superposition, Kip is given by 


KK GER, K <0 (10.6.2) 


tip app! 


1.e. Kip< Kappi. 
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The effect of the shielding stress intensity factor on the v-K-curves was outlined in 
[10.16]. As shown by eq.(1.1.9), the volume strain ev is related to the hydroxyl concentra- 
tion Sby the proportionality £oc.S. Consequently, it holds 


SE 
K,- vo (10.6.3) 


again with «=0.97. The hydroxyl concentration Sin the crack-tip stress field is thermally 
activated as eq.(1.1.3) implies. Since also the diffusivity is thermally activated, the water 
diffusion zone, w=f(VD Ð, increases with temperature due to the increase of the diffusivity. 
When in a subcritical crack growth test the shielding stress intensity factor equals — Kipyi, 
the total stress intensity factor Kip disappears followed by abruptly decreasing crack rates, 


i.e. the threshold values Ku of the v- K curve is reached. This is equivalent to 


(10.6.4) 


Threshold values for silica were evaluated by Sglavo et al. [10.17] on the basis of litera- 
ture results, [10.18,10.19]. The threshold stress intensity factor as a function of tempera- 
ture at «2 kPa partial water pressure was described by an Arrhenius temperature depend- 


ency 


Ka = Kao oa (10.6.5) 


Ka,970.638 MPa\m, g=1.83 kJ/mol. The value Kno in eq.(10.6.5) depends on the water 
partial pressure, for which unfortunately no measurements are currently available For this 
reason, the procedure described here is not suitable for real predictions, but may be suita- 
ble for quailitative interpretations. Equations (10.6.4) and (10.6.5) make it possible to 
calculate the shielding stress intensity factor Ks. and via eq.(10.6.2) the total stress inten- 
sity factor Kip acting at the crack tip. This value is responsible for the crack-tip stresses 


and subcritical crack growth rates. 


Finally, Fig. 10.6 gives the curves of Fig. 10.5 vs. the difference Kiypi- Ku. Under the as- 
sumption that Ksn#f( Kip), the curves in Fig. 10.6 agree with the v- Kip curves. The dia- 
grams of Fig. 10.6 show that, for the same water concentration, plotting against the stress 


intensity factor Kip (actually prevailing at the crack tip) leads to normal subcritical crack 


growth in both cases, namely increasing crack rates with increasing temperature. 
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Fig. 10.6 — Crack-growth rates from Figs.10.5a and 10.5b plotted vs. stress intensity factor Kippi- Ku. 


10.6.2 Variable shielding 


Deviations from the assumption of a constant shielding stress intensity factor Ki have to 
be expected. With increasing crack-growth rates the time for the water entrance via diffu- 
sion decreases and, consequently, the zone height o in eq.(10.6.3) becomes smaller. This 
yields in decreasing shielding stress intensity factor As with increasing v or increasing 
load. A computation of Ka(Kxpp, 7) was performed in [10.6] resulting in rather compli- 


cated equations and should not be duplicated here. 


An approximate method may be suggested to estimate Ki or, equivalently, Ki as a func- 
tion of Kapp. For this purpose, let us use the results of Section 8.2 for the case that shield- 
ing Ki disappears at fracture toughness, Fig. 8.4c. This situation is shown in Fig. 10.7. 
The total or tip-stress intensity factor Kir Kip holds according to eq.(8.2.9) 


Ku Kag (10.6.6) 


tip 
Ke ~ K, 


Figure 10.8 represents the v- Kip results for the data in Figs. 10.5a and 10.5b. Even from 


these plots normal subcritical crack growth can be stated. With increasing temperature 
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the crack growth rates for identical water concentrations at the glass surfaces increase, 


too. This holds for tests in liquid water as well as tests in humid air. 


Finally, Fig. 10.9a shows the stress intensity factor Kip necessary for reaching a constant 
normalized crack velocity v/CT. There seems to exist a discontinuity in Fig. 10.9a be- 
tween 75°C and 150°C. Figure 10.9b shows the applied stress intensity factor necessary 
to reach a crack velocity of v=10" m/s at a partial water vapour pressure of pu»o-2100- 
2200 as a function of temperature. In both representations, there seems a change of the 
crack-growth behaviour visible at about 140?C. A possible reason might be capillary 
condensation at crack tips. A zone of liquid water can be identified at a crack tip at low 
temperatures in humid air as had been shown by Ciccotti et al. [10.20], Grimaldi et al. 
[10.21], and Pallares et al [10.22]. At the high temperatures this is no longer possible. So 


one would expect a discontinuity between tests in “water” and tests in humid air. 


In(v/vo) 


Fig. 10.7 Suggestion for the estimation of the v- Kip curve. 


For the slopes of the v-K-curves one would expect from the theoretical analysis by 
Wiederhorn et al. [10.9] according to eq.(10.3.1) 


Bade (10.6.7) 


logv « BK, , RT 
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10.6 Effect of shielding stress intensity factor 


i.e. a decrease of B with increasing reciprocal absolute temperature 7: From the subcriti- 
cal crack growth rates of Fig. 10.8 the slopes, B, were determined. Additionally the pa- 
rameter B for a representation by log( WB Kappi is introduced for the data by Wiederhorn 
and Bolz [10.1]. The result is illustrated in Fig. 10.10a where the parameter B from 
eq.(10.6.7) is plotted. The squares show the results obtained from the data in [10.5] and 


the circles from data in [10.1]. 
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Fig. 10.8 — Crack-growth rates from Figs. 10.5a and 10.5b plotted vs. total stress intensity factor Kip. 


The values of B increase with decreasing temperature. This applies to all results with the 
exception of the data by Suratwala and Steele [10.5] for temperatures below 150?C, see 
also the jump in Fig. 10.9a. Below this temperature, the effect seems to be slightly re- 
versed. In Fig. 10.10b the parameter 5 defined in eq.(10.3.1) is shown. The parameter 5 
is for the representation via Kup a constant. This holds for the data from [10.1] and the 
data above ~100°C from [10.5] but decreases for lower temperatures as is included in 
Fig. 10.10b as the dashed line. The representation by eq.(10.6.7) yields (Standard Devia- 


tions in brackets) 


b=0.112[0.0017](m*’/mol) for [10.1] (10.6.8a) 
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b =0.158[0.0073](m*’/mol) for [10.5] and 2150*C (10.6.8b) 


For a representation by log( v)ocB Kapp, Wiederhorn and Bolz [10.1] suggested b=0.216 
m??/mol, introduced in Fig. 10.10b as the dash-dotted line. 
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Fig. 10.9 a) Stress intensity factors Kip for constant normalized subcritical crack growth rates v/(C 7) asa 
function of temperature, red: results in humid air by Suratwala and Steele [10.5], blue: results in 
liquid water by Wiederhorn and Bolz [10.1], b) applied stress intensity factors necessary for a 
crack-growth rate of v-107m/s at p=2100-2200 Pa pressure. 


10.6.3 Expectations for other glasses 


In the previous considerations, only silica was considered. Here, we will look on other 
glasses, too. Density measurements on several glass compositions were reported by 
Scholze [23]. The materials were: (B) B20:-glass, (K) K2O-SiO2-glass (20-80 mol-%), 
(Ca) Na?20-CaO-SiO»-glass (16-10-74 mol-%), (Al) Li20-Al20:-SiO»-glass (20-5-75 mol- 
96), and (Si) silica. 


The related volume strains £v of all the glasses according to Section 2 in [10.10] are plot- 


ted in Fig. 10.11. Whereas the B2Os:-glass shows strong volume shrinking due to water 


110 


10.6 Effect of shielding stress intensity factor 


uptake, silica shows clear volume expansion. The effects on the other glasses are rather 


small. Materials Ca and AI also swell. K shrinks slightly. 


In contrast to the crack-tip shielding for silica, it has to be expected that the B2O3-glass B 
should exhibit an increase of total stress intensity factors. Slight shielding effects may 


occur in the soda-lime glass as well as in the Al-glass. 
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Fig. 10.10 a) Slopes B of subcritical crack growth curves v( Kip) as a function of temperature, b) parameter b 
= BxRT, hatched area: scatter ranges of +1 Standard Deviation, dash-dotted line: suggestion by 
Wiederhorn and Bolz [10.1]. 


Conclusions: We studied subcritical crack growth results measured by Suratwala and 
Steele [10.5] in water vapour and by Wiederhorn and Bolz [10.1] measured in liquid 
water under saturation pressure. From the results shown before, the anomalous subcritical 
crack growth in the experiments by Suratwala and Steele [10.5] is a consequence of the 
decreasing surface water solubility with increasing temperature and the plotting against 
the applied stress intensity factor Kapp instead of the physically acting value Kip that 


governs the singular stress field at crack tips. With other words: Subcritical crack growth 
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10 Anomalous subcritical crack growth 


in silica is normal since for constant water concentration at crack tips the crack-growth 


rate at constant, Ki», increases with increasing temperature. 


It is also noteworthy that if the temperature dependence of the saturation pressure is elim- 
inated exclusively in the tests in liquid water [10.1], anomalous crack growth behaviour 


would appear (Fig. 10.3). 


This suggests that the representation of v-K-curves for constant pressure instead of con- 
stant water concentration in the glass and by plotting the crack-growth rates vs. the ap- 
plied stress intensity factor (instead of the physically relevant crack-tip stress intensity 


factor, Kip) is responsible for the anomalous behaviour of silica. 


Finally, it can be shown that the parameter bin eq.(10.3.1) is a constant if the crack-tip 


stress intensity factor Kip is used for K. 
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Fig. 10.11 Volume strains for several glasses from [10.10]. B: B2Os-glass, K: K20-SiO»-glass (20-80 mol-%), 
Ca: Na3O-CaO-SiO»-glass (16-10-74 mol-%), Al: LizO-Al203-SiO2-glass (20-5-75 mol-%), Si: Sili- 


ca 
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Consequences of hydroxyl 
generation by the 
silica/water reaction 


Part Ill: Damage and Young's 


Modulus 


11 Damage and Young's modulus 


11.1 Damage by hydroxyl generation 


When a hydroxyl has been formed, the initial silica ring is broken and the mechanical 
cohesion is weakened as is illustrated in Fig. 11.1. Such *defects" in the glass structure 
can be treated by using the damage variable D of continuum damage mechanics 
(Kachanov [11.1], Lemaitre [11.2]). This parameter is proportional to the density of 


micro-defects. 


According to the postulate of strain equivalence by Lemaitre [11.3], the effective elastic 


modulus, £p, decreases with increasing damage 
Ey - E,(1- D) (11.1.1) 


where £ is the modulus of virgin glass. The damage variable D can be determined from 
module measurements via eq.(11.1.1). The occurrence of this effect is in principle known 
from measurements on NaO-SiO; glass by Ito and Tomozawa [11.4]. To the authors’ 
knowledge, so far no comparable uni-axial measurements on fused silica are available. 
Therefore, we temporarily consider the damage D as a certain function of the hydroxyl 


concentration: 
D= f(S) (11.1.2) 


Apart from the equi-triaxial loading case with ox = ©y = o; (including the case of dis- 
appearing stresses), the elastic modulus must become a tensor with components 


depending on the degree of loading multiaxiality. 


Since this possibility would make the further treatment very difficult [11.3] and non- 
transparent, we assume in the following considerations that the damage remains isotropic 
and is considered to be of scalar nature. This is equivalent to the assumption of randomly 
orientated defects. Then also Æ remains isotropic. We assume that the defects in SiO», 


caused by hydroxyl generation, might behave like normal nano-pores. 
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Fig. 11.1 Volume element of silica showing damage by bond breaking due to the water/silica reaction, third 


dimension ignored. 


11.2 Experimental evidence for modulus reduction in 
silica 


In literature, there is experimental evidence for modulus decrease with increasing 
hydroxyl content. This can be seen from measurements of sound velocity as a function of 
water content. Measurements on longitudinal sound velocities in silica specimens with 
different water content were reported by Fraser [11.5] and Le Parc et al. [11.6]. 


Individual least-squares fits were made resulting for the data set by Fraser [11.5] 
v=5974(1-4.1855) in (m/s) (11.2.1) 
and the set by LeParc et al. [11.6] 
v =5959(1—5.34S) in (m/s) (11.2.2) 


When we normalize the results of the two test series on their individual values for S=0, 


we get the representation in Fig. 11.2a. A common straight-line fit of these data yields 


Y -1-BS (11.2.3) 
4 


0 


with the parameter 
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B=5.04 [4.23, 5.85] 


(90%-Confidence Intervals, CI, in brackets). The dependency of eq.(11.2.3) is introduced 
in Fig. 11.2a as the straight line. The longitudinal sound velocity depends on Young's 
modulus Fand density p via v oc JE/p : 

Whereas the densities for the measurements by Fraser are compiled in [11.5], we used for 
the results by LeParc et al [11.6] S according to the measurements by Shelby [11.7] (see 
also Section 2.1 in [11.8]): 


£u = pi pseu iv iA (11.2.4) 
E, 
with A=10.6 [8.7, 12.5]. It has to be noted that this value holds for isotropic damage since 
the natural OH-content doesn't show any preference for a special direction. The Young's 
modulus resulting from the data by Fraser [11.5] is shown in Fig. 11.2b. 
Figure 11.2b represents in addition stiffness coefficients Ci; by the triangles as were 
reported by Kushibiki et al. [11.9]. 
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Fig. 11.2 a) Longitudinal sound velocity in silica with different OH-content (blue circles: results by Fraser 
[11.5], red circles: results by LeParc et al. [11.6]), b) Young's modulus for the data of Fig. 11.2a 
together with stiffness results by Kushibiki et al. [11.9] as the triangles. 
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11.3 Modelling of damage by spherical pores 


Analytical computations on the reduction of Young's modulus with porosity were carried 
out by Wang [11.10] for spherical pores of equal size and constant spacing. 
Unfortunately, the results were given as tabulated data. They were the basis of a number 
of fitting relations, mostly using exponential functions, which trivially must fail for large 
porosities since they could not represent the requirement of Z=0 for a finite critical 


porosity. An overview can be obtained in [11.11, 11.12]. 


To overcome such problems, we fitted the analytical data by Wang [11.10] according to 


Éo = exp AP A p (11.3.1) 


with Pnax=0.4764 from [11.10] and the free fitting parameters A; and A; resulting as 
A, = 3.580 [3.565, 3.595] , A, = 3.747 [3.699, 3.795] 


with very narrow 90% confidence intervals indicating the good agreement with the 
numerical data. Figure 11.3a shows good agreement between eq.(11.3.1) and the 
numerical data entered by the circles. The dashed line indicates the initial slope that for 


small porosities reads 


En ~1_AP (11.3.2) 
E, 


and intersects the abscissa at a characteristic porosity P. of 
P, =1/ A, =0.278P, x (11.3.3) 


In order to give a better resolution for the agreement at low modules, Fig. 11.3b repre- 
sents the same results in logarithmic ordinate scaling. For our purpose we rewrite 
eq.(11.3.1) by normalizing the porosities on the maximum value Pmax and then replace P 
by $ 


2 
E B S 
—= exp -B (S/S ax) > 1/8 (11.3.4) 
E, -S/S max) Da 
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with parameters B\=1.667, B;-0.933. From eq.(11.2.4) and the initial slope of eq.(11.3.4) 


we obtain via By/Smax=A the maximum value of hydroxyl concentration 


max 


S = 0.157 [0.133, 0.192] (11.3.5) 


A simpler pore model was proposed by Phani and Niyogi [11.13] with the result of 


Fo LL apy (11.3.6) 


0 


This type of equation applied to the hydroxyl damage gives good agreement with the 


Wang solution over a large range of hydroxyl concentrations for the parameter set 
n=2, a=N225.3 (11.3.7) 
The modulus ratio can be described by 


Eo =(1-S/ Sau) (11.3.8) 


0 


The related dependency is shown in Fig. 11.3c by the red curve that gives disappearing of 
Young's modulus at Snax — 0.188 [0.16, 0.23]. 


baa aire ii pa a au Ea La 1 EP ee eee EN 
01 02 03 04 0 04 02 03 04 
Porosity P P max Porosity P Pmax 
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Fig. 11.3 a) Effect of porosity on Young's modulus for spherical pores; analytical results from Wang [11.10] 
(circles) compared with suggested fitting relation eq.(11.3.1) (curve); b) logarithmic ordinate 


scaling, c) Young's modulus as a function of hydroxyl concentration. 
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Fig. 11.4 a) Results of porosity on Young's modulus for silica by Adachi and Sakka [11.14], curve: best fit 


1 


according to eq.(11.3.6), b) comparison with results from longitudinal sound velocity in silica for 


different OH-content (red circles: results by LeParc et al. [11.6] and Fraser [11.5]). 
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11.4 Effect of damage in thin surface layers 


The fact that eq.(11.3.6) is an appropriate description for silica glass in the region of low 
moduli Ep£o—0, as can be seen from measurements by Adachi and Sakka [11.14], 
plotted by the circles in Fig. 11.4a. By fitting the data, we obtain the best representation 


for the parameter set: 
12, a=1.19 (11.3.9) 


introduced in Fig. 11.4a by the circles. Figure 11.4b shows all measurements from Figs. 
11.2 and 11.4a in a normalized representation. The dashed line illustrates a linear 


extension of the small-concentration data. 


11.4 Effect of damage in thin surface layers 


11.4.1 Global stress 


If the damage is concentrated on thin surface layers of a thickness d that is small to the 
undamaged bulk thickness W, d«« W, the strain in the damaged surface layer must equal 
the strain in the bulk: 


p- (11.4.1) 


The reduced average stress in the damaged surface, op, is therefore according to 
eq.(11.1.1) 


o,=0,1-D) (11.4.2) 


(o= stress in the undamaged bulk). 


External loading in z-direction will increase the distance between atoms of the ring 
structure, resulting in the strain £, in loading direction. Failure of the material is reached, 
when the strain reaches a critical value s.. The damage in the material affects also the 
load-bearing capacity of any area element since it reduces its net cross section and, 
consequently, the strength of the damaged material o;. If Ofo denotes the strength in the 


absence of damage, the strength of the damaged material is 


0,5 =O, - D) (114.3) 
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The strength of water-damaged silica fibres is strain-controlled for d«« W. This means 
that the strain at failure £. is the same in the surface layer and the bulk. This condition 


yields 


29,55 O9;,0-D) Fro (11.4.4) 
"Ey. E-D) E, 


Equation (11.4.4) tells us that the damage problem can be handled simply by using bulk 
material parameters (Æo, 0:0) since the damage by broken bonds cancels out in (11.4.4). 


For some more details on fiber strength, see Section 11.7. 


11.4.2 Effect on crack resistance and stress intensity factors 


The damage of the initial ring structure of silica, Fig. 11.1, must also affect the resistance 
against crack propagation. The crack growth resistance G, represents the energy ne- 
cessary to split all bonds that are broken when the crack has passed the considered 
volume element. Application of the damage variable yields for the crack growth resis- 


tance of the damaged material, Gp, 
G.» 2G. (1- D) (11.4.5) 


(G.o- crack resistance for the undamaged material, D-0) or in terms of stress intensity 


factors 


Ks = /E Gap -K,(-D) (11.4.6) 


where K.» is the reduced toughness of the damaged material and Ki. the fracture 
toughness measured in fracture mechanics tests on the undamaged material. 

In the special case of a crack fully embedded in a material of reduced but constant 
modulus, Æp, the problem simplifies strongly. The applied stress intensity factor is 
generally given by 

=0,,,F,/7a (11.4.7) 


appl appl 


K 


(F= Fracture mechanics geometric function, a-crack length, Fig. 11.5). 
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Using the actual stress from eq.(11.4.2), Kappi reads for the crack fully embedded in the 


damaged zone 


K spt > = Cap 0 Fra L- D) = Ky) o(1—D) (11.4.8) 
——————À 


Kappi,o 


where Kappio Stands for the applied stress intensity factor formally computed with stresses 


Oo as present in the bulk material 


Roe name (11.4.9) 


appl,0 


From (11.4.6), (11.4.8) and (11.4.9) it results equivalently to eq.(11.4.4) 


K oppı,o a T D) = K sooo (1 1.4.10) 
c,D K,(1- D) K; 


K appi, D AA 


K 


Also in terms of stress intensity factors, surface cracks fully embedded in the surface 
layer can be handled simply by using bulk material parameters (Kapp, Ki.) instead of the 


unknown parameters in the layer. 


11.4.3 Effect on pre-existing dimples 


By removing the plastic protective coating on fibers with an acid, defects will be 
generated that can be described as hemispherical pits [11.15]. Depending on their 


geometry, such defects can reduce the surface strength clearly. 


At the deepest point of a pit surrounded by the hydroxyl-damaged material under tensile 
stress Op, the stress is increased by the stress-concentration factor ß, so that the maximum 


stress at the deepest point (A), Fig. 11.5, is 
Gun =P Op (11.4.1 1a) 


For a half-sphere it is for instance B=2. Before the hydroxyl-damaged layer was gene- 


rated, it did hold equivalently for the undamaged material under stress oo 


Oo = P9 (114.11b) 
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11 Damage and Young's modulus 


From (11.4.11a) and (11.4.11b) we have to conclude that in the case of failure the 
strengths in presence and absence of hydroxyl damage are reduced both by the common 
factor 1/B. When the strengths in presence of the pit are denoted by orpip and Ofpito, it 
holds 


1 
O f, pit,o Tp t (11.4.12a) 
=e (11.4.12b) 
OF, pit, D ~ Bot D 


Apart from the fact that the strengths are reduced, it can be seen from the combined 
eqs.(11.4.12a) and (11.4.12b) 


Of pit, _ Of, pio (11.4.13) 
Orp 9,0 


with the consequence that also in this case the hydroxyl damage has no influence on 


strength. 


Crack 


Hemispherical 
pit 


1 
1 


Fig. 11.5 Cross-section through a crack and a partially embedded hemispherical pit in a hydroxyl-damaged 
surface layer of thickness &<W (We radius of the fiber) and a<<d. 
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11.4 Effect of damage in thin surface layers 


11.4.4 Effect of swelling stresses 


If the volume swelling occurs in a constrained manner, such as at the glass surface, then 
the swelling causes a build-up of compressive stresses in the surface region, which can 


strengthen the glass significantly as has been outlined in [11.16]. 


The compressive swelling stresses in axial and circumferential directions, y, z, 
respectively, are generally proportional to the product of volume swelling strain £, and 


Young's modulus £ 


0, 0, © EE (11.4.14) 


sw,y? SW,Z 


For the most simple case of isotropic swelling, the factor of proportionality is for instance 
1/(3(1-v)), where v is Poisson’s ratio. 


The related volume swelling strain, £, 70, can be expressed according to eq.(1.1.9) by 


€,=KxS (11.4.15) 


with the coefficient k=0.97. 


When for example 6,,;0 denotes the axial swelling stress in the undamaged material, it 
holds as a consequence of eq.(11.1.1) that the swelling stress in the hydroxyl-damaged 


surface is 


Osw,z,D = w04- D) (11.4.16) 


As an example, this dependency may be applied to the crack-like surface defects. 
The so-called shielding stress intensity factor Ka«0 that shields the crack tip partially 
from the applied load is in the damaged material 


Kap 70, ,5F4£a 260, , F4za(1- D) - K,,(1—-D) (11.4.17) 


sw,z,D 


where Kio stands for the shielding stress intensity factor in the undamaged material. 
In presence of crack-tip shielding, failure occurs, when the total stress intensity factor Kip 
reaches fracture toughness. In the damaged material the principle of superposition results 


in 


«0 (11.4.18) 
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11 Damage and Young's modulus 


Using eqs.(11.4.6), (11.4.8) and (11.4.17) gives the same result in terms of the stress 


intensity factors in the undamaged material 


Kup -(K *KS3(-D) -K,,ü- D) (11.4.19) 


appl,0 


The equal sign on the right of (11.4.19) makes clear, that the shielding effect can be de- 
scribed by stress intensity factors in the undamaged material. 

For the case of the surface additionally damaged by pit-like surface defects, the tensile 
stress Op, has to be replaced by the total stress op*o;,;p resulting, equivalently to 
eq.(11.4.19), in 


O mnax, D = B(op + OD) = P(o + ou, A D) = Oso- D) (11.4.20) 


Again from the equal sign on the right, it becomes obvious that treatment with damaged 
quantities gives the same results as are obtained when hydroxyl damage is completely 


ignored. 


11.5 Example of application: Uniaxial stress-strain 
curve 


11.5.1 Swelling excluded 
As has been outlined in [11.17], the hydroxyl concentration is stress-dependent. In order 


to allow a transparent derivation, swelling stresses may be neglected first. The hydroxyl 


concentration as a function of stress can be written 
AV. 
S = S, exp o —41- (11.5.1) 
RT 


At temperatures >500°C, the parameter A Vef reads for uniaxial loading [11.17] 
AV „= 14.4 cm?/mol (11.5.2) 


The total strain ¢ under the applied stress is 
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11.5 Example of application: Uniaxial stress-strain curve 


_ E api (11.5.3) 
E,(S) 


where for &p(.S) the equations (11.3.4) or (11.3.8) have to be introduced. The maximum 


stress, asymptotically reached for se, is 


" RT Log S, (11.5.4) 
SEU HAV S 


max 


In Fig. 11.6 the o-g-curves are shown for 650°C (355 Torr water vapour pressure) and for 
500°C with parameters So (6 Torr pressure) taken from [11.17]. The dashed lines indicate 


the two asymptotes for £0 and &> of the solid curves. 


2000 
Oo 
1500 
(MPa) 
1000 
0=650°C 
500 So=0.65wt% 
8=500°C 


So=0.16wt% 


4 1 1 L 1 L 


0.1 0.2 0.3 


Fig. 11.6 — Stress-strain curves for 650°C, 355 Torr and 500°C 6 Torr, hydroxyl concentrations So from 
[11.17]; Approximate description of the stress-strain curves by eq.(11.5.5) (dashed lines). 


For simpler mechanics computations, the stress-strain curves can be approximated by the 
two asymptotes (dashed lines in Fig. 11.6) where the rising part represents Young's 
modulus and the horizontal line the saturation stress given by eq.(11.5.4), i.e. a uniaxial 


"yield stress" oy =Omax. Then it holds 
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11 Damage and Young's modulus 


[07 


Ee for e< = 
" E, (11.5.5) 
o, 
Oy for g2— 
E 


0 


2000 r 9-500'C 2000 
So=0.16wt% 
O | oO 
1500 f 1500 
(MPa). (MPa) 
0=650°C 
1000 So=0.65wt% 1000 
500 500 
Without swelling 
Swelling included 
1 L i- L 1 L L L L 
0.1 0.2 0.3 0.1 0.2 0.3 
E€ AE=E—Ewz,0 


Fig. 11.7 Stress-strain curves, curve for 650°C and 500°C, black: without swelling, red: swelling included, 


a) total strain as the abscissa, b) strain by external loading as abscissa. 


11.5.2 Swelling included 


When swelling strains are taken into consideration, eq.(11.5.3) has to be replaced by 


pac ee (11.5.6) 


where £,,,; is the swelling strain in axial direction. According to [11.17], the eq.(11.5.1) 


yields 


Mer] (115.7) 


Ew ZAKS - a KS, explo appi PT 
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11.6 Deviating behaviour under compressive loading 


with «=0.97, a=1.92 [11.17]. The stress-strain curve obtained via eqs. (11.5.6) and 
(11.5.7) is shown in Fig. 11.7a as the red curves together with the curves by eq.(11.5.3) as 
the black ones. The offset of the red curves comes from the fact that even in the absence 
of stresses a swelling strain occurs that is caused by the initial hydroxyl concentration Sp. 
When E£sw.z0 is the swelling strain caused by the initial concentration So, the strain increase 


under load is A£-£—-£,,;o. This is shown in Fig. 11.7b. 


11.6 Deviating behaviour under compressive loading 


In the previous considerations, we applied the model of pore-like defects. This model de- 
scribes symmetrical material response under tension and compression loading. For spe- 
cial cases of applications, the description by crack-like defects may be of advantage. This 
model allows introducing non-symmetry of deformation in tension and compression. Due 
to the fact that the crack surfaces cannot interpenetrate each other under compressive nor- 
mal stress, the reduction of Young's modulus is less strong than in tension, represented 


by 


(11.6.1) 


E,A-D) fo 020 
E,A-hD) for o«0 


with the crack closure parameter Ox x1, often found to be A=0.2 [11.18]. 


11.6.1 Shielding-stress intensity factor for damaged material 


A consequences of damage by hydroxyl generation is the reduction of Young's modulus 
E. When Æp is the modulus in the damaged state and Æ the value for undamaged silica, 


we could derive the relation [11.19, 11.20] (see Section 11.3) 


= =(1-44S) =(1-S/S,,,)° (11.6.2) 


0 


with A=10.6 [8.7, 12.5] (90%-CI in brackets). The hydroxyl concentration at which the 
Young's modulus disappears is Smax =2/A= 0.188 [0.16, 0.23]. 
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11 Damage and Young's modulus 


Including the varying modulus in the computation of the shielding stress intensity factor 
requires a modification of the equations in Section 9.2.2. Now the stress intensity factor 


results from 


K,-—— [(s ZEE (11.6.3) 
S(O)E, } 0o 


The concentration profiles by Lechenault et al. [11.21], Fig. 9.4, clearly show the charac- 
teristic shape with a plateau as would be expected from the Irwin and Dugdale crack-tip 


zone models for "plastic flow" behaviour. 


11.6.2 Shielding stress intensity factor for asymmetric damage 


In the case of disappearing total stress intensity factor at the crack tip, Kiota = Kapp + Kup <0, 
all volume elements in the water-affected zone are under compression due to volume 
swelling by hydroxyl generation [11.22]. In the previous considerations, we applied the 
model of pore-like defects. This model describes symmetrical material response under 
tension and compression loading. The Young's modulus in compression is not necessarily 


identical in tension and compression. 


Figure 11.8b again shows the shielding stress intensity factor for the hydroxyl concentra- 
tion according to eq.(11.6.3) for different values of the non-symmetry coefficient or crack 
closure parameter 0</A<1. The asymmetry gives reason to define an upper limit for the 
shielding stress intensity factor. The zone in front of the crack tip is under tension for a 
load of Ki, 70, Fig. 11.8a. The crack flanks, on the other hand, are stress-free with regard 
to the applied load and are under pressure due to the swelling stresses. In this region the 
lower expression in (11.6.1) holds. As the limit approximation, the reduced modulus in 
the zone ahead of the tip is assumed to be free of any stress (£50) and the crack-face 
region is assumed to exhibit the initial modulus, i.e. #0. This approximation is de- 
scribed by using the value y — 0.37 according to McMeeking and Evans [11.23]. Finally, 
Fig. 11.8b represents results from [11.19] for an erfc-shaped hydroxyl distribution with a 


surface value S(O) and a parameter b of 4.95 nm (solid curves), according to 


S - S(0) LEA (11.6.4) 
2b 
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11.7 Global strength of a fiber 


The dashed line in Fig. 11.8b shows the upper limit solution for the same S-distribution 


addressed before. 
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Fig. 11.8 a) Negative and positive contributions to the total shielding stress intensity factor, b) shielding 


stress intensity factors for an erfc-shaped hydroxyl distribution, eq.(11.6.4), as a function of the 
surface concentration .S(0). Solid curves: from [11.21], dashed curve: via eq.(9.2.2) with y=0.37, 
EB. 


11.7 Global strength of a fiber 


The hydroxyl-damage in the glass also affects the load-bearing capacity of any area 


element, since it reduces its net cross-section, Ap 
Ap = A1- D) (11.7.1) 


where Ao denotes the total geometrical cross section subsuming damaged and undamaged 


regions. 


When again Oto denotes the strength in the absence of hydroxyl damage, and orp the 


strength in the hydroxyl-damaged state, the strength results as given in eq.(11.4.3). 
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11 Damage and Young's modulus 


Let us now consider the strength of silica fibers. First, the water distribution may be 
replaced by C=Co=const. over the thickness b and C=0 in the bulk indicated in the inset 


of Fig. 11.9. The total cross section area A of the fiber with radius Ris 


A- Rz (11.7.2a) 
that of the inner circle of radius R-5 
Am=(R-b)’r (11.7.2b) 
and the ring-shaped area is 
A,=R’r - (R- bz - (2R- byrb (11.7.2c) 


Failure in the undamaged fiber occurs at the surface, when the stress reaches the strength 


Oro. The force A at failure is then 
P, -g, RU (11.7.3) 


When the surface is damaged by an amount of D, the strength in the layer 5 decreases to 
Orp according to eq.(11.4.3). Since the Young's modulus decreases in the same way, 
eq.(11.1.1), it is clear that the critical strain €. remains unaffected as is visible from 
eq.(11.4.4). 


Fig. 11.9 Failure load of a fiber for differently damaged surface regions as a function of layer thickness b. 
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11.7 Global strength of a fiber 


The Bernoulli Theorem tells us that plane cross sections must remain plane. It implies 


that the strain in the damaged and undamaged parts of the fiber must be the same. 


This again implies that at failure o=o¢o is in the inner area Ainm and o-orp-or(1-D) in 
the surface layer. The total load Py of the damaged fiber at failure is the sum of the load 


in the inner (undamaged) area and the load in the damaged ring-shaped zone, i.e. 


P, 2c (R-b) z +0, ,(1- D(2R-b)bz (11.7.4) 


The reduction of the fracture load follows from eqs.(11.7.3) and (11.7.4) 


P, _ (R by t (1 D)b(2R b) E 2bD b’D (11.7.5) 
P R? R R 


0 


and is plotted in Fig. 11.9 for several damage values D as a function of the layer thick- 
ness. It becomes clear, that for small ratios of 5/ R«0.02, the strength (computed from the 
load Pp) is hardly affected. 


In the special case of thin layers, b<<R, it results approximately 


Py wp (11.7.6) 
B, 


Figuratively speaking, it can be said that even in the absence of any load-carrying 
contribution from the stresses in the surface layer, D-1, the strength of the entire fiber is 
hardly affected. 
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12 Shielding affected by damage 


12.1 Effect of local swelling at surface cracks 


The reaction of water and silica in the surface diffusion zone, eq. (1.1.1), affects the frac- 
ture mechanics stress intensity factor K at the tips of cracks. As a consequence of the 
high near-tip stresses, very high hydroxyl concentrations must occur in the crack-tip 
region. This allows the assumption that near a crack tip the water is present predominant- 


ly in the form of hydroxyl water. 


Due to restrictions in free expansion, the swelling strains £y result in swelling stresses 
which give rise for the "shielding" stress intensity factor Ka<0. At first loading, the 


swelling zone at the crack tip 7(on,@) is heart-shaped as illustrated in Fig. 12.1a. 
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Fig. 12.1 a) Swelling zone at the tip of an arrested crack under mechanical loading, caused by stress- 


enhanced diffusion, b) zone for a crack grown by Aa. 


The shielding stress intensity factor caused by a swelling zone of height o is [12.1, 12.2] 


E 
K,, =-y Vo (12.1.1) 


where Fis Young’s modulus and v Poisson’s ratio. 
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12 Shielding affected by damage 


12.2 Shielding stress intensity factors for a damaged 
material 
A consequence of damage by hydroxyl generation is the reduction of Young's modulus 


E. When = is the modulus in the damaged state and Æ the value for undamaged silica, 
we could derive the relation (11.3.8), [12.2] 


E» (1.8/5, Y (122.1) 
E, 
with Smax =0.188 [0.16, 0.23]. 
n n 
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Fig. 12.2 a) Water profile on DCDC specimen on the faces of a growing crack fractured in liquid water by 


Tomozawa et al. [12.3]. 


Including the varying modulus into the computation of the shielding stress intensity fac- 


tor requires a modification of eq.(9.2.6). Now the stress intensity factor results from 


Can fisxe( 2) doi (122.2) 
S(0)E, } 8o ), 
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12.3 Shielding stress intensity factor for asymmetric damage 


In order to explain the evaluation of eq.(12.2.2), again the measurements by Tomozawa 


et al. [12.3], plotted already in Fig. 9.5, may be used here again, Fig. 12.2. 


12.3 Shielding stress intensity factor for asymmetric 
damage 


In the case of disappearing total stress intensity factor at the crack tip, Ki = Kappi +Ktip 
X0, all volume elements in the water-affected zone are under compression due to volume 
swelling by hydroxyl generation [12.4]. In the previous considerations, we applied the 
model of pore-like defects. This model describes symmetrical material response under 
tension and compression loading. The Young's modulus in compression is not necessarily 


identical in tension and compression, see Section 11.6. 


Figure 12.3b again shows the shielding stress intensity factor for the hydroxyl concentra- 
tion according to eqs.(12.2.1) and (11.6.1) for different values of the non-symmetry coef- 


ficient or crack closure parameter Oc x1. 


5 10 15 20 
Z (nm) 


Fig. 12.3 Distribution of the product Sx E relevant for eq.(11.2.2), computation for A=1 in eq.(11.6.1), b) 
for the representations by Fig. 12.2a (blue curve) and Fig. 12.2b (red curve). 
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12 Shielding affected by damage 


The product SxZ that governs the shielding stress intensity factor according to 
eq.(12.2.2) is shown in Fig. 12.3. The blue curve shows the result for the representation 
by Fig.12.2a and the red curve by Fig. 12.2b. The region at a depth of z~5-6 nm has the 


strongest influence on shielding. 
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Fig. 12.4 Shielding stress intensity factors for an erfc-shaped hydroxyl distribution, as a function of the 
surface concentration S(O) from Fig. 11.8b. Square: Result for S(0)= 16.5% and 420.05. 


Figure 12.4 again represents the results of Fig. 11.8b for an erfc-shaped hydroxyl distri- 


bution with a surface value S(O) and a parameter b of 4.95 nm (solid curves). 


An effective value of A can be estimated from the experimental result of S(0)=16.5%, a 
crack-growth rate of v z10?m/s, and a related applied stress intensity factor Kappı=0.43 
MPavm (see Fig. 8.4a). According to eq.(8.2.8), we obtain for Ki-=0.75 MPavm and 
Ka=0.3 MPa\m [12.5]: Ka=-0.213 MPa\m as introduced as the square in Fig. 12.4. This 


value is equivalent to an asymetry parameter of 


h= 0.05 (12.3.1) 
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12.4 Mechanical interpretation of the limit case under compression 


12.4 Mechanical interpretation of the limit case under 
compression 


The question remains how to explain such an asymmetry in the modulus of elasticity. A 
possible interpretation of the asymmetry coefficient A is illustrated in Figure 12.5a, where 


a water molecule may attack a load-carrying bond of the silica net structure. 


A SiO»-ring of the silica network may be described by a spherical (nano)-void in a con- 
tinuous material. This allows a fracture mechanics treatment by an FE-study. According 
to the finite crack-tip curvature for cracks, the radius of such a void has to be expected in 
the order of p 2 0.5 nm [12.6], Fig. 12.5a. 


Under a remote tensile load in z-direction, G», the red Si-O-bond may be cracked by 
reaction with molecular water (blue/red). In the tensile case the opened bond cannot 
transfer any significant tensile forces, Fig. 12.5b. Consequently, it would be A>1 in 
eq.(11.6.1). 
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Fig. 12.5 a) Bond breaking under uni-axial tension by reaction with molecular water (attacked Si-O-bond 


in red). 


Due to the initially load-carrying but now fractured bond, the stress concentrations in- 


crease at the ends of the now tablet-shaped void as indicated by the red circles. Due to 
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12 Shielding affected by damage 


the exponential increase of the hydroxyl concentration with stress, newly reacting water 


molecules will predominantly react at these locations. 
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Fig. 12.6 a) Additionally cracked bond at a location of high stress concentration, b) limit case of a wedge- 


like contact of the hydroxyls under compression. 
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APPENDIX 


Al Modifications of shielding by addi- 
tional effects 


A1.1 Deformed steady-state swelling zones 


The computations of the swelling zones are mostly carried out for diffusion in a motion- 
less coordinate system with a fixed crack tip, i.e. for an arrested crack. Whereas for a 
growing crack the transversal diffusion normal to the crack plane is hardly affected by 
the moving crack tip, the diffusion and the crack propagation compete in crack direction 
[Al.1, A1.2]. The consequence for strongly stress enhanced diffusion zones with cordial 
shape is a reduced ratio of the zone size in length direction m to the zone height o (Fig. 
A1.1a), which for an arrested crack is [A1.3] 


n(0 8 (A1.1.1) 


o(0) 427 


In this relation æ(0), @(0) stands for m(v=0) and & v-0), respectively. 


arrested crack crack growing under 
steady state conditions 


Fig. Al.1 a) Swelling zone ahead the tip of an arrested crack, b) zone with reduced length m(v) for a crack 


growing at a constant crack-growth rate v (schematic). 
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Al Modifications of shielding by additional effects 


For the stress intensity factor of the deformed zones one can use some single results from 
literature. Depending on the subcritical crack growth rate, a growing crack can show a 
crack-tip diffusion zone with a reduced value of m/ as studied in [A1.4]. Such a zone is 
illustrated in Fig. A1.1b. The stress intensity factors for three different ratios of r/o are 
available from [A1.3]. These cases are illustrated in Fig. A1.2a. Figure. A1.2b shows the 


related coefficients y for the general representation 


E 
K,- No (A1.1.2) 


by the circles. The interpolating curve through the data points in Fig. Al.2b. Figure 
Al.2can be expressed simply by 


yz o baz "o ) (A1.1.3) 


[2] 


This interpolation relation enables the computation of stress intensity factors for de- 


formed diffusion zones. 


8/4(27) 


lol (A) 


Fig. A1.2 a) Swelling zones with different shape at the zone end, b) shielding stress intensity factors by 


eq.(A1.1.2) for the zones given in a), expressed by the coefficient y. 
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A1.2 Crack leaving a swelling zone 


A1.2 Crack leaving a swelling zone 


At higher temperatures and after long time, also unloaded cracks may develop a swelling 
zone of thickness 5 by water diffusion into the crack faces, Fig. A1.3a. In this case, the 
zone end is modeled by a semicircle. If the swelling strain is constant in the swelling 


Zone, £v-go-const., and &,=0 outside the zone, the shielding stress intensity factor is 


REDEN? (A1.2.1) 
i-v 

with the coefficient y depending on the zone shape in the crack-tip region. If such a 
crack grows by an amount of Aa, it will escape from the initial swelling zone, as illustrat- 
ed in Fig. A1.3a-Fig. A1.3d for different zone ends. The shielding stress intensity factors 
for two limit cases of the zone shape near the crack tip are plotted in Fig. A1.3e. For the 
circular zone end (Fig. Al.3b), the shielding stress intensity factor is given by the 
squares. A straight zone end with m/b=0 yields the dash-dotted curve. The two dependen- 
cies can be described by eq.(A1.2.1) where now y is a function of Aa/b. In order to 
minimize the uncertainty in the true shape of the zone ends it was recommended in 


[A1.5] to use for computations »/b = 0.5 with coefficient y —0.3. 


The stress intensity factor solution for the case zy/b-'^ (Fig. A1.3d) can be obtained by 
linear interpolation of the two limit cases resulting in the dashed curve of Fig. A1.3e. The 


coefficient y may be approximated by 


Bn <i 
yas 0.305 + i; Aa/b for ^a/lb&4i (A122) 
0.34exp[-24Aa/b - 4 ] Aalb>+ 


It should be mentioned that the curves in Fig. Al.3e change abruptly caused by the as- 
sumption of a sharp layer boundary. For an arbitrarily varying strain distribution, the 
weight function or Green’s function method from [A1.6] can be used. The stress intensi- 
ty factor caused by a swelling layer of thickness dd’ in the surface distance Z=b’ (a) un- 


dergoing the strain &(b’) results as 
d(K,,) = K,,(b'+db') - K,,(b') (A1.2.3) 


This gives rise for a Green’s- or “weight function” / defined as 
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= d(K,, ) 
db' 


(A1.2.4) 


For any swelling distribution &(z)=e(b’), the related shielding stress intensity factor is 
then obtained by summing over all the zone increments db multiplied with the local 


swelling strain 


K,- — feo) h(b')ab" £, - &(b'- 0) (A1.2.5) 


-— —. Aalb 


Fig. Al.3 a) Crack with a thick circular swelling layer of constant swelling strain, b) crack after additional 


—n 


extension Aa for a circular zone, c) for a rectangular zone, d) for a zone with n/b=1/2, e) varia- 
tion of the shielding stress intensity factor for a crack growing in and outside the initial layer ac- 
cording to [A1.3] (dash-dotted curve) and [A1.5] (squares); dashed curve: interpolated from the 


continuous and dash-dotted curves for m/b=1/2. 
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A1.2 Crack leaving a swelling zone 


In the absence of stress-enhanced diffusion, the swelling profiles are represented by the 


complementary error function. The volumetric strain at the location z=’ is in this case 


E = &, erfc 2 (A1.2.6) 
D,t 


0 


where & is the strain at the free surface, Do the diffusion coefficient and ¢ the time after 


the first water contact. 


Equation (A1.2.5) has been evaluated for the case of n/b=\. The resulting shielding 
stress intensity factors related to the erfc-shaped swelling profiles are shown in Fig. 


A1.4b. The numerical result for m/b=2 can approximated by equation (A1.2.1) with 


0.305 (A1.2.7) 


2 
3 
b 


Fig. Al.4 a) Green's function as the shielding stress intensity factor of a thin swelling layer of thickness b’, 
b) stress intensity factors according to the dashed curve in (d) for an erfc-shaped strain distribu- 


tion (insert); curves according to Fig. A1.3e. 
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A1.3 Blunt crack with finite notch root radius 


Two crack-tip models are commonly used in fracture mechanics for the description of 
stresses ahead of a crack. Whereas continuum mechanics assumes the existence of a 
mathematically sharp crack, in a micro-structurally motivated approach the slender 
notch with a finite radius at the crack end is considered. Stresses at a sharp crack tip are 
described by the crack-tip stress intensity factor Kip. They become singular if the tip is 
approached, o—>œ for r0. However, the K-description via continuum mechanics must 
fail for very short distances, when the tip distance competes with the microstructure of 
the material. In the micro-structurally motivated approach, the crack tip region can be 
considered as a slender notch with root radius p in the order of the average radius of the 
SiO» rings. Wiederhorn et al. [A1.7] suggest a crack-tip radius p — 0.5 nm for the silica 


network. 


Whereas shielding stress intensity factors are well known for sharp cracks [A1.6], there 
is a lack in multiaxiality behaviour of swelling stresses for the slender notch model. The 
individual swelling stress components and also the hydrostatic stress terms are necessary 
for computations on stress-enhanced equilibrium of the silica/water reaction. In contrast 
to internal swelling, the stress components by externally applied loading are known from 
Creager and Paris [A1.8]. 


A1.3.1 Finite element model 


Finite Element computations in order to determine swelling stresses for the slender notch 
model as a function of the distance from the notch root and differently high swelling 
Zones were carried out in [A1.9]. A single crack embedded in a swelling zone of height 
h and constant volume swelling is shown in Fig. A1.5. The crack length, a, is small 
compared to the specimen dimensions, realizing the edge-cracked semi-infinite body. On 
the other hand, the crack length is at least a factor of 10 larger than the notch radius p so 
that a "slender" notch is sufficiently guaranteed. Zones of different heights were 
modelled in a finite element (FE) study. Solid continuum elements (8-node bi-quadratic) 
were chosen and the computations were carried out with ABAQUS Version 6.8. The 
volume strain was replaced by the equivalent thermal problem by heating the swelling 


zone by A7 keeping 7=0 in the rest of the structure. This results in the strain 
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A1.3 Blunt crack with finite notch root radius 


£, 23e AT (A1.3.1) 


(a=thermal expansion coefficient) with isotropic expansion in x-, y-, and z-directions 


E, =E, +8, +E, (A1.3.2) 


Requirements for the z-coordinate were chosen as “plane strain" as is fulfilled at crack 


tips. 


Fig. Al.5 Slender notch at a surface with swelling zone of thickness A. 


A1.3.2 Results 


Swelling stresses in x-, y-, and z-directions are plotted in Figs. A1.6a-d for zone heights 


of A/p=0.2, 0.5 and 1.0. The black data represent the hydrostatic stress, defined by 


0,=1(0,+0,+0,) (A1.3.3) 


Whereas all stresses are compressive within the swelling zone, they roughly disappear 
outside for the z-direction and for the hydrostatic stress. The x-component remains in 


compression and the y-component becomes tension as a consequence of the equilibrium 


condition. 
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Fig. Al.6  a)-c) Swelling stresses, normalized on the volume expansion strain for 3 different zone heights 


h, d) Effect of zone height A on swelling stresses at the notch root, x=0. 
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A1.3 Blunt crack with finite notch root radius 


The surface stresses (x0) are represented in Fig. Al.6d as a function of the relative 
height A/p of the swelling zone. The squares indicate theoretically known values for an 
infinitely thin surface layer: At a free surface, the stress state is plane stress and, conse- 


quently, also stresses caused by swelling are equi-biaxial (0-0) 


gigs m. (A1.3.4) 
$ 30 -v) 


where Fis Young's modulus and v is Poisson's ratio. For silica it is v=0.17, consequent- 


ly: o/e«E = —0.4016 for the y- and zdirections. The hydrostatic stress reads 


o, 2i(o,*0,)- i m (A1.3.5) 


i.e. ow(£vE)-0.2677. 


In order to characterize the “multiaxiality” of the swelling stresses, we plot in Fig. Al.7a 
the hydrostatic stress normalized on the zcomponent. This parameter is 2/3~0.67 for 
thin swelling layers and decreases slightly for on/oy and increases for on/oz. The varia- 


tion over the swelling zones is negligible as can be concluded from Fig. A1.7b. 


xip 
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Fig. Al.7 a) Effect of zone height 4 on multiaxiality at the notch root, x=0, b) variation of the multiaxiali- 


ty ratio on/oz over the swelling zones. 
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A1.3.3 Effect of finite notch length 


From Fig. Al.6a-d it is visible that the hydrostatic swelling stress shows a small 
variation over the swelling zone. For an understanding of this small effect, we increased 
for b/p=1 the notch length to a/p=30. Figure A1.8a shows the effect for on. In this plot 
the data for a/p=10 are represented by the black symbols, those for a/p=30 by the red 
ones. It is obvious that the variation over the region 0<x<b clearly decreases. In Fig. 
A1.8b the effect on the normal swelling stress oy is shown. Here, the variation remains 


unaffected over the zone size. Only the absolute stress values o/e,F shift slightly. 
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Fig. A1.8 Effect of the notch length on the hydrostatic and the normal swelling stress. 


A1.3.4 Including shielding stress intensity factors 


Stresses at slender notches were given by Creager and Paris [A1.8]. The stress compo- 


nent normal to the crack plane (the stress oy) is in distance x from the notch root 


2K — p*x (A1.3.6) 


aom 
” az(p-2x)p-*2x 
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A1.3 Blunt crack with finite notch root radius 


The other stress components are 


PEE: x (A1.3.7) 
"omo 2x) p*2x 
_ 2K (A1.3.8) 


AZ (p + 2x) 


Consequently, the hydrostatic stress term from the applied stresses results as 


2(l+v)K (A1.3.9) 


o = 
"^^ Sapa 2x) 


Equations (A1.3.6-A1.3.9) are valid for applied stress intensity factors Kapp, shielding 


stress intensity factors Ka, and total stress intensity factor Kip, which are defined as 


Kip = Kao + Ky (A1.3.10) 


tip — “Sapp 
The shielding stress intensity factor for an edge crack of finite length a, embedded in a 
swelling zone of height b, was already computed via FE in [A1.10]. The results can be 


expressed by the relation 


K,,(b/a) = Ka Oep 1.182 (A1.3.11) 
a 
with 
K,,(0) = -y (0) as vb (A1.3.12) 
=V 


and different zone shape [A1.3] 


0.22 for heart shape (A1.3.13) 


0.25 forsemi circle 


y(0)= | 


Now we can use this solution for the computation of the shielding stress intensity factor 
in eqs.(A1.3.6-A1.3.9), resulting in swelling stresses contributed by Ks. The resulting 
stresses are then subtracted from the results in Fig. A1.6d. The difference is the contribu- 


tion by local swelling exclusively, i.e. for Ki;-0. 
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Al Modifications of shielding by additional effects 


Two different procedures are now possible for the computation of swelling stresses: 


e Fitting the dependencies of Fig. A1.9 and adding the swelling stresses caused 
by the shielding stress intensity factor. 


e — Use of the data in Fig. Al.6d directly because these values include already the 
shielding effect. 


Since the resolved data of Fig. A1.9 vary clearly stronger with the distance x from the 
notch root it is recommended to use results of Fig. Al.6d. For the hydrostatic stress 


term, it is suggested for 0« xx 5 


DEE 
a: ; 2)- er [s ora? | (A1.3.14) 
p pj 9ü-v) p 


Fig. A1.9 Swelling stresses for a disappearing shielding stress intensity factor, Kn=0. 
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A2 Finite Element study on heart- 
shaped zones 


A2.1 Zone shape 


The singular hydrostatic near-tip stresses are given as 


o,=:(+v) cos(g!2) (A2.1.1) 


4/2007 


where r and @ are the polar coordinates with the origin at the crack tip, see Fig. A2.1a. 
Equations (9.1.1) and (A2.1.1) imply that in the high crack-tip stress field nearly all water 


is present in form of hydroxyl S. 


The shape 7(@) of the zone contour for constant o; in plane strain results from eq.(A2.1.1) 


as 


pow pens (B9) (A2.1.2) 


33 


with the height œ of the zone for a prescribed hydrostatic stress 


2 
— (+v) [ K, (A2.1.3) 
pev Teles 
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Fig. A2.1 a) Contour for a constant hydrostatic stress ahead of a loaded crack, b) variation of Young's 


modulus in a frontal zone A(x) (dashed curve), replaced by an average (effective) modulus £. 
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A2 Finite Element study on heart-shaped zones 


The contour of constant hydrostatic stress according to eqs.(A2.1.2, A2.1.3) is shown in 
Fig. A2.1a as the heart-shaped curve. This contour also describes the damage zone in 
front of the tip of a crack that has not yet grown by subcritical or stable crack growth. 
The hydroxyl generation by the reaction (1.1.1) causes damage since the originally intact 
silica ring structure is cracked by the water attack. One of the consequences of such 
damage is the reduction of Young's modulus F (see Section 11). 

As Fig. A2.1a shows, the tip of the crack is not embedded in the heart-shaped zone but 
terminates at =O. This point in principle belongs to both the frontal zone and the un- 
damaged surrounding material. A description by eq.(11.3.8) leads to a variable module 
distribution F(x), which is replaced here by an average constant value for reasons of 
simplicity, Fig. A2.1b. The module Æ prevails outside the crack-tip zone. Another reason 
for the decrease in the elasticity module in the "frontal process zone" is microcracking in 


polycrystalline ceramics as was studied by Evans and Faber [A2.1]. 


A2.2 Crack terminating at ọ=r of the heart-shaped 
zone 


A2.2.1 Constant modulus over the zone size 


First we studied the case that a constant modulus is prescribed over the whole zone. In 
the heart-shaped zone the Young's modulus is F< and outside the zone Æ. FE 
calculations were performed with ABAQUS 6.9. In order to obtain normalized stress 
intensity factors K according to (A2.2.2a), @=1, Oappi=l, and Æ=1 were chosen. Table 
A2.1 shows the J-integrals according to Rice [A2.2] and the stress intensity factors 
obtained from the first contours for the lowest Young's modulus of £E-0.01 Æ. Apart 
from the first contour, the results are almost constant. Due to symmetry, no mode-II 
stress intensity factor could appear. The average mode-I stress intensity factor over 
contour 2-8 is 1.974 and the Standard Deviation 0.00125 (0.06%). The interrelation 
between crack opening displacement v and stress intensity factor Ki is for the near-tip 


field along the symmetry line (here: -x=r<<o) in plain strain with £'—E/(1-v?) 
Kee (eve (A2.2.1) 
8 Jr 
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A2.2 Crack terminating at p=n of the heart-shaped zone 


Normalized displacements v represented in Fig. A2.2 are defined as 


= E, (A2.2.1a) 
Op OL-V*) 


The stress intensity factors obtained from the least-squares fit of the displacements in Fig. 
A2.2 are compiled in Table A2.2, Column 3. In this context, it has to be noted that crack 
opening displacements COD are defined as the distance from the symmetry line to one of 
the surfaces. Unfortunately, the (total) crack opening, i.e. the distance between the crack 


flanks, is sometimes also referred to as COD. 


Contour J Kı Kn 
1 370.3 1.953 0 
2 378.0 1.973 0 
3 378.0 1.973 0 
4 378.1 1.973 0 
5 378.6 1.975 0 
6 379.2 1.976 0 
7 379.6 1.977 0 
8 378.8 1.975 0 


Table A2.1 J-integral and stress intensity factors for the first contours (Z=0.01, Æ=1, v=0.17, &=1) in 


normalized form. 


An independent procedure for the stress intensity factor determination is the evaluation of 
the near-tip stresses ahead of the crack tip, i.e. the stresses along the prospective crack 


propagation plane. For the normal component on this plane it holds for the first terms 


PIE vee US (A2.2.2) 


?w J2zXx 


A regression analysis of the data in Fig. A2.3 provides the stress intensity factors. The 
results obtained from the data fit of Fig. A2.3 are compiled in Table A2.2, Column 4 in 


the dimensionless form of 
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A2 Finite Element study on heart-shaped zones 


ce (A2.2.2a) 


Fig. A2.2  FE-results of near-tip CODs fitted according to eq.(A2.2.1). 
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Fig. A2.3 Normalized stress Oyy /Oappı along the prospective crack propagation line (@=0), lines: fitting curves 
according to eq.(A2.2.2), yielding the stress intensity factor for the data points in the range of 
0.0005<x/m<0.004. 
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A2.2 Crack terminating at ọ=r of the heart-shaped zone 


A2.2.2 Step-shaped modulus distribution 


So far only cases with a constant modulus were studied, Fig. A2.4a. An additional FE- 
evaluation was made for the case of a variable modulus in the zone ahead the tip. As Fig. 
A2.4b shows, the modulus was varied in 4 steps. In the inner zone of Fig. A2.4b, the 
modulus Z was £E-0.25 E, then E-0.57 E», E-0.72 E», and outside Z=%. The result from 
the J-Integral was (for Æ=1, v=0.17, Oap=l1 and w=1): K’=14.83, hardly different from 
the case of constant modulus £-0.25 A. Also this value is introduced in the second 
Column of Table A2.2. 


K, eq.(A2.2.2a) 

E! Eo via J via COD via Oyy 

1 29.21 29.11 28.78 
0.5 21.82 21.81 21.58 
0.25 15.43 15.40 15.21 

0.25/0.57/0.72/1. 14.83 - a 

0.1 9.039 8.871 8.892 
0.03 4.141 4.153 4.065 
0.01 1.974 1.997 1.943 


Table A2.2 Normalized stress intensity factors K from J-Integral, COD, and near-tip stress oyy in the form of 
eq.(A2.2.2a). 


The K-values for the single zones of Fig. A2.4a, compiled in Table A2.2, are plotted in 
Fig. A2.5a by the symbols (for the geometric data see the insert). In this plot the stress 
intensity factors were normalized on the J-based X value at E/Ej-1. The red triangle 
corresponds to the varying Young's modulus shown in Fig. A2.4b. Figure A2.5b repre- 
sents the ratio of the data for COD- and o,,-based stress intensity factors normalized on 
the stress intensity factors from the J-Integral. The hatched area indicates the region of 


data scatter with an average value and "scatter" of about 


K(B)/Ky(£)=0.993+0.015 
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Fig. A2.4 a) Crack-tip zone with a Young's modulus Æ deviating from the bulk modulus Fo, E<Fo, b) 


stepwise decreasing modulus in zones of height 1/3, 2/3, and 3/30. 
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Fig. A2.5 a) Stress intensity factors from J-Integral, COD, and near-tip stress Oyy, (K-values normalized on 
the result from J-Integral for Æ Fo-1), red triangle for step-shaped £-distribution, b) comparison of 
K from COD and oyy with the individual results from J-Integral 


A2.2.3 Conclusions from the FE-analysis 


From the results of Fig. A2.5 it becomes evident that 


e the stress intensity factors are roughly independent of the computation method, 


e the stress intensity factor is only affected by the modulus closest to the crack tip, 
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A2.2 Crack terminating at p=n of the heart-shaped zone 


e the straight-line behaviour of Fig. A2.5a suggests the representation by 


E 
KERG I (A2.2.3) 
0 


The FE-results confirm eq.(A2.2.3) as had to be expected from the theoretical analysis by 


Merkle [A2.3] on slender notches embedded in the reduced modulus at the notch-root. 


From the agreement between the theoretical solution by Merkle [A2.3] and the FE-results 
for 0.1<E/A<1, it can be concluded that our crack-zone modelling and the accuracy of 
our FE-mesh is sufficient at least for this range. For the lowest value of E/£5x0.01, the 
FE-program indicated convergence problems. With regard to the data accuracy it should 
be noted that the theoretical analysis was carried out for an infinite body (w/a, w/H, 


œ/ W)—oo, while the FE calculation was carried out on a finite one. 


For conditions comparable with “small-scale yielding" in elastic-plastic fracture mecha- 
nics, i.e. o << (a, W-a), the effect of the zone with reduced E on Kapı becomes a “null- 
effect" [A2.2]. Then it holds 


Ka = oyF(a/W,H IW) fra (A2.2.4) 


appl 


with the geometric function F(a/ W, HI W) and the remote normal tractions o» at the plate 


ends. 


We used in our study a ratio of /a=1/30 and expected that the zone would be sufficiently 
small to fulfill geometry conditions for small-scale yielding. Despite this relatively small 
zone, influences from this zone can still be seen in Fig. A2.5a. This may at least partially 
be caused by the influence of an effective zone length a..=atA@, A20 with A=0 for the 
case of homogeneous material. 

In the case of A>0, eq.(A2.2.4) reads 


Kopi = FoF (Qs /W,H IW) raus (A2.2.5) 


The stress intensity factors for E/E5-0.25 and w/a=0, 1/30, and 1/90 are plotted in Fig. 
A2.6. A clear trend with @/a is visible, showing an increase of the stress intensity factor 


with increasing relative zone sizes. 


169 


A2 Finite Element study on heart-shaped zones 


0.54r 
K(E) 


Ko) | 


0.537 


Fig. A2.6 Stress intensity factor from J-Integral for (Z/E=0.25). 


A2.3 Derivation of J-Integral for slender notches by 
Merkle 


Merkle [A2.3] showed that the relation (A2.2.3) is exact in the case of a slender notch 
ending in a region of Young's modulus Æ different from the bulk material Zo. In his 
analytical derivation he assumes a half-infinite crack in an infinite body loaded by remote 
tension in y-direction resulting in the applied stress intensity factor Kapp. Our FE- 
computations were carried out on a finite rectangular plate (inset in Fig. A2.5a). 

The stress intensity factors were evaluated via the J-integral defined by 


J 2 (way - T2 as (A2.3.1) 
: Ox 
id 


where Wis the stain energy density, I the path around the crack tip, T the tractions on 
the path, uthe displacement in y-direction, and ds a length increment of T. 

For reasons of simplicity we used the usual path along the surfaces represented in Fig. 
A2.7 with the line segments (1)-(8). Along the contours (3), (5) and (7) no tractions act 
and on (4) and (6) constant tensile stresses are prescribed. 


The J-integral along the surface segment (1) is 
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A2.3 Derivation of J-Integral for slender notches by Merkle 


E fw » (A2.3.2) 


n 


In order to simplify the path-integration, the slender notch was approximated by a narrow 


elliptic notch with parabolic end shape (Fig. A2.8). 


Fig. A2.7 Integration path for J-evaluation (red arrows) 


An arbitrary point on the curved surface of Fig. A2.8c is described by the geometrical 


interrelations 


OD NONE L (A2.3.3) 


y-rsing, ly 
1+ coso 1+cos@ 


First it has to be noted that in the J-Integral only the modulus at the surface, A, enters and 
the variation with the distance from the surface doesn’t play any role. In a subcritically 
grown crack the reduced value £i is constant along the whole notch surface. 


Merkle [A2.3] showed that the elastic strain energy along the crack surfaces (1, 2, 8) is 


W= o; _K* cos’ (p/2) (A2.3.4) 
2E. E zr 


S S 
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A2 Finite Element study on heart-shaped zones 


with the tangential stresses o. From the geometrical relations in eq.(A2.3.3) it follows 


after a simple mathematical manipulation 


w- K G+ cose)’ (A2.3.5) 
E. 2np 


Consequently, the J-integral along path parts (1, 2, 8) reads 


K? 


J= 
E 


1 


17 K? 
ode cos g)d o > J, Ds (A2.3.6) 


S 


This of course holds for large zones, @>p to œ>>p, for which @— 7 can be sufficiently 
fulfilled. Since the results of eq.(A2.3.6) and (A2.3.7) no longer depend on the notch 
radius, they should also apply to the sharp crack by taking p—0. The remote (or 
"applied") stress Oapp at lines (4) and (6) cause the applied stress intensity factor Kapp) and 


the related “applied J-integral” 


js Koi (A2.3.7) 


a) 


Fig. A2.8 a) Contour at the tip of a subcritically grown crack, b) detail of a), c) crack surface approximated 


by an elliptical notch showing the same curvature radius p, schematic. 


172 


A2.4 Results for extended zones 


The J-integral along a closed path must disappear as shown by Rice [A2.2], i.e. 


K? 2 
J =J] om eae |Ë (A2.3.8) 


It should be mentioned, that the derivation by Merkle [A2.3] deals with the case of a zone 


of different Young's modulus ahead of a crack tip with finite notch root radius p>0. 


When the radius p becomes smaller in the limit-case consideration, the geometry of Figs. 


A2.8b and A2.8c remains self-similar and not necessarily represents the sharp crack. 


A2.4 Results for extended zones 


While in Section A2.2 only zones in front of the crack tip were considered, in this Section, 
we additionally compile results for zones extending over the entire crack length. In this 
case, the crack tip is fully embedded in the zone of reduced module. The plate 
dimensions and the crack length were increased by a factor of about 3 in order to reduce 


any influence of w/a furthermore. This is shown in the plots of Fig. A2.9. 


1 W=200 o 
D 


K(E) | K from J 
A O Fig. A2.5a 
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1 
Kj/Kg. 


0.5 


0.2 


0.1 


0.05 


0.02 


0.01 0.02 0.00 01 02 0.5 1 


Fig. A2.9 a) Stress intensity factors from J-Integral, Ks (normalized on the value Ko obtained for F=) for 
embedded crack tips (squares) compared with the K;-results of Fig. A2.5a (plotted by circles and 
triangle), solid curve: Approximation via eq.(A2.4.1), b) Stress intensity factors from J-Integral, 


dashed line: asymptote with slope=1/2, dash-dotted line: asymptote with slope=1. 


The stress intensity factors exclusively on the basis of the J-Integral are given in Fig. 
A2.9a. The squares for Æ Æ <0.25 were fitted by 


K(E) _ 2E (A2.4.1) 
K(E,) E+E, 


This expression is introduced by the solid curve. 
We also evaluated COD-results as done in Section A2.2 and obtained the near-tip profiles 


given in Fig. A2.10. Stress intensity factors Kcop were computed from the Irwin-parabola 


v- lÈ Koon Jr (A2.4.2) 
a El(l-v^) 


where vis the opening displacement and v Poisson’s ratio (0.17 for silica). The displace- 
ments in the crack wake are given for the two geometries @/a=30 and 100. In both cases 


al We al H-0.5 was kept constant. The plot v vs. Vr showed very straight lines. 
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A2.4 Results for extended zones 
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Fig. A2.10 Crack-opening profiles close to the tip versus (7/0)? for a crack of /@=100. 


The stress intensity factors from eq.(A2.2.1) are plotted in Fig. A2.11a as the solid circles 
together with the results from the J-Integral in Fig. A2.9a as the squares. In this represen- 
tation no differences of the two K-values can be seen. Figure A2.11b shows the ratio of 
the two stress intensity factors. Note the significantly higher resolution of the ordinate. 
The mean value and the standard deviation SD (in brackets) are 


Ton — 0.999987 [0.001134] 


J 


Although in principle no theoretical facts can be proven by numerical calculations, this 


good agreement of the X factors speaks for the validity of K; using different methods. 


Conclusions: From the FE results and the theoretical solution by Merkle [A2.3] it was 


obtained 


appl 
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A2 Finite Element study on heart-shaped zones 


When comparing the results for the extended areas of reduced module (here: length L~a) 
with the solution according to Merkle [A2.3], no agreement can be expected a priori. The 
zones on which the FE calculations are based do no longer meet the basic requirements of 


the analytical calculation, namely (œ, L)<<(a, W,A). 


1 F Oe fully embedded crack 1.01r 
K/Ko | © heart-shaped tip-zone Kcop L 
| L +1SD 
0.5 K; t 
0.2 F 
0.99 F 
Oir f b) 
0.05f aab 
0.02 B 
4 L ji dB 1 ll L ipu eal. 0.97 L Lr Ppp ı Pp a | I 1 
0.01 0.02 0.05 01 02 0.5 1 0.2 0.4 0.6 0.8 1 
1/2 
E/E, (E/E) 


Fig. A2.11 a) Ki vs. E/ Fo as squares, results from Fig. A2.5 as open circles, Kcop from Fig. A2.10 via 
eq.(A2.2.1) as solid circles, b) ratio Kcop/ Kr. 


A2.5 Crack-Opening Displacements (COD) 


We evaluated COD-results in Section A2.2 in order to compute the stress intensity factors 
for different Young's moduli and obtained the near-tip profiles. Stress intensity factors 
Kcop were computed from the Irwin-parabola, eq.(A2.4.2). Again v is the opening 
displacement in the distance r from the crack tip and v=0.17 is Poisson's ratio. 


Normalized displacements v represented in Fig. A2.12 are defined by 


E 
ie 0 (A2.5.1) 
“KE, od v") 
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A2.5 Crack-Opening Displacements (COD) 


where A(£) denotes the stress intensity factor for the case of homogeneous material. In 


the following COD in larger distance are considered. 


A2.5.1 COD for heart-shaped zones 


Crack-opening displacements for the crack with a heart-shaped zone ahead of the tip are 
shown in Fig. A2.12a in the form v’=f(7/o) and in Fig. A2.12b in the form of (v’)’=f(7/o). 
From the diagram of (v’)?=f(7/@) it becomes obvious that the squares of CODs tend to 
straight lines for about 71.5. The straight lines for E/Ej«1 show fictive crack origins 
as indicated by the intersections of the dashed extensions of the crack profiles with the 
abscissa. In the case of very short crack-tip distance the COD is similar to crack-tip 
blunting in elastic-plastic fracture mechanics although the material behaviour is purely 
linear-elastic. In larger distances of about 7/0»1.5, the straight line for E/E5-0.01 shows 
nearly the same slope as obtained for E/£Ej-1 and seems to be shifted by the constant 
amount (v (0)} (see the inset in Fig. A2.13). 
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Fig. A2.12 Influence of reduced Young's modulus in the heart-shaped crack-tip zone, a) linear representation 


of the normalized COD according to eq.(A2.5.1), b) square of the COD as the ordinate. 
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A2 Finite Element study on heart-shaped zones 


The values of the straight lines v at 7-0 are called in elastic-ideal plastic fracture 
mechanics the “Crack-Tip-Opening Displacements (CTOD)" and used as material 
property for yielding. Figure A2.13 again shows the COD-profiles for the extremely 
different Young's moduli ratios £/E;-1 and 1/100. The apparent crack-tip blunting is 
illustrated by the inset in Fig. A2.13, when the COD is simplified by the thick red 
contours, resulting in an apparent crack-tip opening and an apparent increase Aa of the 


fictive crack length. 


EE, d( v ?/d(/o) (V0) Ado 
1 2.436 0 0 
0.5 2.433 1.396 0.573 
0.25 2.494 2.746 1.101 
0.1 2.488 3.966 1.594 
0.01 2.511 5.197 2.067 


Table A2.3 Parameters of straight lines in Fig. A2.12b: Slopes of straight lines, opening displacements 
(V (7=0))? at the crack tip, fictive crack-length increase Aa as a function of Young's modulus ratio 
El Eo for v=0.17. 
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Fig. A2.13 Representation of near-tip COD by the crack-tip opening displacement. 
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A2.5 Crack-Opening Displacements (COD) 


A2.5.2, COD for zones extended in the crack wake 


For the zone extending over the whole crack wake, we performed the same evaluations as 


done for the heart-shaped zone. Figure A2.14 and Table A2.4 show the results. 


El Eo d(v)"/d(ro) (v)-0) Aala 
1 2.58 0 0 
0.5 2.69 0.47 0.17 
0.25 2.76 1.32 0.48 
0.1 2.81 2.25 0.80 
0.05 2.83 2.74 0.97 
0.01 2.86 3.19 1.12 


Table A2.4 Parameters characterizing the dashed straight lines in Fig. A2.14. 
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Fig. A2.14 Influence of reduced Young's modulus in the extended zone, square of the COD as the ordinate. 
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A2 Finite Element study on heart-shaped zones 


A2.6 Comparison of parameters for the two zone 
shapes 


Figure A2.15 shows again the parameters from Tables A2.3 and A2.4. In Fig. A2.15a the 
slopes are given and in Fig. A2.15b the offsets of the straight line asymptotes at 7/«—0. It 
is obvious that the slopes hardly vary noticeably with the modulus of elasticity. Only very 


slight changes with increasing modulus can be stated. 


Summary on COD: We studied the COD-behaviour of cracks in a bi-material. The 


following points may be emphasized: 


e From our FE-computations on crack-opening displacements it can be stated that the 


effect of a reduced Young's modulus is visible only in a crack-tip distance of /o<1.5. 


e The slopes of the asymptotes, reached for distances /@>1.5, are within about +5% 


independent of the modulus. 


e The fictive crack lengths defined by the intersections of the straight asymptotes with 


the abscissa is in the order of Aa/m « 0-2. 


e The steeply rising displacements at 7/0 may suggest to the observer an occurrence 


of "crack-tip blunting", although purely linear-elastic material behaviour is present. 
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Fig. A2.15 a) Comparison of slopes and b) crack-tip displacements from Tables A2.3 and A2.4. 
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A3 Hydrostatic stress at crack tips 


A3.1 Hydrostatic stress under hydroxyl damaging 


The hydroxyl generation by the reaction eq.(1.1.1) causes damage since the originally 
intact silica ring structure is cracked by the water attack. One of the consequences of such 
damage is the reduction of Young's modulus Æ. In order to describe this E-decrease, we 


used in [A3.1] the rather simple damage model proposed by Phany and Niyogi [A3.2]. 


When Æp is the modulus in the damaged state and Æ the value for undamaged silica, we 


could derive the relation (see Section 11 and [A3.1]) 


Po = (1-8 / Spa) (A3.1.1) 


0 


with Smax=0.188. Under the condition that the damage is isotropic, the multiaxiality in 
damaged states remains unchanged. In a volume element that is hydrostatically loaded by 


tensile stresses 6,=oy =0,, the externally applied hydrostatic stress is 
Oy api = $(0,* 0, + 0) (A3.1.2) 


The hydroxyl concentration as a function of stress, reads by using eq.(A3.1.1) 


AV AV 
S=S, LE A S, © Onam (1-S/S,.,.)° | (A3.1.3) 


In Fig. A3.1a, the hydroxyl concentration is plotted vs. the applied hydrostatic stress for 
So=0.1 wt-%, and 0-90?C. The dash-dotted line indicates the maximum hydroxyl con- 
centration of Smax=18.8 wt-%. Figure A3.1b shows the related hydrostatic stress in the 
damaged state, Orp, vs. the applied hydrostatic stress HOnapp. In the notation of the 


coordinate axis the abbreviation p is defined as 


uet. (A3.1.4) 
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A3 Hydrostatic stress at crack tip 


The maximum hydrostatic stress, asymptotically reached for On.app >%, is 


HO}, p max x LE | (A3.1.5) 


max 


The curves 6 np(Onappi) of Fig. A3.1b may be simplified by a bi-linear description as intro- 
duced by the dash-dotted lines in Fig. A3.1b. 


In this case, the maximum possible stress in the damaged state is depending on the 


temperature. In the sense of an upper limit solution, the stresses can be approximated as 


On appl for On appl < OD, h,max (A3.1.6) 


[07 = 
h,D 
o h,D,max else 
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Fig. A3.1 a) Hydroxyl concentration vs. applied stress under condition of constant strains, b) hydrostatic 
stress in the damaged material on,» as a function of the applied hydrostatic stress On.appı and the 
initial hydroxyl concentration So with the abbreviation u defined by eq.(A3.1.4), limits of 


hydrostatic stress and approximation by an upper bound description (dash-dotted straight lines). 
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A3.2 Irwin and Dugdale model 


A3.2 Irwin and Dugdale model 


Due to the damage via hydroxyl generation, the stress-strain behavior at a crack tip is no 
longer linear. By application of the bi-linearized dependency eq.(A3.1.6), fracture 
mechanics problems can approximately be solved applying for instance the well-known 
Irwin and Dugdale models, originally developed for elastic/plastic material behaviour. 
The Irwin [A3.3] and the Dugdale model [A3.4] describe p/astic deformations in a zone 
ahead of the crack. Whereas in the Irwin model the size of the "plastic zone" was com- 
puted from “yield”-condition Oappiy=oy and by force equilibrium along the prospective 
crack plane, the size of the Dugdale zone results from the condition that the stress 
intensity factor for the fictive crack in Fig. A3.2a disappears. So the shape of the zones is 
not fixed, especially not for the Dugdale zone. The latter is mostly concentrated along the 
X-axis. 

According to these models the “length” L of the “plastic” zone, is under plane stress 


conditions 
KY Ko 
T 1 T =| Kan! (A3.2.1) 
A 


where oy is the “yield stress" for the case of a linear-elastic/ideal-plastic material. The 
opening at the tip of the physical crack, called Crack Tip Opening Displacement 
(CTOD), à, results from [A3.3, A3.4] as 


2 2 
"uu Es (A3.2.2) 
ü rwin cT Eo, t ug: ale Eo, 


The CTOD, 6, is illustrated in Figs. A3.2b and A3.2c. Whereas the Irwin model [A3.3] is 
preferred for circular "plastic zones", the Dugdale model [A3.4] describes also plastic 
deformations in a narrow strip in front of the crack. In the derivation of the Dugdale 
formula, it was assumed that there exists an internal through-the-thickness crack in an 
infinite plate and that the yield stress is constant over the zone size L regardless of the 
locally varying strains. This provided the analytical solutions (A3.2.1) and (A3.2.2) for 
the zone size and the CTOD. 


Remark: It should be emphasized here once again that it is not claimed that the glass 


behaves in an elastic/ideally plastic manner. This approach only leads to a simple fracture 
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A3 Hydrostatic stress at crack tip 


mechanical treatment of the crack tip problem for the "deformation" curves represented 


by the limit case eq.(A3.1.6). 


a) 


fictive crack tip 
with K=0, (DUGDALE) 


notational crack tip, (IRWIN) 


Fig. A3.2 a) Stress distribution curve near a crack tip (black curve) approximated by an elastic-“plastic” 
material behaviour (red curve) in the crack-tip region (zone shape is not necessarily a circle), b) 


crack tip opening displacement (CTOD) &, c) observable crack profile. 
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A4 Swelling stresses from crack- 
terminating angles 


A4.1 Estimation of the magnitude of residual stresses 


A4.1.1 Shielding stress intensity factors and residual stresses 


Volume expansions occur by swelling of the glass volume due to hydroxyl generation 
[A4.1, A4.2, A4.3]. The results from [A4.1] can be expressed by eq.(1.1.9) in terms of 
the hydroxyl concentration Sin mass units. 

Due to the mechanical boundary conditions, the swelling strains result in swelling stress- 
es proportional to the amount of hydroxyl concentration as were described by 
eq.(1.1.11). A volume element in a thick plate that undergoes swelling cannot freely 
expand. If the diffusion zone is small compared to the component dimensions, expansion 
is completely prevented in the plane of the surface and can only take place normal to the 


surface plane. 


In Section 2 we determined the shielding stress intensity factor from the externally ap- 
plied stress intensity factor Kapp: and the total stress intensity factor Ki. In the case of 
silica, Ki, is for spontaneous and stable crack extension identical with the toughness Ki 
and was in our rather fast subcritical crack growth tests on silica about Kip=0.55-0.6 
MPam. Therefore, we used an average of Kiip=0.575 MPa\m. 

The crack-shielding stress intensity factor must be proportional to the swelling stresses 


Osw acting normal on the crack plane and the square-root of the layer thickness 
K,, = F(o,b)o, zb (A4.1.1) 


The fracture mechanics geometric function for this loading over a small region, A, b), is 
so far unknown, but may be estimated by not too large effort. Very often the stresses 
caused by diffusion processes (here: water diffusion into silica) and the stress distribution 


can be described by 


187 


A4 Swelling stresses from crack terminating angles 


O w = O w.o erfc ea (A4.1.2) 
: 2b 


where Osw,o is the surface value, z=0, of the stresses and bis the diffusion depth at which 
the water concentration and consequently the swelling stresses are reduced to =50 % of 


the surface value as is illustrated in Fig. A4.1a. 
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b z 
Fig. A4.1 a) Stresses in the surface region according to eq.(A4.1.2), b) semi-circular surface crack loaded 


by the erfc-shaped stresses. 


A4.1.2 Stress intensity factor solutions for residual stresses 


In order to give an estimation for the maximum swelling stresses in eq.(A4.1.2) on the 
basis of the experimentally obtained shielding stress intensity factor, the geometric func- 
tion Họ,b) must be known. Following an approximate procedure proposed by Under- 


wood [A4.4], the interrelation between shielding and applied stress intensity factors reads 


— K,,(90°,b) (A4.1.3) 
K pp (90 ) 


K, (p, b) - 4 
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A4.1 Estimation of the magnitude of residual stresses 


However, this procedure can yield larger errors if the shapes of the two cracks to be com- 
pared are significantly different. In the original paper by Underwood [A4.4] it was shown 
that the stress intensity factors for complicated cracks (semi-elliptical surface cracks in 
thick-walled tubes) under complicated stresses (tension and bending superimposed) were 
always smaller than those for the reference crack (continuous crack in a plate under pure 
tension), compare e.g. Fig. 6 in [A4.4]. This fact may be indicated by the parameter Ax1 
in the following eq.(A4.1.3). The equal sign, A=1, trivially holds when the crack to be 
predicted is identical with the reference case. From the special result of [A4.4], we can 


conclude for the coefficient A, introduced in eq.(A4.1.3), 1«2/3. 


The shielding stress intensity factor is available for the special case of a crack terminating 
with @=90°. For this purpose let us replace the straight crack by a semi-circular surface 
crack of depth a that might be large compared to the layer thickness, a>>b. The stress 
intensity factor at the surface point (B), Fig. A4.1b, obtained from [A4.3], is 


Kap 21290, , Ja aman 2 + oa (A4.1.4) 
í i a a 


A series expansion of the FE-results in [A4.3] for b/a—0 results in the simple expression 


Ke = 1.720 bec umb (A4.1.5) 


When calculating K-factors, it does not matter whether the applied stresses are externally 
applied or internal ones. Therefore, both must show the same type of angle dependencies. 
In our first derivation in [A4.5] it was implicitly assumed that the coefficient A can be 
regarded as a constant value. From the comparison with experimental results in Fig. 4 of 
[A4.5], however, a dependence on the terminating angle emerged, A=A(@). This finding 
should be taken into account in the current evaluation. The shielding stress intensity 


factor Kx(q, 5) can now be obtained as 


1 
Ka (g,b) = A(p O woy Tb (A4.1.6) 
N) 0.00658 9 +7.66x10 gf 7^ 


from which the swelling stress oswo can be estimated. 


The results for Kx in Table 2.1 shall be used here to estimate at least the order of magni- 


tude of the residual stresses. 
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A4 Swelling stresses from crack terminating angles 


From the shielding stress intensity factors we estimated the residual stresses [A4.5]. For 
the tests by Schell et al. [A4.6] and Haranoh et al. [A4.7], crack extension is governed by 
Kiip=Kic. Fracture toughness of soda-lime glass is according to Wiederhorn [A4.8]: 
Ki-=0.75 MPavm and Poisson's ratio v=0.225. 


From [A4.5] it follows for the applied stress intensity factor of the crack-terminating 


angle @ (soda-lime glass v=0.225) 


1 
K.lob)=A RR (A4.1.7) 
(95) (9) 00606 9 1 8.54x 10 7 go 1 


A4.2 Discussion of stresses 


All residual stresses in Table A4.1 were found to be negative indicating that the wa- 
ter/silica reaction, the ion exchange, and chemical toughening are accompanied by vol- 
ume expansion. The absolute values of the compressive stresses resulting by the exten- 
sion of the Underwood procedure depend on the coefficient A that is not yet known. As 
can be concluded from Underwood [A4.4], the value of X decreases with increasing devi- 
ation from the reference crack. A rough estimation of X as a function of may be sug- 


gested based on soda-lime glass undergoing ion-exchange. 


In [A4.6] the stresses by ion exchange were concluded from strength measurements as 
Ores = —2425 MPa. A theoretical value of 2300 MPa was determined in [A4.9, A4.10] for 
the volume change by the ion-exchange fully transformed in stresses. As an average val- 
ue, Ores = —2.4 GPa may be used to compute the value of A. From the data of —AOres in 
Table A4.1 and the related terminating angles @, o in Table 2.1, the red circles in Fig. 
A4.2a were obtained. In this plot, the terminating angles ọ were normalized on @=90° 


representing the reference case (black circle). 


The solid curve represents a quadratic fit based on the results for the soda-lime glass with 


the “weight” for the 1-week result (in parentheses) described by 


A(g) x 1-9 c p ) (A4.2.1) 
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A4.2 Discussion of stresses 


In the case of swelling stresses in silica, the surface stresses at 07250?C can be computed 
from the hydroxyl concentration at the free surface as shown in [A4.11]. The swelling 


stress components in the surface plane are 


o =~ 31375MPa exp(0.00868 0) (A4.2.2) 
SW,X SW, 2 
x 5 + lexp(Q/RT) 


with A=32.3 and Q=10.75 kJ/mol, 7-04273*. 
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Fig. A4.2 Coefficient X of eq.(A4.1.3) as a function of the crack-terminating angle q; red circles: results for 
soda-lime glass, red line: related fitting curve for reduced weight of the data point at @=71° (in 
parentheses); blue circles: results for silica, b) swelling stress under saturation pressure from 
crack-terminating angles compared with results from disk experiments, and water concentrations 


by Zouine et al [A4. 12], full symbols: liquid water, open symbols: water vavour. 


Equation (A4.2.2) yields for 02250?C: 654, =—209 MPa. The results for silica predicted 
with eq.(A4.2.1), are introduced in Fig. A4.2a by the blue circles, showing compressive 
stresses in the range of -130 MPa to -170 MPa. 
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A4 Swelling stresses from crack terminating angles 


Applying eq.(A4.2.1) on the chemically toughened soda-lime glass gives a residual stress 
Of Ores=—1980 MPa. This result is in agreement with a surface stress of “more than —1000 
MP2" as had been reported by Haranoh et al. [A4.7]. 


Glass Treatment Ao: (MPa) | Gres (MPa) 
Soda-lime Chem. toughened -174 -1980 
Ion exchange | week -1009. «2400 
5 weeks -976. «2400 
6 weeks -968. «2400 
Silica 192h H>O-vapour —52.4 -128 
48h H20-liquid —62.1 -151 
192h H>O-liquid —69.7 -170 


Table A4.1 Estimates of residual stresses, Ores=Osw, in soda-lime glass and silica, predicted residual stresses 


for silica and chemically toughened soda-lime glass by use of eq.(A4.2.1). 


Figure A4.2b gives a comparison between the swelling stresses from crack terminating 
angles with data from disk bending experiments, [A4.2], and hydroxyl concentrations via 
eqs.(1.1.11) and (3.1.4). The results for liquid water and 192 h soaking time are about 
15% below the curve predicted by eq.(A4.2.2). For the shorter soaking time, the devia- 
tion is somewhat larger. Vapour soaking results in clearly deviations from the curve and 
show the same effect as found for the disk tests, in both cases an indication for the fact, 


that saturation was not reached. 


The silica/water reaction at silica surfaces exposed to water causes volume expansion in 
the water diffusion layer and as the consequence of mechanical boundary conditions 
compressive stresses. Evaluation of crack-terminating angles resulted in strongly nega- 
tive shielding stress intensity factors of about K-—2.5 MPa\m (see Table 2.1), a clear 
evidence for compressive stresses in the water-affected surface layer. Consequently, the 
existence of swelling stresses due to the silica/water reaction is proven. For silica, heat- 
treated in humid environments, compressive stresses in the order of about -130 MPa to - 
170 MPa are obtained for 250?C. 
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loading. In the present booklet we deal with results and interpretations of strength 
increase in the absence of applied stresses. In detail the following topics are discussed: 


- Estimation of swelling effects on crack-terminating angles. 
- Strength in absence and presence of humid environments. 
- Crack-tip shielding. 

- Anomalous subcritical crack growth. 
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